data(sim.rasch)
# estimate Rasch model
mod1 <- tam.mml(resp=sim.rasch)
# standard error estimation
se1 <- tam.se( mod1 )
## > print(se1$xsi, digits=3 )
## xsi.index xsi.label est se
## 1 1 I1 -1.9590 0.0654
## 2 2 I2 -1.8570 0.0639
## 3 3 I3 -1.7444 0.0623
## 4 4 I4 -1.6407 0.0610
Run the code above in your browser using DataLab