Ich habe das Mäusepaket untersucht und noch keine Möglichkeit gefunden, die mehreren Imputationen zu verwenden, um ein Cox-Modell zu erstellen und dieses Modell dann mit der Funktion des Effektivpakets zu validieren validate(). Hier ist ein Beispielcode von dem, was ich bisher habe, unter...
Als Folgemaßnahme zum Modellvorschlag für eine Cox-Regression mit zeitabhängigen Kovariaten ist hier das Kaplan-Meier-Diagramm, das die zeitabhängige Natur von Schwangerschaften berücksichtigt. Mit anderen Worten, der Datensatz wird jetzt in einen langen Datensatz mit mehreren Zeilen entsprechend...
Hat jemand ein kurzes Beispiel, wie man neuronale Netze ( nnetzum Beispiel in R) zum Zwecke der Vorhersage verwendet? Hier ist ein Beispiel in R einer Zeitreihe T <- seq(0,20,length=200) Y <- 1 + 3*cos(4*T+2) +.2*T^2 + rnorm(200) plot(T,Y,type="l") Dies ist nur ein Beispiel, aber ich habe...
Ich verwende derzeit R, um eine Zeitreihe mit den folgenden Anweisungen vorherzusagen: X <- ts(datas, frequency=24) X.arima <- Arima(X, order=c(2,1,0), seasonal=c(1,1,1)) pred <- predict(X.arima, n.ahead=24) plot.ts(pred$pred) Wie Sie sehen, habe ich stündlich Daten und habe den...
Ich versuche, eine latente Klassenwachstumsanalyse (LCGA) und / oder Wachstumsmischungsmodelle (GMMs) in R durchzuführen. Die Daten, die ich verwende, sind eine zunehmende Anzahl von Gabeln von Git-Repositories (diskrete Variable, nicht kategorisch), wie Sie können siehe in diesem Datensatz . Ich...
Die Daten bestehen aus optischen Spektren (Lichtintensität gegen Frequenz), die zu unterschiedlichen Zeiten aufgenommen wurden. Die Punkte wurden in einem regelmäßigen Raster in x (Zeit), y (Frequenz) erfasst. Um die zeitliche Entwicklung bei bestimmten Frequenzen zu analysieren (ein schneller...
Ich habe Probleme, eine Lösung für die Durchführung eines Post-hoc-Tests (Tukey HSD) nach einer ANOVA mit 2 Faktoren (beide innerhalb der Probanden) mit wiederholten Messungen in R zu finden. Für die ANOVA habe ich die aov-Funktion verwendet: summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)),...
Beim Umgang mit Daten mit Faktoren kann R verwendet werden, um die Mittelwerte für jede Gruppe mit der Funktion lm () zu berechnen. Dies gibt auch die Standardfehler für die geschätzten Mittelwerte an. Dieser Standardfehler unterscheidet sich jedoch von dem, was ich aus einer manuellen Berechnung...
Ich verwende das lme4Paket in R, um eine logistische Modellierung mit gemischten Effekten durchzuführen. Mein Verständnis war, dass die Summe aller zufälligen Effekte Null sein sollte. Wenn ich mit Spielzeug lineare gemischte Modelle mache lmer, sind die zufälligen Effekte normalerweise < was...
Ich habe Daten, die ungefähr so aussehen: ID Status 01 A 02 G 03 E ... ... 100 G Sie haben die Idee, denke ich. Ich habe diese Daten aus zwei verschiedenen Populationen (Kohorten) und möchte die Verteilung der Statusvariablen in einer Population mit der Verteilung in einer anderen vergleichen....
Ich habe eine Frage, wie ein Statistiker normalerweise eine Anova-Ausgabe interpretieren würde. Angenommen, ich habe eine Anova-Ausgabe von R. > summary(fitted_data) Call: lm(formula = V1 ~ V2) Residuals: Min 1Q Median 3Q Max -2.74004 -0.33827 0.04062 0.44064 1.22737 Coefficients:...
Der Versuch, die Ergebnisse des kürzlich veröffentlichten Artikels zu replizieren, Aghion, Philippe, John Van Reenen und Luigi Zingales. 2013. "Innovation und institutionelles Eigentum." American Economic Review, 103 (1): 277-304. (Daten- und Datencode finden Sie unter
Wir wissen, dass wenn ein Schätzer ein unverzerrter Schätzer von Theta ist und seine Varianz gegen 0 tendiert, während n gegen unendlich tendiert, er ein konsistenter Schätzer für Theta ist. Dies ist jedoch eine ausreichende und keine notwendige Bedingung. Ich suche ein Beispiel für einen Schätzer,...
Ich versuche, ein logarithmisch lineares Modell an eine große Anzahl von Variablen aus Umfragedaten anzupassen. Es gibt einige Gründe, warum es möglicherweise vorzuziehen ist, stattdessen logistische Regressionen an diese Daten anzupassen. Mehrere Behörden schlagen vor, dass diese gleichwertig...
Ich weiß, dass R für die Analyse großer Datenmengen nicht besonders hilfreich ist, da R alle Daten in den Speicher lädt, während SAS eine sequentielle Analyse durchführt. Es gibt jedoch Pakete wie bigmemory, mit denen Benutzer die Analyse großer Datenmengen (statistische Analysen) in R effizienter...
Unter Verwendung des hervorragenden Prognosepakets von Rob Hyndman stieß ich auf die Notwendigkeit, nicht nur Vorhersageintervalle zu haben, sondern auch eine Reihe zukünftiger Pfade zu simulieren, wenn man frühere Beobachtungen einer Zeitreihe mit komplexen Saisonalitäten berücksichtigt. Es gibt...
Ich habe einen variablen Satz von Antworten, die als Intervall ausgedrückt werden, wie im folgenden Beispiel. > head(left) [1] 860 516 430 1118 860 602 > head(right) [1] 946 602 516 1204 946 688 Dabei ist links die Untergrenze und rechts die Obergrenze der Antwort. Ich möchte die Parameter...
Ich versuche diese hintere Verteilung zu
Bei der Vorhersage von Zeitreihen mit verschiedenen Modellen wie AR, MA, ARMA usw. konzentrieren wir uns normalerweise auf die Modellierung der Daten im Zeitwechsel. Wenn wir jedoch zwei Zeitreihen haben, bei denen der Pearson-Korrelationskoeffizient zeigt, dass sie stark korreliert sind, ist es...
Ich muss die Dirichlet-CDF berechnen , kann aber nur Implementierungen des PDF finden . Kennt ihr eine Bibliothek (vorzugsweise in R), die sie