Non Linear, Regression, R, Modelling.
much use. Let's incorporate the hierarchy in the model by adding a random intercept. gamm_fit2<-mgcv::gamm(D9std~s(time), data=data_grouped, random=list(Big=~1)) plot(gamm_fit2$gam,pages=1) big <- data_grouped$Big rm(big) par(mfrow = c(2, 2)) mgcv::gam.check(gamm_fit2$gam) Nothing really changed to