Gemischte (auch als mehrstufige oder hierarchische) Modelle sind lineare Modelle, die sowohl feste als auch zufällige Effekte enthalten. Sie werden verwendet, um longitudinale oder verschachtelte Daten zu modellieren.
Kürzlich wurde mir gesagt, dass es nicht möglich ist, zeitvariable Kovariaten in longitudinale gemischte Modelle einzubeziehen, ohne eine Zeitverzögerung für diese Kovariaten einzuführen. Können Sie dies bestätigen / ablehnen? Haben Sie Hinweise zu dieser Situation? Ich schlage eine einfache Situation zur Klärung vor. Angenommen, ich habe quantitative Messungen (etwa 30 …
Ich versuche, mehrere Interaktionstests zwischen beiden lmund lmerwiederholten Messungen (2x2x2) zu reproduzieren. Der Grund, warum ich beide Methoden vergleichen möchte, ist, dass das GLM von SPSS für wiederholte Messungen genau die gleichen Ergebnisse liefert wie der lmhier vorgestellte Ansatz. Am Ende möchte ich SPSS mit R-lmer vergleichen. Bisher habe ich …
Ich versuche, die Zufallseffektvorhersagen aus einem linearen gemischten Modell von Hand zu berechnen, und unter Verwendung der von Wood in Generalized Additive Models bereitgestellten Notation : eine Einführung mit R (S. 294 / S. 307 von PDF), bin ich verwirrt über die einzelnen Parameter repräsentiert. Unten finden Sie eine Zusammenfassung …
Als Gegenstück zu diesem Beitrag habe ich daran gearbeitet, Daten mit kontinuierlichen Variablen zu simulieren und mich für korrelierte Abschnitte und Steigungen zu eignen. Obwohl es auf der Website und außerhalb der Website großartige Beiträge zu diesem Thema gibt , hatte ich Schwierigkeiten, ein Anfang-zu-Ende-Beispiel mit simulierten Daten zu finden, …
In der lmerFunktion in lme4in Rgibt es einen Aufruf zum Erstellen einer Modellmatrix von Zufallseffekten , wie hier auf den Seiten 7 bis 9 erläutert .ZZZ Die Berechnung von beinhaltet KhatriRao- und / oder Kronecker-Produkte aus zwei Matrizen, J_i und X_i . ZZZJiJiJ_iXiXiX_i Die Matrix JiJiJ_i ist ein Schluck: "Indikatormatrix …
Ich führe Post-hoc-Tests an einem linearen Mischeffektmodell in R( lme4Paket) durch. Ich verwende multcompPaket ( glht()Funktion), um die Post-Hoc-Tests durchzuführen. Mein experimenteller Entwurf besteht aus wiederholten Messungen mit einem zufälligen Blockeffekt. Die Modelle sind wie folgt spezifiziert: mymod <- lmer(variable ~ treatment * time + (1|block), data = mydata, REML …
Bei drei Variablen yund x, die positiv stetig sind und zdie kategorisch sind, habe ich zwei Kandidatenmodelle, die gegeben sind durch: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) und fit.fe <- lm( y ~ 1 + x ) Ich hoffe, …
Bevor ich diese Frage stellte, habe ich unsere Website durchsucht und viele ähnliche Fragen gefunden (wie hier , hier und hier ). Ich bin jedoch der Meinung, dass diese verwandten Fragen nicht gut beantwortet oder diskutiert wurden, und möchte diese Frage daher erneut stellen. Ich denke, es sollte eine große …
Ich analysiere die Ergebnisse eines Reaktionszeitversuchs in R. Ich führte eine ANOVA mit wiederholten Messungen durch (1 Faktor innerhalb des Subjekts mit 2 Stufen und 1 Faktor zwischen den Subjekten mit 2 Stufen). Ich habe ein ähnliches lineares gemischtes Modell ausgeführt und wollte die früheren Ergebnisse in Form einer ANOVA-Tabelle …
Im gemischten Modell nehmen wir an, dass die zufälligen Effekte (Parameter) Zufallsvariablen sind, die Normalverteilungen folgen. Es sieht der Bayes'schen Methode sehr ähnlich, bei der alle Parameter als zufällig angenommen werden. Ist das Zufallseffektmodell also ein Sonderfall der Bayes'schen Methode?
Auf R-sig-ME gibt es mehrere Themen zum Erhalten von Konfidenzintervallen für Vorhersagen mit lme4und nlmein R. Zum Beispiel hier und hier im Jahr 2010, einschließlich einiger Kommentare von Dougals Bates, einem der Autoren beider Pakete. Ich zögere, ihn wörtlich zu zitieren, aus Angst, dass sie aus dem Zusammenhang gerissen werden, …
Mir wurde immer beigebracht, dass zufällige Effekte nur die Varianz (Fehler) beeinflussen und dass feste Effekte nur den Mittelwert beeinflussen. Aber ich habe ein Beispiel gefunden, bei dem zufällige Effekte auch den Mittelwert beeinflussen - die Koeffizientenschätzung: require(nlme) set.seed(128) n <- 100 k <- 5 cat <- as.factor(rep(1:k, each = …
Ich arbeite an einem Datensatz. Nachdem ich einige Modellidentifikationstechniken angewendet hatte, kam ich mit einem ARIMA (0,2,1) -Modell heraus. Ich habe die detectIOFunktion im Paket TSAin R verwendet, um bei der 48. Beobachtung meines ursprünglichen Datensatzes einen innovativen Ausreißer (IO) zu erkennen . Wie kann ich diesen Ausreißer in mein …
Ich verstehe die Auswirkung eines kategorialen Zufallseffekts auf ein Modell mit gemischten Effekten darin, dass es eine teilweise Zusammenfassung der Beobachtungen nach Ebene des zufälligen Effekts durchführt, wobei effektiv angenommen wird, dass die Beobachtungen selbst nicht unabhängig sind, sondern nur ihre Teilpools. Nach meinem Verständnis überwiegen in einem solchen Modell …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.