mixed-models

Extract prediction band z LME fit

Mam następujący model x <- rep(seq(0, 100, by=1), 10) y <- 15 + 2*rnorm(1010, 10, 4)*x + rnorm(1010, 20, 100) id &l ... owiedź populacji. Jak wyodrębnić i wykreślić 95% przedziały ufności / pasma Przewidywania z obiektu nlme dla całej populacji?

lme4:: raport lmer "fixed-effect model matrix is rank deficient", Czy potrzebuję poprawki i jak to zrobić?

Próbuję uruchomić model o mieszanych efektach, który przewiduje F2_difference z resztą kolumn jako predyktorami, ale dostaję ... iwić rzeczywiste uruchomienie modelu w lmer dla czytnika. Ale to nie jest zbyt duży problem. To jest nadal bardzo dobry post!