Ich habe kürzlich 4 multiple Regressionsmodelle für dieselben Prädiktor- / Antwortdaten angepasst. Zwei der Modelle passen zur Poisson-Regression. model.pois <- glm(Response ~ P1 + P2 +...+ P5, family=poisson(), ...) model.pois.inter <- glm(Response ~ (P1 + P2 +...+ P5)^2, family=poisson(),...
Ich reproduziere ein Beispiel aus verallgemeinerten, linearen und gemischten Modellen . Mein MWE ist unten: Dilution <- c(1/128, 1/64, 1/32, 1/16, 1/8, 1/4, 1/2, 1, 2, 4) NoofPlates <- rep(x=5, times=10) NoPositive <- c(0, 0, 2, 2, 3, 4, 5, 5, 5, 5) Data <- data.frame(Dilution,...
Ich arbeite an einem Datensatz. Nachdem ich einige Modellidentifikationstechniken angewendet hatte, kam ich mit einem ARIMA (0,2,1) -Modell heraus. Ich habe die detectIOFunktion im Paket TSAin R verwendet, um bei der 48. Beobachtung meines ursprünglichen Datensatzes einen innovativen Ausreißer (IO)...
Ich verstehe, dass ich, wenn ich zwei Modelle A und B habe und A in B verschachtelt ist, anhand einiger Daten die Parameter von A und B mithilfe von MLE anpassen und den verallgemeinerten Log-Likelihood-Ratio-Test anwenden kann. Insbesondere sollte die Verteilung des Tests mit Freiheitsgraden sein,...
Nach dem, was ich gelesen habe, sind unter anderem auf der Website der UCLA-Statistikberatungsgruppe Likelihood-Ratio-Tests und Wald-Tests ziemlich ähnlich, wenn getestet wird, ob zwei glm-Modelle einen signifikanten Unterschied in der Passform für einen Datensatz aufweisen (entschuldigen Sie, wenn...
Ich habe einen Datensatz, in dem ich auf signifikante Unterschiede zwischen drei Populationen in Bezug auf etwa 50 verschiedene Variablen teste. Ich mache dies einerseits mit Kruskal-Wallis-Tests und andererseits mit Likelihood-Ratio-Tests für verschachtelte GLM-Modellanpassungen (mit und ohne...
Ich benutze die glmfitFunktion in MATLAB. Die Funktion gibt nur die Abweichung und nicht die Protokollwahrscheinlichkeit zurück. Ich verstehe, dass die Abweichung im Grunde doppelt so groß ist wie der Unterschied zwischen den Log-Wahrscheinlichkeiten der Modelle, aber was ich nicht bekomme, ist,...
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...
Ich arbeite derzeit daran, eine Methode zu implementieren, die in einem populären Artikel mit dem Titel "Ich habe gerade zwei Millionen Regressionen durchgeführt" verwendet wird. Die Grundidee dahinter ist, dass es bestimmte Fälle gibt, in denen nicht klar ist, welche Steuerelemente in das Modell...
Ich habe Daten mit einem Doppelpeak, die ich zu modellieren versuche, und es gibt genügend Überlappungen zwischen den Peaks, sodass ich sie nicht unabhängig behandeln kann. Ein Histogramm der Daten könnte ungefähr so aussehen: Ich habe dafür zwei Modelle erstellt: eines verwendet zwei...
Sei und zwei unabhängige Zufallsvariablen mit entsprechenden PDFs:XXXYYY f(x;θi)={1θie−x/θi0<x<∞,0<θi<∞0elsewheref(x;θi)={1θie−x/θi0<x<∞,0<θi<∞0elsewheref \left(x;\theta_i \right) =\begin{cases} \frac{1}{\theta_i} e^{-x/ {\theta_i}} \quad...
Ich analysiere einen bestimmten Datensatz und muss verstehen, wie ich das beste Modell für meine Daten auswähle. Ich benutze R. Ein Beispiel für Daten, die ich habe, ist das Folgende: corr <- c(0, 0, 10, 50, 70, 100, 100, 100, 90, 100, 100) Diese Zahlen entsprechen dem Prozentsatz der richtigen...
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 habe zwei unabhängige Poisson-Zufallsvariablen, und X 2 , mit X 1 ∼ Pois ( λ 1 ) und X 2 ∼ Pois ( λ 2 ) . Ich möchte H 0 testen :X1X1X_1X2X2X_2X1∼Pois(λ1)X1∼Pois(λ1)X_1 \sim \text{Pois}(\lambda_1)X2∼Pois(λ2)X2∼Pois(λ2)X_2 \sim \text{Pois}(\lambda_2) gegenüber der Alternative H 1...
Ich arbeite an einem Projekt, in dem wir das Verhalten einer Aufgabe beobachten (z. B. Reaktionszeit) und dieses Verhalten als Funktion mehrerer experimentell manipulierter Variablen sowie mehrerer beobachteter Variablen (Geschlecht der Teilnehmer, IQ der Teilnehmer, Antworten auf eine Folge)...
Die Software, die ich derzeit zum Erstellen eines Modells verwende, vergleicht ein "aktuelles Lauf" -Modell mit einem "Referenzmodell" und meldet (falls zutreffend) sowohl einen Chi-Quadrat-p-Wert basierend auf Likelihood-Ratio-Tests als auch AIC-Werte für jedes Modell. Ich weiß, dass ein Vorteil...