R: извлечь матрицу дизайна из объекта MixedClass, используя пакет mirt
Я использую mixedmirt функция из пакета R mirt оценивать модели смешанных эффектов с ковариатами элементов как фиксированные эффекты (по существу, модель IRT с объяснением элементов или LLTM). Я хочу знать, хранится ли матрица дизайна (матрица модели) где-либо в объекте MixedClass, и как ее извлечь.
Пример кода из смешанной виньетки
##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 В любом случае для выполнения оценки необходимо сгенерировать матрицу проектирования, но хранится ли она где-нибудь? И как мне его извлечь?