Я пытаюсь подобрать следующий блеск: разложение (пропорция) помета в зависимости от ковариат. Фиксированными ковариатами являются тип чая (категориальный с двумя уровнями), микробитат (категориальный с тремя уровнями), MAP (непрерывный) и MAT (непрерывный). Термины взаимодействия представляют собой MAT x MAP, MAT x microhabitat и MAP x microhabitat. Мы используем график, вложенный в сайт как случайный эффект, чтобы учесть зависимость наблюдений на графиках внутри сайтов.
MAP1 <-MAP / 100 </p>
M2<-glmer(Fraction_decomposed~MAT*MAP1*Treatment+Tea+(1|Site/Plot),family = binomial(link = "logit"))
Но я получил следующие предупреждающие сообщения:
Warning messages:
1: In eval(family$initialize, rho) :
non-integer #successes in a binomial glm!
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
unable to evaluate scaled gradient
3: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge: degenerate Hessian with 1 negative eigenvalues
введите описание изображения здесь
Я также попытался преобразовать переменную ответа и выполнил lmer, но я прочитал, что преобразования не рекомендуется
Fraction_decomposed1 <-lo git (Fraction_decomposed) </p>
M3 <-lmer (Fraction_decomposed1 ~ MAT * MAP1 * Лечение + Чай + (1 | Участок / Участок)) isSingular (M3) Резюме Anova (M3) (M3) </p>