Non Linear, Regression, R, Modelling.
be honest. To be expected, the variance was not really at the start. Let's add a random slope as well. gamm_fit2.1<-mgcv::gamm(D9std~s(time), data=data_grouped, random=list(Big=~time)) anova(gamm_fit$lme, gamm_fit2$lme, gamm_fit2.1$lme) Now, surprise surpise, it does not seam that the random slope