Estou ajustando algumas linhas e sinto que estou dizendo ao R exatamente como ajustá-las, mas sinto que há algo (algum fator ou efeito) que desconheço que está impedindo um bom ajuste.

Minha unidade experimental é "plot", como na plotagem de campo, que sinto muito é confusa.

Os dados podem ser encontrados:https://www.dropbox.com/s/a0tplyvs8lxu1d0/rootmeansv2.csv . com

df$plot.f<-as.factor(df$plot)
dfG<-groupedData(mass ~ year|plot.f, data=df)
dfG30<-dfG[dfG$depth == 30,]

Simplesmente, tenho massa ao longo do tempo e a encaixo em cada unidade experimental com o modelo:

fit <- lme(mass ~ year , random = ~ 1 | plot, data = df)

e complot (augPred(fit)) Recebo esses ajustes para cada unidade experimental ("plot"):

O que preciso fazer para permitir que a inclinação varie mais entre as unidades experimentais? Não estou interessado nisso de uma perspectiva estatística, mas de uma perspectiva preditiva - para que qualquer coisa no modelo possa ser manipulada para que essas linhas se movam.