Question

我正在使用R包mixedmirt中的mirt函数来估计带有项目协变量作为固定效果的混合效果模型（本质上是不言自明的IRT模型或LLTM）。我想知道设计矩阵（模型矩阵）是否存储在MixedClass对象的任何位置，以及如何提取它。

mixedmirt小插图的代码示例

##LLTM, and 2PL version of LLTM
data(SAT12)
data <- key2binary(SAT12,
                   key = c(1,4,5,2,3,1,2,1,3,1,2,4,2,1,5,3,4,4,1,4,3,3,4,1,3,5,1,3,1,5,4,5))
model <- 'Theta = 1-32'

itemdesign <- data.frame(itemorder = factor(c(rep('easier', 16), rep('harder', 16))))

LLTM <- mixedmirt(data, model = model, fixed = ~ 0 + itemorder, itemdesign = itemdesign,
   SE = TRUE) # SE argument ensures that the information matrix is computed accurately
summary(LLTM)

mirt不具有允许您从LLTM自动提取加权项目难度的功能，但是您可以通过提取系数并将其内积与设计矩阵相乘来手动计算它们（例如，通过使用原始项目上的model.matrix函数协变量）。

coef(LLTM)[[1]][1,1:2]%*%t(model.matrix(~itemorder+0, itemdesign))

但是，我想做一个通用函数，这样我就可以使用许多不同的模型来进行此操作，因此model.matrix方法对于不止一个解释性预测变量来说是不切实际的（仅需考虑单个因素））。如果您可以直接从MixedClass对象中提取矩阵，这将容易得多，例如您可以使用lme4命令对getME(LLTM,"X")包进行处理。

mixedmirt必须始终生成设计矩阵以进行估算，但是将其存储在任何地方吗？以及如何提取它？

R：使用包mirt

mixedmirt小插图的代码示例

0 个答案: