Я пытаюсь понять, почему результаты запуска модели вида ~inner|outer не сопоставимы с результатами запуска ~1 | inner ,~ 1 | outer в некоторых случаях
В примере simple , результаты сопоставимы , как и ожидалось
# Simple example
case <- data.frame(Dataset= c("a","b","c","d"),Cohort=c("c1","c1","c2","c3"), Tech= c("a1","a2","a1","a1"),Effect_size=c(-1.5,-3,1.5,3),Standard_error=c(.2,.4,.2,.4))
res1 = rma.mv(Effect_size, Standard_error^2, random = list(~ 1 | Dataset,~ 1 | Cohort), data=case)
res2=rma.mv(Effect_size, Standard_error^2, random = ~ Dataset | Cohort, data=case)
В более сложных случаях, однако [см. пример ниже] результаты двух моделей не сопоставимы На самом деле набор данных ~ | Когорт оценивает значение в 0,02, а список (~ 1 | Набор данных, ~ 1 | Когорт) оценивает значение в .2
case <- data.frame(Dataset= c("Dt1","Dt2", "Dt3", "Dt4","Dt5","Dt5","Dt6", "Dt7","Dt8","Dt9"),
Cohort=c("c1","c2","c3","c3","c3","c3","c3","c4","c4","c5"),
Effect_size=c(-0.002024454, -0.003915314, -0.042282757,-1.43826175,-0.045423574,-0.17682309,-21.72691245, -2.559727204,-0.091972279,-0.763332081),
Standard_error=c(0.15283972,0.117452325,0.262002289,0.555230971,0.708917912,0.682989908,2.704749864,1.40514335,0.735696048, 0.713557015))```