# NOT RUN {
if(require(TAM)){
set.seed(2524)
diff <- rnorm(10)
N <- 500
th <- rnorm(N)
probs <- 1 / (1 + exp(-outer(th, diff, "-")))
dat <- apply(probs, 2, function(p) as.numeric(p > runif(N)))
# fit the model
mod <- TAM::tam.mml(resp = dat, irtmodel = "1PL")
get_thresh(mod$xsi$xsi, itemtype = "1PL", item_ind = 1:10)
}
# }
Run the code above in your browser using DataLab