Ich mache eine Simulationsstudie, die Bootstrapping-Schätzungen erfordert, die aus einem verallgemeinerten linearen gemischten Modell erhalten wurden (tatsächlich das Produkt von zwei Schätzungen für feste Effekte, eine von einem GLMM und eine von einem LMM). Um die Studie gut durchzuführen, wären...
Ich habe Google und diese Website durchsucht und bin immer noch verwirrt über die lmer-Funktion in der lme4-Bibliothek. Ich habe einige Daten von verschiedenen psychiatrischen Stationen gesammelt, die eine mehrstufige Struktur haben. Zur Vereinfachung werde ich zwei Variablen der Ebene 2 und zwei...
Stellen Sie sich ein Experiment mit mehreren menschlichen Teilnehmern vor, die jeweils mehrfach unter zwei Bedingungen gemessen wurden. Ein Modell mit gemischten Effekten kann (unter Verwendung der lme4- Syntax) wie folgt formuliert werden : fit = lmer( formula = measure ~ (1|participant) +...
Ich habe zwei logistische Regressionsmodelle in R erstellt mit glm(). Beide verwenden dieselben Variablen, wurden jedoch mit unterschiedlichen Teilmengen einer Matrix erstellt. Gibt es eine einfache Möglichkeit, ein Durchschnittsmodell zu erhalten, das die Mittelwerte der Koeffizienten angibt, und...
Ich untersuche das Zusammenspiel zweier Variablen ( und ). Zwischen diesen Variablen besteht eine große lineare Korrelation mit . Aus der Natur des Problems kann ich nichts über die Ursache sagen (ob verursacht oder umgekehrt). Ich möchte Abweichungen von der Regressionslinie untersuchen, um...
Wenn die Kontingenztabelle Nullen enthält und wir verschachtelte Poisson / Loglinear-Modelle (unter Verwendung der R- glmFunktion) für einen Likelihood-Ratio-Test anpassen, müssen wir die Daten vor dem Anpassen der glm-Modelle anpassen (z. B. 1/2 zu allen hinzufügen) die zählt)? Natürlich können...
Warum verwenden die Standardeinstellungen in R qqplot(linear model)die standardisierten Residuen auf der y-Achse? Warum verwendet R nicht die "regulären"
Ich wurde gebeten, einige Daten aus einer klinischen Studie zu analysieren, in der zwei Methoden zur Messung des Blutdrucks untersucht wurden. Ich habe Daten von 50 Probanden mit jeweils zwischen 2 und 57 Messungen mit jeder Methode. Ich frage mich, wie ich am besten vorgehen soll....
Angenommen, in einer Studie mit 15 Probanden wird die Antwortvariable (res) mit zwei erklärenden Variablen modelliert, wobei eine (Ebene) mit 5 Ebenen kategorisch ist und die andere (Antwortzeit: RT) kontinuierlich ist. Mit lmer im lme4-Paket von R habe ich: fm1 <- lmer(res ~ level * RT +...
Bei folgendem Versuchsaufbau: Von einem Probanden werden mehrere Proben entnommen und jede Probe wird auf verschiedene Arten behandelt (einschließlich einer Kontrollbehandlung). Was hauptsächlich interessant ist, ist der Unterschied zwischen der Kontrolle und jeder Behandlung. Ich kann mir zwei...
Diese Frage mag dumm klingen, aber ... ist es richtig, dass zufällige Effekte nur für kategoriale Variablen gelten können (wie individuelle ID, Populations-ID, ...), z. B. ist eine kategoriale Variable:xixix_i ~ β x iyiyiy_iβxiβxi\beta_{x_i} ~ N o r m ( μ , δ 2
Wikipedia schlägt vor, dass eine Möglichkeit, die Zuverlässigkeit zwischen Bewertern zu untersuchen, darin besteht, ein Zufallseffektmodell zur Berechnung der Korrelation zwischen Klassen zu verwenden . Das Beispiel der Intraclass-Korrelation spricht vom Betrachten
Ich habe folgende Art von Daten. Ich habe 10 Personen bewertet, die jeweils 10 Mal wiederholt wurden. Ich habe eine 10x10 Beziehungsmatrix (Beziehung zwischen allen Kombinationen der Individuen). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10,...
Mein Kollege und ich passen eine Reihe von linearen und nichtlinearen Mischeffektmodellen in R an. Wir werden gebeten, eine Kreuzvalidierung der angepassten Modelle durchzuführen, damit überprüft werden kann, ob die beobachteten Effekte relativ verallgemeinerbar sind. Dies ist normalerweise eine...
Wie spezifiziere ich in R ein früheres Modell ohne globalen festen Effekt? Zum Beispiel, wenn ich so etwas sage lmer(y ~ (1 | group) + (0 + x | group), data = my_df) das angepasste Modell wird sein yi j= a + αich+ βichxi jyij=a+αi+βixijy_{ij} = a + \alpha_{i} + \beta_{i} x_{ij} Wie passe ich...
Ich schaue, ob Überfluss mit Größe zusammenhängt. Die Größe ist (natürlich) kontinuierlich, jedoch wird die Häufigkeit auf einer solchen Skala aufgezeichnet, dass A = 0-10 B = 11-25 C = 26-50 D = 51-100 E = 101-250 F = 251-500 G = 501-1000 H = 1001-2500 I = 2501-5000 J = 5001-10,000 etc... A...
Angenommen, ich generiere die Wahrscheinlichkeit eines Ergebnisses basierend auf einem bestimmten Faktor und zeichne die Kurve dieses Ergebnisses. Gibt es eine Möglichkeit, die Gleichung für diese Kurve aus R zu extrahieren? > mod = glm(winner~our_bid, data=mydat,
Was ist der Unterschied zwischen einem GLM-Modell (logistische Regression) mit einer binären Antwortvariablen, die Subjekt und Zeit als Kovariaten enthält, und dem analogen GEE-Modell, das die Korrelation zwischen Messungen zu mehreren Zeitpunkten berücksichtigt? Mein GLM sieht aus...
Ich habe meine eigene, leicht verbesserte Version des Termplots erstellt, die ich in diesem Beispiel verwende. Sie finden sie hier . Ich habe zuvor auf SO gepostet, aber je mehr ich darüber nachdenke, desto mehr hängt dies wahrscheinlich mit der Interpretation des Cox Proportional Hazards-Modells...
Wie würde ich die Ausgabe einer Hauptkomponentenanalyse (PCA) in einem verallgemeinerten linearen Modell (GLM) verwenden, vorausgesetzt, die PCA wird für die Variablenauswahl für das GLM verwendet? Erläuterung: Ich möchte PCA verwenden, um die Verwendung korrelierter Variablen im GLM zu vermeiden....