Ich habe eine zufällige logistische Intercept-Regression (aufgrund wiederholter Messungen) und möchte einige Diagnosen durchführen, insbesondere in Bezug auf Ausreißer und einflussreiche Beobachtungen. Ich habe mir die Residuen angesehen, um festzustellen, ob es Beobachtungen gibt, die auffallen....
Bei der Suche nach Informationen über das Randmodell und das Zufallseffektmodell und wie man zwischen ihnen wählt, habe ich einige Informationen gefunden, aber es war mehr oder weniger eine mathematisch abstrakte Erklärung (wie zum Beispiel hier: https: //stats.stackexchange .com / a / 68753/38080...
Ich suche nach Ratschlägen zur Analyse komplexer Umfragedaten mit Mehrebenenmodellen in R. Ich habe das surveyPaket verwendet, um ungleiche Auswahlwahrscheinlichkeiten in Einstufenmodellen zu gewichten, aber dieses Paket verfügt nicht über Funktionen für die Mehrebenenmodellierung. Das lme4Paket...
Betrachten Sie ein lineares Modell mit zufälligen Abschnitten. Dies entspricht einer linearen GEE-Regression mit einer austauschbaren Arbeitskorrelationsmatrix. Angenommen , die Prädiktoren sind und x 3 und die Koeffizienten für diese Prädiktoren sind β 1 , β 2 und β 3 . Was ist die Interpretation...
Ich habe einige Daten mithilfe der linearen Mischeffektmodellierung in R analysiert. Ich habe vor, ein Poster mit den Ergebnissen zu erstellen, und ich habe mich nur gefragt, ob jemand, der Erfahrung mit Mischeffektmodellen hat, vorschlagen könnte, welche Diagramme zur Veranschaulichung der...
Ich beende eine Analyse eines großen Datensatzes. Ich möchte das im ersten Teil der Arbeit verwendete lineare Modell verwenden und es mithilfe eines linearen gemischten Modells (LME) neu anpassen. Die LME wäre sehr ähnlich, mit der Ausnahme, dass eine der im Modell verwendeten Variablen als...
Ich hoffe aufrichtig, dass ich diese Frage so formuliert habe, dass sie definitiv beantwortet werden kann - wenn nicht, lassen Sie es mich bitte wissen und ich werde es erneut versuchen! Ich sollte auch beachten, dass ich R für diese Analysen verwenden werde. Ich habe mehrere Maßnahmen, von plant...
Mathematisch gesehen wird oft gesehen, dass Ausdrücke und Algorithmen für die Erwartungsmaximierung (Expectation Maximization, EM) für gemischte Modelle oft einfacher sind, aber es scheint, dass fast alles (wenn nicht alles), was mit EM gelöst werden kann, auch mit MLE gelöst werden kann (z. B....
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...
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...
Gibt es einen intuitiven Grund dafür, dass zufällige Effekte im allgemeinen linearen gemischten Modell auf ihren erwarteten Wert geschrumpft
Datenstruktur > str(data) 'data.frame': 6138 obs. of 10 variables: $ RT : int 484 391 422 516 563 531 406 500 516 578 ... $ ASCORE : num 5.1 4 3.8 2.6 2.7 6.5 4.9 2.9 2.6 7.2 ... $ HSCORE : num 6 2.1 7.9 1 6.9 8.9 8.2 3.6 1.7 8.6 ... $ MVMNT : Factor w/ 2 levels "_Withd","Appr": 2 2 1 1 2 1 2 1...
Ich habe ein GLMM der Form: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Wenn ich benutze drop1(model, test="Chi"), erhalte ich andere Ergebnisse als wenn ich Anova(model, type="III")aus dem Autopaket oder benutze summary(model). Diese...
Ich versuche, ein multivariates (dh Mehrfachantwort-) Mischmodell einzubauen R. Abgesehen von den ASReml-rund SabreR-Paketen (für die externe Software erforderlich ist) scheint dies nur in möglich zu sein MCMCglmm. In dem Papier , das die begleitet MCMCglmmPaket (pp.6) beschreibt Jarrod Hadfield...
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, diese Modelle vergleichen zu können, um festzustellen, welches...
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...
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...
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...
Was ist die Standardvarianz-Kovarianz-Struktur für zufällige Effekte im glmeroder lmerim lme4Paket? Wie spezifiziert man eine andere Varianz-Kovarianz-Struktur für zufällige Effekte im Code? Ich konnte diesbezüglich keine Informationen in der lme4Dokumentation
Ich analysiere einige Verhaltensdaten mit lme4in R, hauptsächlich nach Bodo Winters hervorragenden Tutorials , aber ich verstehe nicht, ob ich mit Interaktionen richtig umgehe. Schlimmer noch, niemand anderes, der an dieser Forschung beteiligt ist, verwendet gemischte Modelle. Ich bin also ein...