Ich passe einige Linien an und habe das Gefühl, ich sage R genau, wie sie zu passen sind, aber ich habe das Gefühl, dass es etwas gibt (einen Faktor oder Effekt), von dem ich nicht weiß, dass es eine gute Passform verhindert.

Meine experimentelle Einheit ist "Handlung" wie in der Feldhandlung, was leider verwirrend ist.

Die Daten können gefunden werden:https://www.dropbox.com/s/a0tplyvs8lxu1d0/rootmeansv2.csv . mit

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

Ich habe einfach eine Masse über der Zeit und passe sie mit dem Modell an jede experimentelle Einheit an:

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

und mitplot (augPred(fit)) Ich bekomme diese Anpassungen für jede experimentelle Einheit ("Handlung"):

Was muss ich tun, um die Steigung zwischen den experimentellen Einheiten stärker variieren zu lassen? Das interessiert mich nicht aus statistischer Sicht, sondern aus prädiktiver Sicht - also kann alles im Modell manipuliert werden, um diese Linien in Bewegung zu versetzen.