Ich versuche, Werte über die Zeit (Tage in x-Achse) für ein glmeres Modell vorherzusagen, das mit meinen Binomialdaten ausgeführt wurde. Total Alive und Total Dead sind Zähldaten. Dies ist mein Modell und die entsprechenden Schritte unten.

full.model.dredge<-glmer(cbind(Total.Alive,Total.Dead)~(CO2.Treatment+Lime.Treatment+Day)^3+(Day|Container)+(1|index),
                         data=Survival.data,family="binomial")

Wir haben eine Überdispersion berücksichtigt, wie Sie im Code sehen können (1: Index).

Anschließend verwenden wir den Befehl Dredge, um die am besten angepassten Modelle mit den Haupteffekten (CO2-Behandlung, Kalk-Behandlung, Tag) und den entsprechenden Wechselwirkungen zu bestimmen.

dredge.models<-dredge(full.model.dredge,trace=FALSE,rank="AICc")

Dann machte eine Arbeitsbereichsvariable für sie

my.dredge.models<-get.models(dredge.models)

Anschließend haben wir einen Modelldurchschnitt erstellt, um die Koeffizienten für die am besten geeigneten Modelle zu mitteln

silly<-model.avg(my.dredge.models,subset=delta<10)

Jetzt möchte ich ein Diagramm mit der Gesamtlebensdauer auf der Y-Achse und den Tagen auf der X-Achse und einer angepassten Linie in Abhängigkeit von der Ausgabe des Modells erstellen. Ich verstehe, dass dies schwierig ist, weil das Modell Total.Alive und Total.Dead verkettet hat (siehecbind(Total.Alive,Total.Dead) im Modell.

Wenn ich versuche, einen Vorhersagebefehl auszuführen, wird der Fehler angezeigt

# 9: In UseMethod("predict") :
#   no applicable method for 'predict' applied to an object of class "mer"