我正在使用mirt
包来制作2PL
模型。我试图为每个估算参数获取standard error
。
使用mirt(mydata, 1, itemtype = '2PL',SE=TRUE)
后,我得到以下结果,我感到困惑,哪些信息显示标准错误。
仅供参考,mydata
是稀疏和二进制的。另外,我已经阅读了this page,但我没有得到答案。
第59项的结果样本
$`59`
a1 d g u
par 36.859 -12.266 0 1
CI_2.5 -905.576 -445.648 NA NA
CI_97.5 979.295 421.117 NA NA
哪些信息显示标准错误? a1
,d
,g
或u
。
欣赏!
答案 0 :(得分:1)
直接回答您的问题,该表中的任何信息都不是标准错误。
要从适合的mirt对象获取标准错误,您需要使用以下内容:
mydata.mirt <- mirt(mydata, model=1, itemtype="2PL", SE=TRUE)
mydata.coefficients <- coef(mydata.mirt, printSE=TRUE)
这应该给你一个类似的列表输出到你在问题中粘贴的内容,但是你应该获得每个参数估计的SE,而不是置信区间。
根据您发布的表格,我预计这些SE会非常高。可能值得检查数据和模型以确保没有出错,因为参数估计值令人惊讶且置信区间很大。