library(lme4) library(car) library(bbmle) A <- read.table(file = "AREA_corr.txt", header = TRUE, dec = ".") A$Treatment = factor(A$Treatment) A$Day = as.numeric(A$Day) A$Label = factor(A$Label) A$small = factor(A$small) # Remove the individuals that remained small during the whole experiment aa = subset(A, small =="0") str(aa) a = na.omit(aa) str(a) ############ # ANALYSES ################################################################## ############ size.x2 = lmer((Area) ~ Treatment+strain+Day+ I(Day^2) +Treatment:strain + Treatment:Day+strain:Day+Treatment:I(Day^2)+strain:I(Day^2) + Treatment:strain:Day+Treatment:strain:I(Day^2)+ (Day|Label), data = a) #full model, random slope size.x3 = lmer((Area) ~ Treatment+strain+Day+ I(Day^2) +Treatment:strain + Treatment:Day+strain:Day+Treatment:I(Day^2)+strain:I(Day^2) + Treatment:strain:Day+ (Day|Label), data = a) #removing one 3-way interaction size.x4 = lmer((Area) ~ Treatment+strain+Day+ I(Day^2) +Treatment:strain +Treatment:Day+strain:Day+Treatment:I(Day^2) +strain:I(Day^2) + (Day|Label), data = a) #removing both 3-way AICctab(size.x2,size.x3,size.x4) #model selection Anova(size.x4,3) #To ensure that the pattern is not driven by the 20C treatment, we rerun the analysis without this treatment. # Without 20C a.no20<-subset(a,a$Treatment !=20) size.no20.x2 = lmer((Area) ~ Treatment+strain+Day+ I(Day^2) +Treatment:strain + Treatment:Day+strain:Day+Treatment:I(Day^2)+strain:I(Day^2) + Treatment:strain:Day+Treatment:strain:I(Day^2)+ (Day|Label), data = a.no20) #full model, random slope size.no20.x3 = lmer((Area) ~ Treatment+strain+Day+ I(Day^2) +Treatment:strain + Treatment:Day+strain:Day+Treatment:I(Day^2)+strain:I(Day^2) + Treatment:strain:Day+ (Day|Label), data = a.no20) #removing one 3-way interaction size.no20.x4 = lmer((Area) ~ Treatment+strain+Day+ I(Day^2) +Treatment:strain +Treatment:Day+strain:Day+Treatment:I(Day^2) +strain:I(Day^2) + (Day|Label), data = a.no20) #removing both 3-way AICctab(size.no20.x2,size.no20.x3,size.no20.x4) #model selection Anova(size.no20.x4,3)