Wir wissen, dass ein gepaarter t- Test nur ein Sonderfall von Einweg-ANOVA mit wiederholten Messungen (oder innerhalb des Subjekts) sowie eines linearen Mischeffektmodells ist, das mit der Funktion lme () des nlme-Pakets in R demonstriert werden kann Wie nachfolgend dargestellt. #response data from...
Eine Mitarbeiterin analysiert einige biologische Daten für ihre Dissertation mit böser Heteroskedastizität (Abbildung unten). Sie analysiert es mit einem gemischten Modell, hat aber immer noch Probleme mit den Residuen. Durch die Protokolltransformation der Antwortvariablen werden die Dinge...
Ich habe mehrere Imputationen verwendet, um eine Reihe vollständiger Datensätze zu erhalten. Ich habe bei jedem der vervollständigten Datensätze Bayes'sche Methoden verwendet, um die posterioren Verteilungen für einen Parameter zu erhalten (ein zufälliger Effekt). Wie kann ich die Ergebnisse für...
Ich versuche, lmeaus dem nlmePaket Ergebnisse von aovANOVAs für wiederholte Messungen zu replizieren . Ich habe dies für ein Einzelfaktor-Experiment mit wiederholten Messungen und für ein Zweifaktor-Experiment mit einem Zwischensubjektfaktor und einem Innensubjektfaktor getan, aber ich habe...
Sagen wir, wir müssen GLMMs mod1 <- glmer(y ~ x + A + (1|g), data = dat) mod2 <- glmer(y ~ x + B + (1|g), data = dat) Diese Modelle sind nicht im üblichen Sinne verschachtelt: a <- glmer(y ~ x + A + (1|g), data = dat) b <- glmer(y ~ x + A + B + (1|g), data = dat) Also können wir...
Ich versuche, ein Modell zu erstellen, für das ich eine Antwortvariable habe, die einen Anteil zwischen 0 und 1 hat. Dazu gehören einige 0s und 1s, aber auch viele Werte dazwischen. Ich denke darüber nach, eine Beta-Regression zu versuchen. Das Paket, das ich für R (betareg) gefunden habe, erlaubt...
Betrachten Sie die folgenden Daten aus einer wechselseitigen Betrachtung innerhalb des Probandendesigns: df <- "http://personality-project.org/r/datasets/R.appendix4.data" df <- read.table(df,header=T) head(df) Observation Subject Task Valence Recall 1 1 Jim Free Neg 8 2 2 Jim Free Neu 9 3 3...
Bitte beachten Sie diese Daten: dt.m <- structure(list(id = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12), occasion = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("g1", "g2"), class =...
Als Antwort auf diese Frage, ob mein Entwurf, in dem ich Teilnehmer nach dem Zufallsprinzip mit Bildern aus verschiedenen Kategorien präsentierte, ein Beispiel dafür war, dass ich eine ANOVA mit wiederholten Maßen verwenden sollte, erhielt ich die Antwort, dass ich stattdessen ein gemischtes Modell...
Angenommen, wir arbeiten an einem Zufallseffektmodell einiger Zählungsdaten im Zeitverlauf und möchten einige Trends überprüfen. Normalerweise würden Sie etwas tun wie: lmer(counts ~ dependent_variable + (1+t+I(t^2)|ID), family="poisson") eine quadratische Form für einschließen t. Ist es möglich,...
Grundsätzlich frage ich mich, wie unterschiedliche Kovarianzstrukturen erzwungen werden und wie die Werte in diesen Matrizen berechnet werden. Funktionen wie lme () erlauben es uns, die Struktur auszuwählen, die wir möchten, aber ich würde gerne wissen, wie sie geschätzt werden. Betrachten Sie das...
Es gibt eine Unterscheidung, die mich mit gemischten Modellen auslöst, und ich frage mich, ob ich Klarheit darüber bekommen könnte. Nehmen wir an, Sie haben ein gemischtes Modell der Zähldaten. Es gibt eine Variable, von der Sie wissen, dass Sie sie als festen Effekt haben möchten (A) und eine...
Kann ich Hilfe bei diesem vorläufigen (in Bearbeitung befindlichen) Versuch erhalten, mich auf die Äquivalente von ANOVA und REGRESSION zu beziehen? Ich habe versucht, die Konzepte, die Nomenklatur und die Syntax dieser beiden Methoden in Einklang zu bringen. Es gibt viele Posts auf dieser Site...
Wie sollen insbesondere die Standardfehler der Fixeffekte in einem linearen Mischeffektmodell (im häufigeren Sinne) berechnet werden? Ich bin zu der Annahme , dass die typischen Schätzungen ( ), wie sie in Laird und Ware [1982] vorgestellt wurden, SE's dazu geben werden werden in ihrer Größe...
Ich weiß, dass einer der Vorteile gemischter Modelle darin besteht, dass sie die Angabe einer Varianz-Kovarianz-Matrix für die Daten ermöglichen (zusammengesetzte Symmetrie, autoregressiv, unstrukturiert usw.). Die lmerFunktion in R ermöglicht jedoch keine einfache Angabe dieser Matrix. Weiß...
Hintergrund: Hinweis: Mein Datensatz und R-Code sind unter dem Text enthalten Ich möchte AIC verwenden, um zwei Modelle mit gemischten Effekten zu vergleichen, die mit dem lme4-Paket in R erstellt wurden. Jedes Modell hat einen festen und einen zufälligen Effekt. Der festgelegte Effekt...
Hallo, ich habe zwei Probleme, die wie natürliche Kandidaten für mehrstufige / gemischte Modelle klingen, die ich nie benutzt habe. Die einfachere Variante, die ich als Einführung versuchen möchte, sieht wie folgt aus: Die Daten sehen aus wie viele Zeilen des Formulars x y innergroup outergroup...
Ich verstehe, dass wir Modelle mit zufälligen Effekten (oder gemischten Effekten) verwenden, wenn wir glauben, dass einige Modellparameter über einen Gruppierungsfaktor zufällig variieren. Ich möchte ein Modell xanpassen, bei dem die Reaktion über einen Gruppierungsfaktor normalisiert und zentriert...
Ich habe mehrere Bücher durchgesehen (Raudenbush & Bryk, Snijders & Bosker, Gelman & Hill usw.) und mehrere Artikel (Gelman, Jusko, Primo & Jacobsmeier usw.), und ich habe meinen Kopf immer noch nicht richtig umwickelt Die Hauptunterschiede zwischen der Verwendung von gruppierten...
Was sind einige der Hauptprobleme bei der Verwendung linearer Modelle mit gemischten Effekten? Was sind die wichtigsten Dinge, auf die Sie achten müssen, um die Angemessenheit Ihres Modells zu beurteilen? Nach welchen wichtigen Dingen muss beim Vergleich von Modellen desselben Datensatzes gesucht...