Die Frage zum Lebenslauf Ich versuche, (eine) detaillierte und präzise mathematische Darstellung (en) eines gemischten Effektmodells zu geben. Ich verwende das lme4Paket in R. Was ist die richtige mathematische Darstellung für mein Modell? Die Daten, die wissenschaftliche Frage und der...
Ich experimentiere mit dem Algorithmus der Gradientenverstärkungsmaschine über das caretPaket in R. Unter Verwendung eines kleinen Datensatzes für Hochschulzulassungen habe ich den folgenden Code ausgeführt: library(caret) ### Load admissions dataset. ### mydata <-
Bei der Angabe eines Modells mit gekreuzten gemischten Effekten versuche ich, Wechselwirkungen einzubeziehen. Es wird jedoch folgende Fehlermeldung angezeigt: Error in lme.formula(rate ~ nozzle, random = ~nozzle | operator, data = Flow) : nlminb problem, convergence error code = 1 message =...
Kann mir bitte jemand sagen, wie R den Bruchpunkt in einem stückweisen linearen Modell (als fester oder zufälliger Parameter) abschätzen soll, wenn ich auch andere zufällige Effekte abschätzen muss? Im Folgenden ist ein Spielzeugbeispiel aufgeführt, das eine Hockeyschläger- /...
Ich überprüfe gerade einige Arbeiten und bin auf Folgendes gestoßen, was mir falsch erscheint. Mit lmer werden zwei gemischte Modelle (in R) montiert. Die Modelle sind nicht verschachtelt und werden durch Likelihood-Ratio-Tests verglichen. Kurz gesagt, hier ist ein reproduzierbares Beispiel für...
Ich habe einen Datensatz, in dem die Variable, die ich als Zufallseffekt verwenden möchte, für einige Ebenen nur eine einzige Beobachtung enthält. Aufgrund der Antworten auf vorherige Fragen habe ich festgestellt, dass dies im Prinzip in Ordnung sein kann. Kann ich ein gemischtes Modell mit...
Angenommen, ich habe eine Antwortvariable yijyijy_{ij} , die vom jjj ten Geschwister in der iii ten Familie gemessen wurde . Darüber hinaus sind einige Verhaltensdaten xijxijx_{ij} in der gleichen Zeit von jedem Probanden erhoben wurden. Ich versuche die Situation mit dem folgenden linearen...
Nehmen wir an, wir haben ein Modell mod <- Y ~ X*Condition + (X*Condition|subject) # Y = logit variable # X = continuous variable # Condition = values A and B, dummy coded; the design is repeated # so all participants go through both Conditions # subject = random effects for different subjects...
Angenommen , ich habe Teilnehmer, von denen jeder eine Antwort gibt mal 20, 10 in einem Zustand und 10 in einem anderen. Ich passe ein lineares Mischeffektmodell an, das in jeder Bedingung vergleicht. Hier ist ein reproduzierbares Beispiel, das diese Situation anhand des Pakets in simuliert :Y...
Ich habe ein Experiment durchgeführt, bei dem ich verschiedene Familien aus zwei verschiedenen Bevölkerungsgruppen großgezogen habe. Jede Familie erhielt eine von zwei Behandlungen. Nach dem Experiment habe ich mehrere Merkmale an jedem Individuum gemessen. Um die Wirkung einer Behandlung oder...
Ich habe einen unsymmetrischen Datensatz mit wiederholten Messungen zur Analyse und ich habe gelesen, dass die Art und Weise, wie die meisten statistischen Pakete mit ANOVA umgehen (dh die Summe der Quadrate des Typs III), falsch ist. Daher würde ich gerne ein gemischtes Effektmodell verwenden, um...
Ich habe eine Frage zu meiner Verwendung eines gemischten Modells / Modells. Das Grundmodell ist folgendes: lmer(DV ~ group * condition + (1|pptid), data= df) Gruppe und Bedingung sind beide Faktoren: Die Gruppe hat zwei Ebenen (GruppeA, GruppeB) und die Bedingung hat drei Ebenen (Bedingung1,...
Douglas Bates gibt an, dass die folgenden Modelle äquivalent sind, "wenn die Varianz-Kovarianz-Matrix für die vektorwertigen Zufallseffekte eine spezielle Form hat, die als zusammengesetzte Symmetrie bezeichnet wird" ( Folie 91 in dieser Präsentation ): m1 <- lmer(y ~ factor + (0 +...
Ich führe ein logistisches Regressionsmodell in der folgenden Form aus: lmer(response~1+(1|site), family=binomial, REML = FALSE) Normalerweise würde ich den ICC aus dem Achsenabschnitt und den Restvarianzen berechnen, aber die Zusammenfassung des Modells enthält keine Restvarianzen. Wie berechne...
Ich habe über die Berechnung von R2R2R^2 -Werten in gemischten Modellen gelesen und nach dem Lesen der R-sig-FAQ, anderer Beiträge in diesem Forum (ich würde ein paar verlinken, aber ich habe nicht genug Ruf) und einiger anderer Referenzen, die ich unter Verwendung von verstehe 2R2R2R^2 Werte im...
Ich möchte 95% -Konfidenzintervalle für die Vorhersagen eines nichtlinearen gemischten nlmeModells erhalten. Da dies innerhalb von nichts Standardmäßigem vorgesehen ist, habe nlmeich mich gefragt, ob es richtig ist, die Methode der "Bevölkerungsvorhersageintervalle" zu verwenden, die in Ben Bolkers...
Meine Daten werden hier beschrieben. Was kann dazu führen, dass beim Anpassen einer ANOVA mit wiederholten Messungen ein "Error () - Modell ist ein singulärer Fehler" auftritt? Ich versuche, den Effekt einer Interaktion zu sehen, indem lmerich folgenden Basisfall benutze : my_null.model <-...
Ich habe folgendes Modell: > model1<-lmer(aph.remain~sMFS1+sAG1+sSHDI1+sbare+season+crop +(1|landscape),family=poisson) ... und das ist die Zusammenfassung. > summary(model1) Generalized linear mixed model fit by the Laplace approximation Formula: aph.remain ~ sMFS1 + sAG1 + sSHDI1 +...
TL; DR: Die lme4Optimierung scheint standardmäßig in Bezug auf die Anzahl der Modellparameter linear zu sein und ist viel langsamer als ein äquivalentes glmModell mit Dummy-Variablen für Gruppen. Kann ich irgendetwas tun, um es zu beschleunigen? Ich versuche, ein ziemlich großes hierarchisches...
Ich überarbeite ein Papier zur Bestäubung, in dem die Daten binomial verteilt sind (Obst reift oder nicht). Also habe ich glmermit einem zufälligen Effekt (einzelne Pflanze) und einem festen Effekt (Behandlung) gearbeitet. Ein Gutachter möchte wissen, ob die Pflanze einen Einfluss auf den...