Я пытаюсь подобрать логистическую модель случайного перехвата, используя функцию glmer из пакета lme4. К сожалению, я получаю следующие предупреждающие сообщения и явно неправильные результаты (для коэффициентов).
Warning messages:
1: In vcov.merMod(object, use.hessian = use.hessian) :
variance-covariance matrix computed from finite-difference Hessian is
not positive definite: falling back to var-cov estimated from RX
2: In vcov.merMod(object, correlation = correlation, sigm = sig) :
variance-covariance matrix computed from finite-difference Hessian is
not positive definite: falling back to var-cov estimated from RX
Проведя небольшое исследование, я обнаружил, что glmer создает объекты merMod. Опираясь на это: http://jaredknowles.com/journal/2014/5/17/mixed-effects-tutorial-2-fun-with-mermod-objects и воспроизведя раздел «Изучение внутреннего устройства объекта merMod», я получил следующие результаты. что касается моей модели:
### [1] "standardGeneric"
###attr(,"package")
###[1] "methods"
что явно отличается от
### [1] "lmerMod"
### attr(,"package")
### [1] "lme4"
как указано в учебнике.
Мой вопрос, почему мой объект не относится к классу merMod? Связаны ли приведенные выше предупреждения с этим и как я могу их исправить?
Это код, который я использовал для создания фрейма данных и запуска модели.
diagn00<- rep(0,240)
drug00<- rep(0,240)
time00<- c(rep(0,80),rep(1,80),rep(2,80))
response00<- c(rep(0,39),rep(1,41),rep(0,33),rep(1,47),rep(0,26),rep(1,54))
patients00<- rep(1:80,3)
test<- data.frame(patients00,diagn00,drug00,time00,response00)
diagn01<- rep(0,210)
drug01<- rep(1,210)
time01<- c(rep(0,70),rep(1,70),rep(2,70))
response01<- c(rep(0,33),rep(1,37),rep(0,15),rep(1,55),rep(0,2),rep(1,68))
patients01<- rep(81:150,3)
diagn10<- rep(1,300)
drug10<- rep(0,300)
time10<- c(rep(0,100),rep(1,100),rep(2,100))
response10<- c(rep(0,79),rep(1,21),rep(0,72),rep(1,28),rep(0,54),rep(1,46))
patients10<- rep(151:250,3)
diagn11<- rep(1,270)
drug11<- rep(1,270)
time11<- c(rep(0,90),rep(1,90),rep(2,90))
response11<- c(rep(0,74),rep(1,16),rep(0,45),rep(1,45),rep(0,15),rep(1,75))
patients11<- rep(251:340,3)
diagnosis<- c(diagn00,diagn01,diagn10,diagn11)
diagnosis<- as.factor(diagnosis)
id<- c(patients00,patients01,patients10,patients11)
id<- as.factor(id)
drug<- c(drug00,drug01,drug10,drug11)
drug<- as.factor(drug)
time<- c(time00,time01,time10,time11)
response<- c(response00,response01,response10,response11)
id<- c(patients00,patients01,patients10,patients11)
data<- data.frame(id, response, diagnosis, drug, time)
e<- order(data$id)
d<- data[e,]
library(lme4)
d2<- data.frame(d)
d2$response<- as.factor(d2$response)
d2$time<- as.factor(d2$time)
d2$id<- as.factor(d2$id)
t<- glmer(response ~ diagnosis + drug + time + time:drug +
(1 | id), family=binomial, data=d2)
информация о сеансе
e: [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C
[5] LC_TIME=English_United States.1252
прикрепленные базовые пакеты: [1] stats graphics grDevices utils наборы данных методы базы
другие прикрепленные пакеты: [1] glmmML_1.0 nlme_3.1-111 MASS_7.3-29 lme4_1.1-6 Rcpp_0.11.1 Matrix_1.1-3 MuMIn_1.10.0 [8] gee_4.13-18 geepack_1.1-6
загружено через пространство имен (и не прикреплено): [1] grid_3.0.2 решетка_0.20-29 minqa_1.2.3 RcppEigen_0.3.2.1.1 splines_3.0.2
[6] tools_3.0.2
информация о пакете
packageVersion("lme4")
[1] ‘1.1.6’
control
, я получу ошибку (сбой с неиспользуемым аргументом (maxIter=500)). Если я прокомментирую это, модель работает нормально и дает, по-видимому, разумные результаты. РезультатыsessionInfo()
/packageVersion("lme4")
... ? - person Ben Bolker   schedule 10.07.2014