Wenn bekannte Dichten sind, aus denen ich simulieren kann, dh für die ein Algorithmus verfügbar ist. und wenn das Produkt integrierbar ist, gibt es einen generischen Ansatz, um aus dieser Produktdichte unter Verwendung von zu simulieren Simulatoren aus dem 's?k ∏ i = 1 f i ( x ) α...
Ich verstehe, dass der Unterschied zwischen ihnen damit zusammenhängt, ob die Gruppierungsvariable im Modell als fester oder zufälliger Effekt geschätzt wird, aber mir ist nicht klar, warum sie nicht gleich sind (wenn sie nicht gleich sind). Ich bin speziell daran interessiert, wie dies bei der...
Wir haben viele gute Diskussionen über die perfekte Trennung in der logistischen Regression. So führte die logistische Regression in R zu einer perfekten Trennung (Hauck-Donner-Phänomen). Was jetzt? Das logistische Regressionsmodell konvergiert nicht . Ich persönlich glaube immer noch, dass es...
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...
Ist die Neugestaltung eines multivariaten linearen Regressionsmodells als multiple lineare Regression völlig gleichwertig? Ich beziehe mich nicht einfach laufen separate Regressionen.ttt Ich habe an einigen Stellen gelesen (Bayesian Data Analysis - Gelman et al. Und Multivariate Old School -...
Ich habe es mit einem Bayesian Hierarchical Linear Model zu tun , hier das Netzwerk, das es beschreibt. YYY für den täglichen Verkauf eines Produkts in einem Supermarkt (beobachtet). XXX ist eine bekannte Matrix von Regressoren, einschließlich Preisen, Werbeaktionen, Wochentagen, Wetter,...
Wenn ich die Werte der beobachteten Knoten eines MRF festlege, wird es dann zu einem
Ich benötige einige Ratschläge in Bezug auf zwei Hauptprobleme in meiner Forschung, die eine Fallstudie von drei großen Pharmazeutika und Innovationen ist. Anzahl der Patente pro Jahr ist die abhängige Variable. Meine Fragen sind Was sind die wichtigsten Kriterien für ein gutes Modell? Was ist mehr...
Das scheint so elementar zu sein, aber ich bleibe immer an diesem Punkt stecken ... Die meisten Daten, mit denen ich zu tun habe, sind nicht normal, und die meisten Analysen basieren auf einer GLM-Struktur. Für meine aktuelle Analyse habe ich eine Antwortvariable, die "Gehgeschwindigkeit" (Meter /...
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...
Angenommen, ich habe einen "Nierenkatheter" -Datensatz. Ich versuche, eine Überlebenskurve mit einem Cox-Modell zu modellieren. Wenn ich ein Cox-Modell betrachte: brauche ich die Schätzung der . Mit der eingebauten Paket-R-Funktion kann ich das ganz einfach so
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...
Ich habe eine Tabelle mit einigen Mannschaftswerten. Erste Mannschaft mit 10 Punkten gewinnt. Jedes Team besteht aus 2 Spielern. Die Spieler spielen die ganze Zeit mit verschiedenen Teammitgliedern, obwohl sie nicht vollkommen zufällig ausgewählt werden. Es werden keine Einzelbewertungen...
Ich interessiere mich für das Bewerten und Einordnen von Personen in einer Gruppe, die nur paarweise interagieren / gegeneinander antreten (dh Systeme wie das ELO- Bewertungssystem für Schach). Gibt es Go-to-Methoden oder genauere und fortgeschrittenere Methoden? Gibt es R-Pakete, die die...
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,...
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...
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...
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...
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 =...
Ich habe eine Darstellung der Restwerte eines linearen Modells in Abhängigkeit von den angepassten Werten, wobei die Heteroskedastizität sehr klar ist. Ich bin mir jedoch nicht sicher, wie ich jetzt vorgehen soll, da diese Heteroskedastizität meines Wissens mein lineares Modell ungültig macht. (Ist...