Ich versuche ein lineares gemischtes Modell zu erzeugen. Der R-Code lautet wie folgt. lme (Average.payoff ~ Spiel + Typ + Andere.Typ + Spiel: Typ + Spiel: Andere.Typ + Typ: Andere.Typ, zufällig = ~ 1 | Subjekte, Methode = "REML", Daten = Subjektsm1) -> lme1 Der Antwortterm Average.payoff ist...
Angenommen, ich habe eine Stichprobe von Häufigkeiten von 4 möglichen Ereignissen: Event1 - 5 E2 - 1 E3 - 0 E4 - 12 und ich habe die erwarteten Wahrscheinlichkeiten, dass meine Ereignisse eintreten: p1 - 0.2 p2 - 0.1 p3 - 0.1 p4 - 0.6 Mit der Summe der beobachteten Häufigkeiten meiner vier...
Intro Ich habe Teilnehmer, die unter zwei Bedingungen wiederholt kontaminierte Oberflächen mit E. coli berühren ( A = Handschuhe tragen, B = keine Handschuhe). Ich möchte wissen, ob es einen Unterschied zwischen der Menge an Bakterien auf ihren Fingerspitzen mit und ohne Handschuhe gibt, aber auch...
Das nicht so seltene Auftreten bei komplexen maximal gemischten Modellen (Schätzung aller möglichen zufälligen Effekte für bestimmte Daten und Modelle) ist eine perfekte (+1 oder -1) oder nahezu perfekte Korrelation zwischen einigen zufälligen Effekten. Betrachten wir zum Zweck der Diskussion das...
Ich passe ein glmerModell in das lme4R-Paket. Ich suche nach einer Anova-Tabelle mit dem darin gezeigten p-Wert, kann aber kein passendes Paket finden. Ist es möglich, es in R zu tun? Das Modell, das ich anpasse, hat folgende Form: model1<-glmer(dmn~period*teethTreated+(1|fullName),...
library(lme4) out <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), data = cbpp, family = binomial, contrasts = list(period = "contr.sum")) summary(out) Fixed effects: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.32337 0.22129 -10.499 < 2e-16 *** period1...
Ich habe einige Probleme, gleichwertige Ergebnisse zwischen einem aovModell zwischen wiederholten Messungen und einem lmergemischten Modell zu erhalten. Meine Daten und mein Skript sehen wie folgt
In einem logistischen verallgemeinerten linearen gemischten Modell (Familie = Binomial) weiß ich nicht, wie die Varianz der zufälligen Effekte zu interpretieren ist: Random effects: Groups Name Variance Std.Dev. HOSPITAL (Intercept) 0.4295 0.6554 Number of obs: 2275, groups: HOSPITAL, 14 Wie...
Ich simuliere Bernoulli-Versuche mit einem zufälligen zwischen Gruppen und passe dann das entsprechende Modell mit an das Paket:logitθ ∼ N.( logitθ0, 12)logitθ∼N(logitθ0,12)\text{logit}\, \theta \sim {\cal N}(\text{logit}\, \theta_0, 1^2)lme4 library(lme4) library(data.table) I <- 30 # number...
Ich verwende lmer in R, um die Auswirkung von condition ( cond) auf ein Ergebnis zu überprüfen . Hier sind einige aus Daten, wobei s der Objektidentifizierer und a, bund csind Bedingungen. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30), a = rnorm(30, -2,...
Ich hoffe, jemand kann bei einer meiner Meinung nach relativ einfachen Frage helfen, und ich glaube, ich kenne die Antwort, aber ohne Bestätigung ist sie zu etwas geworden, dessen ich mir einfach nicht sicher sein kann. Ich habe einige Zähldaten als Antwortvariable und möchte messen, wie sich diese...
Ich habe einige Fragen zur Spezifikation und Interpretation von GLMMs. 3 Fragen sind definitiv statistisch und 2 beziehen sich genauer auf R. Ich poste hier, weil ich letztendlich denke, dass das Problem die Interpretation der GLMM-Ergebnisse ist. Ich versuche gerade, ein GLMM zu montieren. Ich...
Ich passe ein Modell mit gemischten Effekten mit einem Spline-Term in einer Anwendung an, bei der der zeitliche Trend bekanntermaßen kurvenlinear ist. Ich möchte jedoch beurteilen, ob der kurvenlineare Trend aufgrund einer individuellen Abweichung von der Linearität auftritt oder ob es sich um...
Ich habe gelesen, dass die Verteilung für den Zufallseffektschätzer in lme4 stark verzerrt ist und aus diesem Grund keine Standardfehler gemeldet werden. Ich frage mich, ob jemand eine Referenz dafür liefern kann. Ich habe Zugang zu dem Buch von Bates und Pinherio, aber nicht zu Raudenbush und...
Angenommen, A ist ein fester Faktor und B ist ein Zufallsfaktor. A und B sind gekreuzt. Ein weiterer Zufallsfaktor C ist in B verschachtelt. Wie sollte ich unter Berücksichtigung der AB- und AC-Wechselwirkungen das Modell mithilfe der lme-Funktion im nlme-Paket anpassen? Vielen
Ich habe 2 Faktoren Aund B(5 × 3) und eine Kovariate Xin einem innerhalb des Subjekts liegenden Design. So spezifiziere ich mein Gesamtmodell: lme.out = lme(y~ A*B*X, random=~1|Subject, data=mydata) Meine Interpretation ist, dass ich ein Diagramm betrachte y~x, in dem sich die Steigung aufgrund...
Ich habe ein Problem mit der Interpretation von 2-Wege- und 3-Wege-Interaktionen in lmer. Mein DV ist die Höhe, die eine kontinuierliche Variable ist. Alle IVs sind kategoriale Variablen. Der erste Faktor ist ein Tier, entweder eine Ratte oder ein Löwe. Der zweite Faktor ist das Geschlecht,...
Ich bin gespannt, wie das lmerTest-Paket in R, insbesondere die "rand" -Funktion, Tests mit zufälligen Effekten verarbeitet. Betrachten Sie das Beispiel aus dem lmerTest-PDF auf CRAN, das den integrierten Datensatz "Karotten" verwendet: #import lme4 package and lmerTest package library(lmerTest)...
Ich verstehe das Problem der Bestimmung der Freiheitsgrade in Mehrebenenmodellen; Daher ist die Entscheidung von Doug Bates et al. p-Werte nicht als Teil des lme4-Pakets in R zu melden. Ganz zu schweigen von der Fülle von Problemen mit und der unangemessenen Konzentration auf p-Werte im...
Was ist der Unterschied zwischen (1|DNA.concentration/mouse.id)und (DNA.concentration|mouse.id)? Was bedeuten die Symbole |und /die Bedeutung der Syntax für den Zufallseffekt?