Включение корреляционной структуры блеска в функцию пакета lme4 - PullRequest
Новая
       15

Включение корреляционной структуры блеска в функцию пакета lme4

0 голосов
/ 18 мая 2018

Я использую функцию glmer пакета lme4 из R. Мой ответ - Бернулли, а мои данные продольные.

В пакете bild я уже нашел статью, использующую Mc1 или MC2 для корреляционных структур.

Можно ли как-нибудь использовать функцию glmer и иметь возможность включать структуру корреляции?Если это невозможно, не могли бы вы сказать мне, какое значение по умолчанию?

Код в R предназначен только для облегчения помощи. 

library(dplyr)
n = 300
xx<-c("r1","r2","r3","r4","r5")
xxx<-c("e1","e2","e3")
p=0.3
df1 <- data_frame(
  xx1 = runif(n, min = 0, max = 10),
  xx2 = runif(n, min = 0, max = 10),
  xx3 = runif(n, min = 0, max = 10),
 School = factor(sample(xxx, n,re=TRUE)),
 yx = as.factor(rbinom(n, size = 1, prob = p))
)
df1
library(lme4)
mm2 <- glmer(yx ~ xx1 + xx2 + xx3 +  (1 | School), data = df1, 
family = "binomial",control = glmerControl(calc.derivs = FALSE))
...