Beim Programmieren in R habe ich das Multicore- Paket einige Male verwendet. Ich habe jedoch noch nie eine Aussage darüber gesehen, wie es mit Zufallszahlen umgeht. Wenn ich openMP mit C verwende, achte ich darauf, ein korrektes paralleles RNG zu verwenden, aber mit R gehe ich davon aus, dass etwas...
Geschlossen. Diese Frage ist nicht zum Thema . Derzeit werden keine Antworten akzeptiert. Geschlossen vor 4 Jahren . Verschlossen . Diese Frage und ihre Antworten sind gesperrt, da die Frage nicht zum Thema gehört, aber von historischer Bedeutung ist. Derzeit...
Bei zwei Arrays x und y, die beide die Länge n haben, passe ich ein Modell y = a + b * x an und möchte ein 95% -Konfidenzintervall für die Steigung berechnen. Dies ist (b - delta, b + delta), wobei b in üblicher Weise und gefunden wird delta = qt(0.975,df=n-2)*se.slope und se.slope ist der...
Ich habe mit glm in R ein logistisches Regressionsmodell erstellt. Ich habe zwei unabhängige Variablen. Wie kann ich die Entscheidungsgrenze meines Modells im Streudiagramm der beiden Variablen darstellen? Wie kann ich beispielsweise eine Abbildung wie die folgende zeichnen:
Ich sammle zurzeit Daten für ein Experiment zu psychosozialen Merkmalen, die mit dem Erleben von Schmerz verbunden sind. Im Rahmen dessen sammle ich von meinen Teilnehmern elektronisch GSR- und BP-Messungen sowie verschiedene Selbstberichte und implizite Maßnahmen. Ich habe einen psychologischen...
Ich benutze GNU R auf einem Ubuntu-Lucid-PC mit 4 CPUs. Um alle 4 CPUs nutzen zu können, habe ich das Paket "r-cran-multicore" installiert. Da das Handbuch des Pakets keine praktischen Beispiele enthält, die ich verstehe, benötige ich Ratschläge zur Optimierung meines Skripts, um alle 4 CPUs nutzen...
Wie würden Sie stündliche Mittelwerte für mehrere Datenspalten für einen täglichen Zeitraum erhalten und die Ergebnisse für zwölf "Hosts" in derselben Grafik anzeigen? Das heißt, ich möchte für Daten im Wert von mehreren Wochen grafisch darstellen, wie ein 24-Stunden-Zeitraum aussieht. Das...
Ich versuche, eine "Änderungspunkt" -Analyse oder eine mehrphasige Regression mit nls()in R zu implementieren . Hier sind einige gefälschte Daten, die ich gemacht habe . Die Formel, die ich verwenden möchte, um die Daten anzupassen, lautet: y= β0+ β1x + β2max ( 0 , x - δ)y=β0+β1x+β2max(0,x-δ)y =...
Ich habe mit einigen Unit-Root-Tests in R herumgespielt und bin mir nicht ganz sicher, was ich mit dem Parameter k lag anfangen soll. Ich habe den erweiterten Dickey-Fuller-Test und den Philipps-Perron-Test aus dem tseries- Paket verwendet. Offensichtlich hängt der voreingestellte Parameter (für )...
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 im vergangenen Jahr . Könnte sich jemand einen R- Code einfallen lassen, um eine Ellipse...
Ich würde gerne Wege in R finden, um ein lineares Modell effizient zu aktualisieren, wenn eine Beobachtung oder ein Prädiktor hinzugefügt wird. biglm kann beim Hinzufügen von Beobachtungen aktualisiert werden, aber meine Daten sind klein genug, um sich im Speicher zu befinden (obwohl ich eine große...
Angenommen, ich habe ein einfaches Experiment mit 2 × 2 Faktoren, an dem ich eine ANOVA durchführen möchte. So zum Beispiel: d <- data.frame(a=factor(sample(c('a1','a2'), 100, rep=T)), b=factor(sample(c('b1','b2'), 100, rep=T))); d$y <- as.numeric(d$a)*rnorm(100, mean=.75, sd=1) +...
Wie kann man diesem Fehler entgegenwirken, wenn Sie bei einem GLM den Fehler "Wegen Singularitäten nicht definiert" in der anova-Ausgabe erhalten? Einige haben vermutet, dass es an der Kollinearität zwischen den Kovariaten liegt oder dass eine der Ebenen im Datensatz nicht vorhanden ist (siehe:...
Ich versuche, die ' Dichte' -Funktion in R zu verwenden, um Kernel-Dichteschätzungen durchzuführen. Ich habe einige Schwierigkeiten, die Ergebnisse zu interpretieren und verschiedene Datensätze zu vergleichen, da die Fläche unter der Kurve nicht unbedingt 1 zu sein scheint. Für jede...
Verschlossen . Diese Frage und ihre Antworten sind gesperrt, da die Frage nicht zum Thema gehört, aber von historischer Bedeutung ist. Derzeit werden keine neuen Antworten oder Interaktionen akzeptiert. Ich habe die Grafik wie folgt: R-Code für die Generierung ist: DF <- data.frame(date =...
Ich habe ein klassisches lineares Modell mit 5 möglichen Regressoren. Sie sind nicht miteinander korreliert und weisen eine relativ geringe Korrelation mit der Antwort auf. Ich bin zu einem Modell gekommen, bei dem 3 der Regressoren signifikante Koeffizienten für ihre t-Statistik haben (p...
Es ist schon eine Weile her, dass ich mir das Teilen von Bäumen angesehen habe. Als ich das letzte Mal so etwas gemacht habe, mag ich Party in R (erstellt von Hothorn). Die Idee der bedingten Folgerung durch Stichproben ist für mich sinnvoll. Aber rpart hatte auch Anklang. In der aktuellen...
Es gab mehrere gute Fragen und Antworten zu Einführungsbüchern oder Lernansätzen, z . B. hier und hier . Aber ich habe ein etwas anderes Problem - die beste Möglichkeit, eine einstündige Sitzung (oder mehrere solcher Sitzungen) in einem Computerraum durchzuführen, damit die Leute mit R beginnen und...
Ich analysiere Daten aus einem unausgeglichenen faktoriellen Experiment mit SASund R. Beide SASund Rliefern eine ähnliche Quadratsumme vom Typ I, aber ihre Quadratsumme vom Typ III unterscheidet sich voneinander. Im Folgenden sind SASund RCodes und Ausgänge. DATA ASD; INPUT Y T B; DATALINES; 20 1 1...
CARET verwendet automatisch ein vorab festgelegtes Abstimmungsraster, um verschiedene Modelle zu erstellen, bevor ein endgültiges Modell ausgewählt wird, und trainiert dann das endgültige Modell mit den vollständigen Trainingsdaten. Ich kann mein eigenes Stimmraster mit nur einer Kombination von...