Die Daten stammen von hier
%Vor% Hier passte ich ein nichtlineares Mixed-Effects-Modell mit nlme im nlme -Paket an. Das Modell hat 4 feste Effekte und 4 zufällige Effekte. Ich habe eine diagonale Struktur in der Varianz-Kovarianz-Matrix angegeben, und jedes patid bildet eine Gruppe.
Hier möchte ich das gleiche Modell mit dem lme4 -Paket anpassen. Aus der Dokumentation scheint, dass die formula für nlmer auch eine Gradientenkomponente haben muss, also habe ich zuerst die Funktion deriv verwendet. Wie kann ich den Rest der Parameter spezifizieren? Das
gibt 4 feste Effekte (im ersten Listenobjekt) und die entsprechenden 4 zufälligen Effekte (im zweiten Listenobjekt) an. Allerdings bin ich nicht ganz sicher, wie man eine diagonale Varianz-Kovarianz-Struktur spezifiziert und stelle sicher, dass die Beobachtungen nach patid gruppiert sind, wie ich es in random = list(patid = pdDiag(list(p1 ~ 1, b1 ~ 1, p2 ~ 1, b2 ~ 1))) mit nlme spezifiziert habe.
Die Standardmethode zum Festlegen der fixierten Effekte FE1 ... FE4 und der unabhängigen Zufallseffekte RE1 ... RE4 wird unten angezeigt
Das nlme -Paket hat eine etwas andere Syntax als lme4 package.
Das heißt, ich bin mir nicht sicher, ob ich die Nuancen Ihrer Frage vollständig verstehe, also ist es möglich, dass Ihre Situation bedeutet, dass dies etwas modifiziert werden muss. Wenn Sie Kommentare abgeben, überarbeite ich gerne meine Antwort nach Bedarf.
Tags und Links r lme4 mixed-models nlme