Statistiken und Big Data

Fragen und Antworten für Personen, die sich für Statistik, maschinelles Lernen, Datenanalyse, Data Mining und Datenvisualisierung interessieren


1
Logistische Regression mit Regressionssplines in R
Ich habe ein logistisches Regressionsmodell entwickelt, das auf retrospektiven Daten aus einer nationalen Traumadatenbank für Kopfverletzungen in Großbritannien basiert. Das Hauptergebnis ist die 30-Tage-Mortalität (als "Überlebensmaß" bezeichnet). Andere Maßnahmen mit veröffentlichten Hinweisen auf signifikante Auswirkungen auf das Ergebnis in früheren Studien umfassen: Year - Year of procedure = 1994-2013 Age …

3
Maximum-Likelihood-Schätzer der gemeinsamen Verteilung bei nur marginalen Zählungen
Sei px,ypx,yp_{x,y} eine gemeinsame Verteilung zweier kategorialer Variablen X,YX,YX,Y mit x,y∈{1,…,K}x,y∈{1,…,K}x,y\in\{1,\ldots,K\} . Angenommen, es wurden nnn Stichproben aus dieser Verteilung gezogen, aber wir erhalten nur die Grenzwerte, nämlich für j=1,…,Kj=1,…,Kj=1,\ldots,K : Sj=∑i=1nδ(Xi=l),Tj=∑i=1nδ(Yi=j),Sj=∑i=1nδ(Xi=l),Tj=∑i=1nδ(Yi=j), S_j = \sum_{i=1}^{n}{\delta(X_i=l)}, T_j = \sum_{i=1}^{n}{\delta(Y_i=j)}, Was ist der Maximum-Likelihood-Schätzer für px,ypx,yp_{x,y} bei Sj,TjSj,TjS_j,T_j ? Ist das bekannt? …

1
Extrahieren von Steigungen für Fälle aus einem Mixed-Effects-Modell (lme4)
Ich möchte die Steigungen für jedes Individuum in einem Mischeffektmodell extrahieren, wie im folgenden Absatz erläutert Modelle mit gemischten Effekten wurden verwendet, um einzelne Änderungspfade in den kognitiven Übersichtsmaßen zu charakterisieren, einschließlich Begriffen für Alter, Geschlecht und Bildungsjahr als feste Effekte (Laird und Ware, 1982; Wilson et al., 2000, 2002c). …
12 r  mixed-model 


1
So finden Sie ein GBM-Vorhersageintervall
Ich arbeite mit GBM-Modellen unter Verwendung des Caret-Pakets und suche nach einer Methode, um die Vorhersageintervalle für meine vorhergesagten Daten zu lösen. Ich habe intensiv gesucht, aber nur ein paar Ideen gefunden, um Vorhersageintervalle für Random Forest zu finden. Jeder Hilfe / R-Code wäre sehr dankbar!

3
Wie kann ich einen Plot des Trimmprozentsatzes im Vergleich zum getrimmten Mittelwert interpretieren?
Für einen Teil einer Hausaufgabenfrage wurde ich gebeten, den getrimmten Mittelwert für einen Datensatz durch Löschen der kleinsten und größten Beobachtung zu berechnen und das Ergebnis zu interpretieren. Der getrimmte Mittelwert war niedriger als der nicht getrimmte Mittelwert. Meine Interpretation war, dass dies daran lag, dass die zugrunde liegende Verteilung …

4
Welches Verhältnis der unabhängigen Verteilungen ergibt eine Normalverteilung?
Das Verhältnis zweier unabhängiger Normalverteilungen ergibt eine Cauchy-Verteilung. Die t-Verteilung ist eine Normalverteilung geteilt durch eine unabhängige Chi-Quadrat-Verteilung. Das Verhältnis von zwei unabhängigen Chi-Quadrat-Verteilungen ergibt eine F-Verteilung. Ich suche ein Verhältnis unabhängiger stetiger Verteilungen, das eine normalverteilte Zufallsvariable mit mittlerem μμ\mu und Varianz ergibt ?σ2σ2\sigma^2 Es gibt wahrscheinlich unendlich viele …

3
Grundlegendes zu Gaußschen Basisfunktionsparametern für die lineare Regression
Ich möchte die Gaußsche Basisfunktion auf eine lineare Regressionsimplementierung anwenden. Leider fällt es mir schwer, ein paar Parameter in der Basisfunktion zu verstehen. Insbesondere und .μμ\muσσ\sigma Mein Datensatz ist eine 10.000 x 31-Matrix. 10.000 Samples und 31 Features. Ich habe gelesen, dass "Jede Basisfunktion den Eingabevektor x in einen Skalarwert …

1
Determinante der Fisher-Information
(Ich habe eine ähnliche Frage auf math.se gestellt .) In der Informationsgeometrie ist die Determinante der Fisher-Informationsmatrix eine natürliche Volumenform auf einer statistischen Mannigfaltigkeit, daher hat sie eine schöne geometrische Interpretation. Die Tatsache, dass es beispielsweise in der Definition eines Jeffreys vorkommt, hängt mit seiner Invarianz unter Reparametrisierungen zusammen, die …

1
t.test gibt einen Fehler zurück "Daten sind im Wesentlichen konstant"
R version 3.1.1 (2014-07-10) -- "Sock it to Me" > bl <- c(140, 138, 150, 148, 135) > fu <- c(138, 136, 148, 146, 133) > t.test(fu, bl, alternative = "two.sided", paired = TRUE) Error in t.test.default(fu, bl, alternative = "two.sided", paired = TRUE) : data are essentially constant Dann …
12 r  t-test 

4
Gute Form, um Ausreißer zu entfernen?
Ich arbeite an Statistiken für Software-Builds. Ich habe Daten für jeden Build auf Pass / Fail und abgelaufene Zeit und wir generieren ~ 200 davon / Woche. Die Erfolgsquote lässt sich leicht zusammenfassen. Ich kann sagen, dass 45% einer Woche vergangen sind. Aber ich möchte auch die verstrichene Zeit zusammenfassen …

2
Warum gibt mir SAS PROC GLIMMIX SEHR unterschiedliche Zufallsneigungen als glmer (lme4) für ein Binomial-Glmm
Ich bin ein mit R vertrauterer Benutzer und habe versucht, zufällige Steigungen (Auswahlkoeffizienten) für ungefähr 35 Personen über 5 Jahre für vier Lebensraumvariablen zu schätzen. Die Antwortvariable gibt an, ob ein Standort als Lebensraum "verwendet" (1) oder "verfügbar" (0) war ("Verwendung" unten). Ich verwende einen Windows 64-Bit-Computer. In R-Version 3.1.0 …

3
Testen von logistischen Regressionskoeffizienten unter Verwendung von Freiheitsgraden
Zusammenfassung: Gibt es eine statistische Theorie, die die Verwendung der Verteilung (mit Freiheitsgraden basierend auf der Restabweichung) für Tests von logistischen Regressionskoeffizienten anstelle der Standardnormalverteilung unterstützt?ttt Vor einiger Zeit habe ich festgestellt, dass beim Anpassen eines logistischen Regressionsmodells in SAS PROC GLIMMIX unter den Standardeinstellungen die logistischen Regressionskoeffizienten unter Verwendung …

1
Erfassen der Ergebnisse der Mediationsanalyse in R
Ich versuche, mich mit dem Mediationspaket in R vertraut zu machen, indem ich die Vignette für das Paket verwende. Ich habe Mühe, die Ausgabe der mediate()Funktion zu verstehen . require("mediation") require("sandwich") data("framing") med.fit <- lm(emo ~ treat + age + educ + gender + income, data = framing) out.fit <- …
12 r  mediation 

Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.