Я запускаю GLMM, используя glmer()
в R:
glmer(survive ~ fyear + site + fyear * site.x + (1|fyear),
family = binomial(link = logexp(shaffer.sub$exposure)),
data = shaffer.sub)
где выживание равно 0 или 1 в зависимости от того, было ли гнездо успешным или нет. Здесь вы можете увидеть, как выглядят данные:
structure(list(id = structure(1:7, .Label = c("1", "2", "3",
"4", "5", "6", "7"), class = "factor"), year.x = structure(c(1L,
1L, 2L, 3L, 3L, 3L, 3L), .Label = c("1994", "1995", "1999"), class = "factor"),
survive = structure(c(1L, 2L, 2L, 2L, 2L, 2L, 1L), .Label = c("0",
"1"), class = "factor"), fyear = structure(c(1L, 1L, 2L,
3L, 3L, 3L, 3L), .Label = c("1994", "1995", "1999"), class = "factor"),
site.x = structure(c(1L, 2L, 1L, 1L, 1L, 2L, 1L), .Label = c("N",
"S"), class = "factor")), .Names = c("id", "year.x", "survive",
"fyear", "site.x"), row.names = c(NA, -7L), class = "data.frame")
но я получаю это предупреждающее сообщение:
*Warning messages:
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 0.0299425 (tol = 0.001, component 12)
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model is nearly unidentifiable: large eigenvalue ratio
- Rescale variables?*
Мне сказали, что я не должен использовать тот же случайный фактор в качестве фиксированного эффекта для той же модели.
В конце я хотел бы получить результат, в котором я могу видеть год, сайт и год взаимодействия: эффекты сайта. Как в таблице ANOVA (возможно ли это? Я пытался использовать summary(aov(model))
, но это не работает; anova(model)
тоже не работает.
Я получаю эту ошибку для команды aov()
:
*Error in summary`(aov(syearXsite))` :
error in evaluating the argument 'object' in selecting a method for function 'summary': Error in if (fixed.only) { : argument is not interpretable as logical*
Как я могу увидеть влияние этих переменных на выживаемость?