R: извлечь матрицу дизайна из объекта MixedClass, используя пакет mirt - PullRequest
Новая
       37

R: извлечь матрицу дизайна из объекта MixedClass, используя пакет mirt

0 голосов
/ 25 февраля 2019

Я использую функцию 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 должен в любом случае генерировать матрицу проектирования, чтобы выполнить оценку, но хранится ли она где-нибудь?И как мне его извлечь?

...