Tag: glm

Für Fragen im Zusammenhang mit verallgemeinerten linearen Modellen. Für die GLM-math-Bibliothek, siehe [tag:glm-math].

Statsmodels Poisson-glm anders als R

Anzahl der Antworten 1 Antworten
Ich versuche, passen einige Modelle (Spatial interaction models) laut einigen code, die in R. ich habe in der Lage zu bekommen, einige der code, um die Arbeit mit statsmodels in einem python-framework, sondern einige von Ihnen überhaupt

Logistische regression - cbind-Befehl in glm

Anzahl der Antworten 1 Antworten
Ich mache Logistische regression in R. Kann jemand klären, was die Unterschiede der Ausführung dieser beiden Zeilen? 1. glm(Response ~ Temperature, data=temp, family = binomial(link="logit")) 2. glm(cbind(Response, n - Response) ~ Temperature, data=temp, family =binomial, Ntrials=n) Den

Berechnen Sie die cross-Validierung für das Verallgemeinerte Lineare Modell in Matlab

Anzahl der Antworten 1 Antworten
Mache ich eine regression mit Verallgemeinerten Linearen Modell.Ich bin gefangen offguard mit der crossVal Funktion. Meine Implementierung so weit; x = 'Some dataset, containing the input and the output' X = x(:,1:7); Y = x(:,8); cvpart =

Polynom-Daten und R glm()

Anzahl der Antworten 2 Antworten
Wie kann man R mit der glm() zu entsprechen Polynom Daten? Ich habe versucht, mehrere Iterationen von " Familie=AAA(link="BBB")' aber ich kann nicht scheinen, um triviale Vorhersagen übereinstimmen. Zum Beispiel, bitte helfen Sie mit R glm passend

Plotten glm Interaktionen: "newdata=" Struktur in der predict () - Funktion

Anzahl der Antworten 1 Antworten
Mein problem ist mit der predict() Funktion, seiner Struktur und der Darstellung der Vorhersagen. Mit den Vorhersagen aus meinem Modell, würde ich mag, um zu visualisieren, wie mein signifikanten Faktoren (und deren Interaktion) Einfluss auf die Wahrscheinlichkeit,

Mit ROCR-Paket Schwierigkeiten

Anzahl der Antworten 2 Antworten
Ich bin nach einer Analyse der deutschen Kredit-Daten und ich bekam eine Fehlermeldung, die ich versage zu korrigieren, weil ich bereits installiert ROCR Paket. Unten ist der code, der verwendet die ROCR: #load library library(ROCR) #score test

predict() mit beliebigen Koeffizienten in r

Anzahl der Antworten 3 Antworten
Habe ich einige Koeffizienten eines logit-Modells, die von einer nicht-r-Benutzer. Ich möchte diese importieren Koeffizienten in r und generieren einige goodness-of-fit Schätzungen, die auf dem gleichen Datensatz (ROC und confusion matrix) vs mein eigenes Modell. Mein Erster

Ändern der Referenz-Gruppe mit binomial-glm mit Familie

Anzahl der Antworten 2 Antworten
Wenn ich eine binomial-regression in R mit einem Faktor unabhängig variable, bestehend aus drei Ebenen der "Höheren" , "Mittleren" und "Unteren", den ich ändern möchten die referenzkategorie mit relevel bekomme ich diesen Fehler: “Error in relevel.ordered(cbsnivcat3, "Lower")

Fehlermeldungen beim ausführen von glmer in R

Anzahl der Antworten 1 Antworten
Ich bin versucht, um zwei ähnliche, verallgemeinerte lineare gemischte Modelle in R. Beide Modelle haben die gleichen input-Variablen für Einflussfaktoren, kovariablen und zufällige Faktoren, jedoch, response-Variablen unterscheiden. Modelle erfordern die lme4 Paket. Das Problem hatte ich mit

Bekommen 95% - Konfidenzintervall mit glm(..) in R

Anzahl der Antworten 2 Antworten
Hier sind einige Daten dat = data.frame(y = c(9,7,7,7,5,6,4,6,3,5,1,5), x = c(1,1,2,2,3,3,4,4,5,5,6,6), color = rep(c('a','b'),6)) und der plot dieser Daten, wenn Sie es wünschen require(ggplot) ggplot(dat, aes(x=x,y=y, color=color)) + geom_point() + geom_smooth(method='lm') Wenn ein Modell mit der

R Fehler, der sagt "Modelle waren nicht alle montiert, um die gleiche Größe des datasets"

Anzahl der Antworten 5 Antworten
Erstellt habe ich zwei generalisierte lineare Modelle wie folgt: glm1 <-glm(Y ~ X1 + X2 + X3, family=binomial(link=logit)) glm2 <-glm(Y ~ X1 + X2, family=binomial(link=logit)) Dann benutze ich die anova Funktion: anova(glm2,glm1) bekomme aber eine Fehlermeldung: "Fehler

wie kann ich ausschließen bestimmter Variablen aus einer glm in R?

Anzahl der Antworten 1 Antworten
Habe ich 50 Variablen. Dies ist, wie ich Sie alle in meine glm. var = glm(Stuff ~ ., data=mydata, family=binomial) Aber ausschließen will ich 2 davon. Also wie kann ich ausschließen 2 in bestimmten? Ich hatte gehofft,

Wie melde Likelihood-Ratio-Test-Ergebnisse

Anzahl der Antworten 1 Antworten
Ich bin mit einem Likelihood-Ratio-Test (in R) zu sehen, die für die wichtigsten Effekte in meinem Modell mit drei festen Faktoren (Standort, Jahr, habitat) wie folgt: model1<-glm(tot.mass~hab, data=biom, family = Gamma(link = "log")) anova(model1, test="Chisq") model2<-glm(tot.mass~year, data=biom,

Mit der predict () - Funktion mit glm

Anzahl der Antworten 2 Antworten
Sagen wir mal so, ich habe folgende Daten eingestellt und führe eine regression mit glm in R. habe ich die Koeffizienten, aber ich möchte, um vorherzusagen, "nächsten Monaten" Wert (Besuche). Wie würde ich mich über, dass in

glm Startwerte nicht akzeptiert log-link

Anzahl der Antworten 2 Antworten
Möchte ich laufen eine Gauß-GLM mit log-link und einen offset. Die folgenden Probleme auftreten: y <- c(1,1,0,0) t <- c(5,3,2,4) Kein problem: exp(coef(glm(y~1 + offset(log(t)), family=poisson))) mit family=gaussian ab Werte müssen angegeben werden, es funktioniert hier: exp(coef(glm(y~1,

Die Entscheidung, die Schwelle für glm logistischen Regressionsmodell in R

Anzahl der Antworten 4 Antworten
Habe ich einige Daten mit Prädiktoren und eine binäre Ziel. Eg: df <- data.frame(a=sort(sample(1:100,30)), b= sort(sample(1:100,30)), target=c(rep(0,11),rep(1,4),rep(0,4),rep(1,11))) Trainierte ich ein Logistik-regresion-Modell mit glm() model1 <- glm(formula= target ~ a + b, data=df, family=binomial) Nun bin ich versucht,

Unterschied zwischen glmnet() und cv.glmnet() in R?

Anzahl der Antworten 2 Antworten
Arbeite ich an einem Projekt, das zeigen würde, den potenziellen Einfluss einer Gruppe von Ereignissen auf ein Ergebnis. Ich bin mit dem glmnet () - Paket, speziell mit Hilfe der Poisson-Funktion. Hier ist mein code: # de

Wie zum extrahieren von einem p-Wert bei der Durchführung der anova() zwischen zwei glm-Modelle in R

Anzahl der Antworten 2 Antworten
So, ich bin versucht, zwei Modelle vergleichen, fit1 und fit2. Erst war ich nur tun anova(fit1,fit2), und dort ergab sich eine Ausgabe, die ich verstanden habe (darunter ein p-Wert). Jedoch, wenn ich wechselte meine Modelle von lm()-based

Vorhersagen.glm nicht Vorhersage der fehlenden Werte in der Antwort

Anzahl der Antworten 3 Antworten
Aus irgendeinem Grund, wenn ich angeben glms (und lm ' s zu, es stellt sich heraus), R ist nicht die Vorhersage der fehlenden Werte der Daten. Hier ist ein Beispiel: y = round(runif(50)) y = c(y,rep(NA,50)) x

glmer - Vorhersagen mit binomiale Daten (cbind count data)

Anzahl der Antworten 1 Antworten
Ich versuche voraussagen, - Werte über die Zeit (Tage x-Achse) für eine glmer-Modell, das war auf meinem binomialen Daten. Total Lebendig und Total Tot sind count data. Dies ist mein Modell, und die entsprechenden Schritte unten. full.model.dredge<-glmer(cbind(Total.Alive,Total.Dead)~(CO2.Treatment+Lime.Treatment+Day)^3+(Day|Container)+(1|index),

Zum Lesen der Werte aus einer glm::mat4

Anzahl der Antworten 1 Antworten
Ich habe ein glm::mat4 matrix und ich brauche, um die Werte in ein double[16] - array. Irgendwelche Ideen auf, wie dieses problem zu lösen?? InformationsquelleAutor Edvin | 2013-09-19

Wie extrahiere ich lmer festen Effekte durch Beobachtung?

Anzahl der Antworten 2 Antworten
Habe ich ein lme-Objekt, gebaut aus einigen Maßnahmen wiederholt Nährstoff-Aufnahme der Daten (zwei 24-Stunden-Aufnahme-Perioden pro RespondentID): Male.lme2 <- lmer(BoxCoxXY ~ -1 + AgeFactor + IntakeDay + (1|RespondentID), data = Male.Data, weights = SampleWeight) und ich kann erfolgreich

Die Auswahl der statistisch signifikanten Variablen in einer R glm-Modell

Anzahl der Antworten 4 Antworten
Habe ich eine outcome-variable, sagen wir Y und eine Liste der 100 Dimensionen, die Einfluss auf Y (X1 sagen...X100). Nach dem laufen meine glm und zum anzeigen einer Zusammenfassung von meinem Modell, ich sehe die Variablen, die

Cross-Validierung für das glm() Modelle

Anzahl der Antworten 2 Antworten
Ich versuche zu tun, eine 10-fold-cross-validation für einige glm-Modelle, die ich gebaut habe früher in R. ich bin ein wenig verwirrt über die cv.glm() Funktion in der boot Paket, obwohl ich gelesen habe, viel Hilfe-Dateien. Wenn ich

Auszug Standardfehler von glm

Anzahl der Antworten 3 Antworten
Habe ich ein glm und ich will einfach nur, um zu extrahieren die Standardfehler der einzelnen Koeffizienten. Ich sah im internet die Funktion se.coef() aber es funktioniert nicht, es gibt "Error: could not find function "se.coef"". Vielleicht

Die Angabe der Formel in R mit glm ohne explizite Deklaration der einzelnen kovariable

Anzahl der Antworten 2 Antworten
Ich möchte zu zwingen, bestimmte Variablen in glm-Regressionen ohne Angabe jedes. Meine real-Daten festgelegt hat ~200 Variablen. Ich habe nicht in der Lage, um Proben von dieser in meinem online-Suche so weit. Beispiel (mit nur 3 Variablen):

Extrahiere pvalue aus Glm

Anzahl der Antworten 4 Antworten
Ich bin mit vielen Rückführungen und bin nur daran interessiert, die Auswirkungen auf die Koeffizienten und p-Wert einer bestimmten Variablen. Also in meinem script, ich möchte in der Lage sein, um Sie einfach extrahieren Sie den p-Wert

Wie schreibe ich eine Formel mit vielen Variablen aus einem Datenrahmen?

Anzahl der Antworten 6 Antworten
Angenommen, ich habe eine Antwort-Variablen und-Daten mit drei kovariablen (wie ein Spielzeug Beispiel): y = c(1,4,6) d = data.frame(x1 = c(4,-1,3), x2 = c(3,9,8), x3 = c(4,-4,-2)) Will ich fit einer linearen regression, zu den Daten: fit

caret train () sagt das sehr unterschiedlich vorher. vorhergesagt.glm ()

Anzahl der Antworten 1 Antworten
Ich versuche zur Schätzung einer logistischen regression unter Verwendung der 10-fold-cross-validation. #import libraries library(car); library(caret); library(e1071); library(verification) #data import and preparation data(Chile) chile <- na.omit(Chile) #remove "na's" chile <- chile[chile$vote == "Y" | chile$vote == "N" ,

Warnung: nicht ganzzahlige #successes in einem Binomial-Glm! (Umfragepakete)

Anzahl der Antworten 3 Antworten
Ich bin mit der twang - Paket zu erstellen, propensity scores, die als GEWICHTE in einem binomial-glm mit survey::svyglm. Der code sieht ungefähr so aus: pscore <- ps(ppci ~ var1+var2+.........., data=dt....) dt$w <- get.weights(pscore, stop.method="es.mean") design.ps <-

Konfidenzintervalle für Vorhersagen aus der logistischen Regression

Anzahl der Antworten 1 Antworten
In R Vorhersagen.lm berechnet Vorhersagen auf der Grundlage der Ergebnisse aus der linearen regression und bietet auch zur Berechnung von Konfidenzintervallen für diese Vorhersagen. Laut Handbuch diese Intervalle werden auf der Grundlage der Fehler-Varianz-fitting, aber nicht auf