Lesen mgcv::gamder Hilfeseite: Vertrauens- / Glaubwürdigkeitsintervalle sind für jede anhand eines angepassten Modells vorhergesagte Menge leicht verfügbar Allerdings kann ich keinen Weg finden, um tatsächlich einen zu bekommen. Ich dachte, ich predict.gamhätte einen type=confidenceund einen...
Ich muss die Dichtefunktion basierend auf einer Reihe von Beobachtungen mit dem Kernel-Dichteschätzer abschätzen. Basierend auf den gleichen Beobachtungen muss ich auch die erste und die zweite Ableitung der Dichte unter Verwendung der Ableitungen des Kerndichteschätzers schätzen. Die Bandbreite...
Ich versuche, ein QQ-Diagramm mit zwei Datensätzen von ungefähr 1,2 Millionen Punkten in R zu zeichnen (unter Verwendung von qqplot und Eingabe der Daten in ggplot2). Die Berechnung ist einfach genug, aber das resultierende Diagramm ist schmerzhaft langsam zu laden, weil es so viele Punkte gibt....
Ich habe eine Reihe von Daten aus zwei Proben (Kontrolle und Behandlung) mit jeweils mehreren tausend Werten, die in R einer Signifikanzprüfung unterzogen werden sollen. Theoretisch sollten die Werte kontinuierlich sein, aber aufgrund der Rundung durch die Messsoftware sind sie t und sie haben...
Ich habe kürzlich einen Kalman-Filter am einfachen Beispiel der Messung einer Partikelposition mit einer zufälligen Geschwindigkeit und Beschleunigung implementiert. Ich habe festgestellt, dass der Kalman-Filter gut funktioniert, habe mich dann aber gefragt, was der Unterschied zwischen diesem und...
Ich habe einen Vektor von Zahlen, den ich hier mit dput hochgeladen habe (... / code / MyData.Rdata). Ich möchte das bca ci bekommen, also habe ich diesen Code geschrieben: my.mean <- function(dat, idx){ return (mean(dat[idx], na.rm = TRUE)) } boot.out<-boot(data=my.data, statistic = my.mean,...
Ich habe randomForest verwendet, um 6 Verhaltensweisen von Tieren (z. B. Stehen, Gehen, Schwimmen usw.) anhand von 8 Variablen (unterschiedliche Körperhaltungen und Bewegungen) zu klassifizieren. Der MDSplot im randomForest-Paket gibt diese Ausgabe aus und es treten Probleme bei der...
Ich führe eine Kreuzvalidierung mit der Methode "Auslassen" durch. Ich habe eine binäre Antwort und benutze das Boot-Paket für R und die cv.glm-Funktion . Mein Problem ist, dass ich den Teil "Kosten" in dieser Funktion nicht vollständig verstehe. Nach meinem Verständnis ist dies die Funktion, die...
Geschlossen. Diese Frage ist nicht zum Thema . Derzeit werden keine Antworten akzeptiert. Möchten Sie diese Frage verbessern? Aktualisieren Sie die Frage so dass es beim Thema für Kreuz Validated. Geschlossen vor 2 Jahren . Ich habe einen ziemlich...
Ich arbeite derzeit an einem Regressionsmodell, bei dem ich nur kategoriale / Faktor-Variablen als unabhängige Variablen habe. Meine abhängige Variable ist ein logit transformiertes Verhältnis. Es ist ziemlich einfach, eine normale Regression in R auszuführen, da R automatisch weiß, wie Dummies...
Ich habe Daten erhalten x = c(21,34,6,47,10,49,23,32,12,16,29,49,28,8,57,9,31,10,21,26,31,52,21,8,18,5,18,26,27,26,32,2,59,58,19,14,16,9,23,28,34,70,69,54,39,9,21,54,26) y = c(47,76,33,78,62,78,33,64,83,67,61,85,46,53,55,71,59,41,82,56,39,89,31,43,29,55,
Die Daten, die ich habe, sind ein Regressionssteigungswert von y ~ time, ein Standardfehler, ein n-Wert und ein p-Wert für eine bestimmte Spezies in zwei verschiedenen Gebieten. Ich möchte prüfen, ob die Regressionssteigung für einen Bereich erheblich von der Regressionssteigung für den anderen...
Geschlossen. Diese Frage ist nicht zum Thema . Derzeit werden keine Antworten akzeptiert. Möchten Sie diese Frage verbessern? Aktualisieren Sie die Frage so dass es beim Thema für Kreuz Validated. Geschlossen vor 5 Monaten . Beim Testen der Annahme der...
Ich bin neu in der R-Sprache. Ich möchte wissen, wie man aus einem multiplen linearen Regressionsmodell simuliert, das alle vier Annahmen der Regression erfüllt. OK danke. Angenommen, ich möchte die Daten basierend auf diesem Datensatz
Ist die folgende Formel richtig, wenn ich den Standardfehler des Medians bei einer kleinen Stichprobe mit nicht normaler Verteilung messen möchte (ich verwende Python)? sigma=np.std(data) n=len(data)
Ich frage mich, ob jemand mich wissen lassen könnte, ob ich das Konfidenzintervall für die Differenz zwischen zwei Anteilen richtig berechnet habe. Die Stichprobengröße beträgt 34, von denen 19 Frauen und 15 Männer sind. Daher beträgt der Unterschied in den Anteilen 0,1176471. Ich berechne das 95%...
Ich suche ein KNN-Anrechnungspaket. Ich habe mir das Imputationspaket angesehen ( http://cran.r-project.org/web/packages/imputation/imputation.pdf) ) angesehen, aber aus irgendeinem Grund scheint die KNN-Impute-Funktion (auch wenn dem Beispiel aus der Beschreibung folge) nur zu funktionieren...
Ich habe einen unsymmetrischen Datensatz mit wiederholten Messungen zur Analyse und ich habe gelesen, dass die Art und Weise, wie die meisten statistischen Pakete mit ANOVA umgehen (dh die Summe der Quadrate des Typs III), falsch ist. Daher würde ich gerne ein gemischtes Effektmodell verwenden, um...
Ich habe eine Zeitreihe, für die ich eine Prognose erstellen möchte und für die ich das saisonale Modell ARIMA (0,0,0) (0,1,0) [12] (= fit2) verwendet habe. Es unterscheidet sich von dem, was R mit auto.arima vorgeschlagen hat (R berechnete ARIMA (0,1,1) (0,1,0) [12] wäre besser, ich nannte es...
Dies ist eine Folgefrage, aber auch eine andere Frage als meine vorherige . Ich habe auf Wikipedia gelesen, dass " ein median-unverzerrter Schätzer das Risiko in Bezug auf die von Laplace beobachtete absolute Abweichungsverlustfunktion minimiert ". Meine Monte-Carlo-Simulationsergebnisse stützen...