Разница в lme и lme4 в R? Отличие от STATA тоже? - PullRequest
Новая
       114

Разница в lme и lme4 в R? Отличие от STATA тоже?

1 голос
/ 30 сентября 2019

Я занимаюсь моделированием смешанных эффектов (используя данные в два момента времени), и мне интересно, почему существуют различия в пакетах R, а также между программами. Вот воспроизводимый пример: 

  USArrests$ID <- as.factor(matrix(1:50))
  USArrests <- USArrests[,-c(2,3,4)]
  USArrests$Murder2 <- USArrests$Murder + 5
  library(reshape)
  USArrests <- melt(USArrests, id.vars=c("ID"))[,-2]
  USArrests$time <- c(matrix(rep(0,50)), matrix(50:99))
  m1 <- lme(value ~ time, data=USArrests, random = ~1|ID)
  print(summary(m1))
  m2 <- lmer(value ~ time +(1|ID), data = USArrests, REML=F)
  print(summary(m2))

Затем я импортирую это в Stata и использую код: 

mixed value time || id:time

Даже в R модели немного отличаются, хотя по сравнению со STATA,есть большая разницаПочему существует эта разница? Это в алгоритмах, используемых для оценки максимальной вероятности?

Спасибо!

...