Statistiken und Big Data r

2

Wie berechnet man das Konfidenzintervall des x-Achsenabschnitts in einer linearen Regression?

Da der Standardfehler einer linearen Regression normalerweise für die Antwortvariable angegeben wird, frage ich mich, wie Konfidenzintervalle in die andere Richtung erhalten werden können - z. B. für einen x-Achsenabschnitt. Ich kann mir vorstellen, was es sein könnte, aber ich bin sicher, dass es einen einfachen Weg geben muss, dies …

9 r regression confidence-interval bootstrap

1

Definition der gewichteten kleinsten Quadratgewichte: R lm Funktion vs.

Kann mir jemand sagen, warum ich durch Rgewichtete kleinste Quadrate und manuelle Lösung durch Matrixoperation unterschiedliche Ergebnisse erhalte ? Insbesondere versuche ich, manuell zu lösen , wobei die Diagonalmatrix für Gewichte ist, die Datenmatrix ist, die Antwort ist Vektor. WAx=WbWAx=Wb\mathbf W \mathbf A\mathbf x=\mathbf W \mathbf bWW\mathbf WAA\mathbf Abb\mathbf b …

9 r regression least-squares weighted-regression weighted-data

2

So verbessern Sie die Laufzeit für die R MICE-Datenimputation

Meine Frage kurz: Gibt es Methoden zur Verbesserung der Laufzeit von R MICE (Datenimputation)? Ich arbeite mit einem Datensatz (30 Variablen, 1,3 Millionen Zeilen), der (ziemlich zufällig) fehlende Daten enthält. Etwa 8% der Beobachtungen in etwa 15 von 30 Variablen enthalten NAs. Um die fehlenden Daten zu unterstellen, führe ich …

9 r multiple-imputation mice

2

Neuere Beobachtungen in der Regression stärker gewichten

Wie ordne ich neueren Beobachtungen in R mehr Gewicht zu? Ich nehme dies als häufig gestellte Frage oder Wunsch an, aber es fällt mir schwer, genau herauszufinden, wie ich dies umsetzen soll. Ich habe versucht, viel danach zu suchen, aber ich kann kein gutes praktisches Beispiel finden. In meinem Beispiel …

9 r time-series regression random-forest

2

Was ist der Unterschied zwischen diesen beiden Breusch-Pagan-Tests?

Wenn ich R für einige Daten verwende und versuche festzustellen, ob meine Daten heteroskedastisch sind oder nicht, habe ich zwei Implementierungen des Breusch-Pagan-Tests gefunden, bptest (Paket lmtest) und ncvTest (Paketauto). Diese führen jedoch zu unterschiedlichen Ergebnissen. Was ist der Unterschied zwischen den beiden? Wann sollten Sie sich für das eine …

9 r regression heteroscedasticity assumptions breusch-pagan

1

Anpassen eines binomischen GLMM (glmer) an eine Antwortvariable, die ein Anteil oder ein Bruchteil ist

Ich hoffe, jemand kann bei einer meiner Meinung nach relativ einfachen Frage helfen, und ich glaube, ich kenne die Antwort, aber ohne Bestätigung ist sie zu etwas geworden, dessen ich mir einfach nicht sicher sein kann. Ich habe einige Zähldaten als Antwortvariable und möchte messen, wie sich diese Variable mit …

9 r binomial proportion glmm lme4-nlme

2

Funktionsauswahl für Clustering-Probleme

Ich versuche, verschiedene Datensätze mit unbeaufsichtigten Algorithmen (Clustering) zu gruppieren. Das Problem ist, dass ich viele Funktionen (~ 500) und eine kleine Anzahl von Fällen (200-300) habe. Bisher habe ich nur Klassifizierungsprobleme gemacht, für die ich Daten immer als Trainingssätze gekennzeichnet hatte. Dort habe ich ein Kriterium (dh random.forest.importance oder …

9 r clustering feature-selection unsupervised-learning

1

Warum ersetzt die Intercept-Spalte in model.matrix den ersten Faktor?

Ich versuche, meine Faktorspalte in Dummy-Variablen umzuwandeln: str(cards$pointsBin) # Factor w/ 5 levels ".lte100",".lte150",..: 3 2 3 1 4 4 2 2 4 4 ... labels <- model.matrix(~ pointsBin, data=cards) head(labels) # (Intercept) pointsBin.lte150 pointsBin.lte200 pointsBin.lte250 pointsBin.lte300 # 741 1 0 0 0 0 # 407 1 1 0 0 …

9 r categorical-data categorical-encoding

1

Wie viele Seiten hat ein Würfel? Bayesianische Folgerung in JAGS

Problem Ich möchte auf ein System schließen, das analog dazu ist, mit einer unbekannten Anzahl von Seiten zu sterben. Der Würfel wird mehrmals gewürfelt, wonach ich eine Wahrscheinlichkeitsverteilung über einen Parameter ableiten möchte, der der Anzahl der Seiten des Würfels entspricht, θ. Intuition Wenn Sie nach 40 Rollen 10 Rot-, …

9 r probability bayesian jags

2

Interpretation und Durchführung von Prognosen mit dem tsoutliers-Paket und auto.arima

Ich habe monatliche Daten von 1993 bis 2015 und möchte diese Daten prognostizieren. Ich habe das tsoutliers-Paket verwendet, um die Ausreißer zu erkennen, aber ich weiß nicht, wie ich mit meinen Daten weiter prognostizieren soll. Das ist mein Code: product.outlier<-tso(product,types=c("AO","LS","TC")) plot(product.outlier) Dies ist meine Ausgabe vom tsoutliers-Paket ARIMA(0,1,0)(0,0,1)[12] Coefficients: sma1 …

9 r time-series forecasting arima outliers

1

So stellen Sie benutzerdefinierte Kontraste mit lmer in R ein

Ich verwende lmer in R, um die Auswirkung von condition ( cond) auf ein Ergebnis zu überprüfen . Hier sind einige aus Daten, wobei s der Objektidentifizierer und a, bund csind Bedingungen. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30), a = rnorm(30, -2, 1), b = rnorm(30, -3, …

9 r lme4-nlme contrasts

2

Was ist VectorSource und VCorpus im Paket 'tm' (Text Mining) in R.

Ich bin mir nicht ganz sicher, was genau VectorSource und VCorpus im 'tm'-Paket enthalten. Die Dokumentation dazu ist unklar. Kann mich jemand in einfachen Worten verstehen lassen?

9 r text-mining

1

Kreuzvalidierende Lasso-Regression in R.

Die R-Funktion cv.glm (Bibliothek: Boot) berechnet den geschätzten K-fachen Kreuzvalidierungs-Vorhersagefehler für verallgemeinerte lineare Modelle und gibt Delta zurück. Ist es sinnvoll, diese Funktion für eine Lasso-Regression (Bibliothek: glmnet) zu verwenden, und wenn ja, wie kann sie ausgeführt werden? Die glmnet-Bibliothek verwendet eine Kreuzvalidierung, um den besten Drehparameter zu erhalten, aber …

9 r regression cross-validation lasso glmnet

1

Schätzen mehrstufiger logistischer Regressionsmodelle

Das folgende mehrstufige Logistikmodell mit einer erklärenden Variablen auf Ebene 1 (Einzelebene) und einer erklärenden Variablen auf Ebene 2 (Gruppenebene): π 0 j = γ 00 + γ 01 z j + u 0 j … ( 2 ) π 1 j = γ 10 + γ 11 z j …

9 r logistic generalized-linear-model simulation multilevel-analysis

1

Wie verwende ich Anova für den Vergleich zweier Modelle?

Wie soll ich das anovaErgebnis beim Vergleich zweier Modelle verstehen ? Beispiel: Res.Df RSS Df Sum of Sq F Pr(>F) 1 9 54.032 2 7 4.632 2 49.4 37.329 0.0001844 *** In der Manpage heißt es: "Berechnen Sie die Analyse von Varianz- (oder Abweichungs-) Tabellen für ein oder mehrere angepasste …

9 r regression anova

Als «r» getaggte Fragen