## Data from 109th US Senate
data(s109)
## Convert data and make starts/priors for estimation
rc <- convertRC(s109)
p <- makePriors(rc$n, rc$m, 1)
s <- getStarts(rc$n, rc$m, 1)
## Conduct estimates
lout <- binIRT(.rc = rc,
.starts = s,
.priors = p,
.control = {
list(threads = 1,
verbose = FALSE,
thresh = 1e-6
)
}
)
## Look at first 10 ideal point estimates
lout$means$x[1:10]
Run the code above in your browser using DataLab