Eine diskrete, univariate Verteilung, die die Anzahl von modelliert B e r n o u l l i (p) Versuchserfolge, bis eine bestimmte Anzahl von Fehlern auftritt.
Ich möchte also ein Negativ-Binomial-Modell mit zufälligen Effekten anpassen. Für ein solches Modell kann STATA potenzierte Koeffizienten erzeugen. Der Hilfedatei zufolge können solche Koeffizienten als Inzidenzratenverhältnisse interpretiert werden. Leider bin ich kein englischer Muttersprachler und verstehe nicht wirklich, welche Inzidenzratenverhältnisse vorliegen und wie ich sie übersetzen könnte. Meine Frage ist …
Ich habe drei Vorschläge zur Behandlung der Überdispersion in einer Poisson-Antwortvariablen und einem Startmodell mit allen festen Effekten gefunden: Verwenden Sie ein Quasi-Modell. Verwenden Sie negatives Binomial-GLM. Verwenden Sie ein gemischtes Modell mit einem zufälligen Effekt auf Subjektebene. Aber was soll man eigentlich wählen und warum? Gibt es ein tatsächliches …
Ich bin auf eine Reihe praktischer Fragen gestoßen, wenn ich Zähldaten aus experimenteller Forschung mithilfe eines subjektinternen Experiments modelliere. Ich beschreibe kurz das Experiment, die Daten und das, was ich bisher gemacht habe, gefolgt von meinen Fragen. Einer Stichprobe von Befragten wurden nacheinander vier verschiedene Filme gezeigt. Nach jedem Film …
Ich habe einen sehr kleinen Datensatz zur Bienenfülle, den ich nur schwer analysieren kann. Es handelt sich um Zählungsdaten, und fast alle Zählungen erfolgen in einer Behandlung, wobei die meisten Nullen in der anderen Behandlung vorkommen. Es gibt auch einige sehr hohe Werte (jeweils einer an zwei der sechs Stellen), …
Ich habe einen Zeitreihendatensatz, an den ich ein Hidden Markov Model (HMM) anpasse, um die Anzahl der latenten Zustände in den Daten abzuschätzen. Mein Pseudocode dafür ist der folgende: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ... optimal_number_of_states = "model with smallest BIC" ... …
Die Frage ist folgende: Eine Zufallsstichprobe von n Werten wird aus einer negativen Binomialverteilung mit dem Parameter k = 3 gesammelt. Finden Sie den Maximum-Likelihood-Schätzer des Parameters π. Finden Sie eine asymptotische Formel für den Standardfehler dieses Schätzers. Erklären Sie, warum die negative Binomialverteilung ungefähr normal ist, wenn der Parameter …
Ich versuche, die mittleren Intensitäten von Parasiten, die einen Wirt in R beeinflussen, mithilfe eines negativen Binomialmodells zu modellieren. Ich bekomme immer wieder 50 oder mehr Warnungen, die besagen: In dpois(y, mu, log = TRUE) : non-integer x = 251.529000 Wie kann ich damit umgehen? Mein Code sieht folgendermaßen aus: …
In diesem Artikel (frei verfügbar über PubMed Central) verwenden die Autoren eine negative binomiale Regression, um die Punktzahl auf einem 10-Punkte-Screening-Instrument mit einer Bewertung von 0 bis 40 zu modellieren. Bei diesem Verfahren werden Zähldaten vorausgesetzt, was hier eindeutig nicht der Fall ist. Ich würde gerne Ihre Meinung dazu erfahren, …
Ich habe ein GLMM der Form: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Wenn ich benutze drop1(model, test="Chi"), erhalte ich andere Ergebnisse als wenn ich Anova(model, type="III")aus dem Autopaket oder benutze summary(model). Diese beiden letzteren geben die gleichen Antworten. Unter Verwendung einer Reihe …
Ich versuche, die Arbeit eines Kollegen zu replizieren, und verschiebe die Analyse von Stata nach R. Die von ihr verwendeten Modelle rufen die Option "Cluster" innerhalb der Funktion nbreg auf, um die Standardfehler zu gruppieren. Unter http://repec.org/usug2007/crse.pdf finden Sie eine ziemlich vollständige Beschreibung des Was und Warum dieser Option Meine …
Ich habe negative Binomial- und Quasi-Poisson-Modelle ausgeführt, die auf einem Ansatz zum Testen von Hypothesen basieren. Meine endgültigen Modelle, die beide Methoden verwenden, haben unterschiedliche Kovariaten und Wechselwirkungen. Es scheint, dass es in beiden Fällen keine Muster gibt, wenn ich meine Residuen zeichne. Daher habe ich mich gefragt, mit welchem …
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 verwende US-Volkszählungsdaten aus der Longitudinal …
Ist es möglich, ein GAMM (Generalized Additive Mixed Model) für Daten ohne Inflation in R anzupassen? Wenn nicht, ist es möglich, ein GAM (Generalized Additive Model) für Daten ohne Inflation mit einer negativen Binomial- oder Quasi-Poisson-Verteilung in R anzupassen? (Ich fand COZIGAM :: zigam- und mgcv: ziP- Funktionen für die …
Ich hatte eine Hausaufgabe, um die negative Binomialverteilung als exponentielle Verteilungsfamilie auszudrücken, da der Dispersionsparameter eine bekannte Konstante war. Das war ziemlich einfach, aber ich fragte mich, warum sie erfordern würden, dass wir diesen Parameter festhalten. Ich stellte fest, dass ich keinen Weg finden konnte, es in die richtige Form …
Ich habe diese Daten: set.seed(1) predictor <- rnorm(20) set.seed(1) counts <- c(sample(1:1000, 20)) df <- data.frame(counts, predictor) Ich habe eine Poisson-Regression durchgeführt poisson_counts <- glm(counts ~ predictor, data = df, family = "poisson") Und eine negative binomiale Regression: require(MASS) nb_counts <- glm.nb(counts ~ predictor, data = df) Dann berechnete ich …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.