2017-09-12 1 views
3

Hallo, alle. Also ich weiß, dass das ziemlich einfach sein muss, aber was ist die richtige, genaue Art, das pdf von einigen Beispieldaten zu zeichnen, die Sie wissen, kommt von etwas Pop. Verteilung, wie wenn Sie es unter Verwendung rnorm() oder rexp() erzeugten?Die richtige Art, PDF einer Stichprobe von Daten zu plotten?

Der Grund, warum ich frage ist, weil ich eine Menge Leute density() nutzen wissen, und geben Sie dann, dass in plot(), aber die density() Funktion scheint zu beliebig genau zu sein; Zum Beispiel ist es ungenau, wenn es einen negativen Wert für Daten, die aus der Exponentialverteilung stammen, approximiert, die neg nicht besitzt. Werte.

Also könnte mir jemand eine feinere Methode empfehlen, Plot-Beispiele zu erstellen?

+1

eine ähnliche Frage (oder zumindest ähnliche), wurde gebeten, auf [stats.se] nehmen würde - mein [answer] (https://stats.stackexchange.com/a/71291/1390) behandelt das Bit über die unplausiblen negativen Werte für eine Variable. –

Antwort

1

Die Funktion density führt eine Kerneldichteschätzung (KDE) durch. Um die beste KDE für Ihr Dataset zu finden, sollten Sie die Bandbreite einstellen (Parameter bw). Hier ist ein Papier, das KDE und Bandbreite Auswahl diskutiert: http://www.stat.washington.edu/courses/stat527/s13/readings/Sheather_StatSci_2004.pdf

Oder für einen einfacheren Ansatz, können Sie verschiedene Bandbreite Methoden ausprobieren passieren zu bw: https://stat.ethz.ch/R-manual/R-devel/library/stats/html/bandwidth.html

Der aktuelle Standard „nrd0“, gibt es für historische Gründe. Ich finde, dass "ucv" und "bcv" besser für meine Datensätze funktioniert haben.

+1

Wenn die Datenwerte nahe bei 0 liegen, wird die Dichte immer noch in negative Werte umgewandelt, es sei denn, Sie haben tatsächlich eine sehr kleine Bandbreite. –

0

ggplot hilft bei negativen Werten, wenn sie nicht geeignet sind. Es kann auf folgende Weise verwendet werden:

ggplot(df, 
     aes(x=contVar, fill = "green")) + 
    geom_density(alpha=.3) 

ich auch einen Blick auf diesen Beitrag in cross validated

+0

Während ich dachte, dass die obige Antwort meine Frage auf eine technischere und gründlichere Weise ansprach, ist Ihr Link eine großartige Ressource und möchte Ihnen dafür danken. – Coolio2654

Verwandte Themen