Извлечение случайных эффектов уровня 2 из вмененного набора данных с помощью MItools - PullRequest
Новая
       125

Извлечение случайных эффектов уровня 2 из вмененного набора данных с помощью MItools

0 голосов
/ 12 июля 2020

Я пытаюсь извлечь случайные эффекты 3-го уровня из 3-уровневой иерархической регрессии, которую я выполняю с lmer. Проблема в том, что мой набор данных исчислен (правдоподобные значения), и я пытаюсь извлечь ТОЛЬКО случайные эффекты L3 для объединения с MItools. Я могу извлечь их на обычной модели lmer только с одним PV 

m0 <- lmer(PV4MATH ~ 1 + (1 | SCHOOLID) + (1 | CNT), data = pisa, REML = F)

ranef (m0) $ CNT

Но у меня не получается извлечь именно те из модели, которую я подхожу с mitools на всех PV, а затем попытайтесь объединить их с MIcombine 

model0<-with(mis, lmer(MATH ~ 1 + (1 | SCHOOLID) + (1 | CNT), REML = F))

Я хотел бы извлечь только случайные эффекты L 2 (= средние значения для всех стран при контроле переменных индивидуального уровня, я знаю, что это пустая модель). Вот как я могу извлечь все ranefs, но тогда у меня проблема с MIcombine. 

betas<-MIextract(model0, , fun = ranef)

Я пробовал разные комбинации $ CNT, но это не сработало. Есть у кого-нибудь идеи? Я получаю сообщение об ошибке, сообщающее мне, что это нечисловая операция c для бинарного оператора (не уверен, что это правильный перевод). У меня нет проблем с фиксированными эффектами.

1 Ответ

0 голосов
/ 12 июля 2020

Нашел решение ... На всякий случай, если кому-то это тоже понадобится: 

betas2 <- list()

for ( i in 1:5){ betas2[[i]] <- betas[[i]]$CNT

}

Извлеките все коэффициенты в списке «вручную», а затем вы можете использовать их. Также работает с отклонениями.

...