Ich habe eine Dichte Objekt wie folgt erstellt DD:aus einer Dichte Objekt mit R
x1 <- rnorm(1000)
x2 <- rnorm(1000, 3, 2)
x <- rbind(x1, x2)
dd <- density(x)
plot(dd)
die eben diese nicht-Gauß-Verteilung erzeugt:
alt text http://www.cerebralmastication.com/wp-content/uploads/2009/09/nongaus.png
würde ich Letztendlich mag es Zufallsabweichungen von dieser Verteilung zu bekommen, ähnlich wie rnorm von einer normalen Verteilung abweicht.
Die Art, wie ich dies zu knacken versuchen ist die CDF meiner Kernel zu bekommen und dann bekommt es mir die Veränderlichen zu sagen, wenn ich ihm eine kumulative Wahrscheinlichkeit (inverse CDF) übergebe. Auf diese Weise kann ich einen Vektor von uniformen Zufallsvariablen in Draws aus der Dichte drehen.
Es scheint so, als ob das, was ich versuche, etwas Grundlegendes sein sollte, was andere vor mir getan haben. Gibt es dafür einen einfachen Weg oder eine einfache Funktion? Ich hasse es, das Rad neu zu erfinden.
FWIW Ich fand this R Help article, aber ich kann nicht was sie tun, und die endgültige Ausgabe scheint nicht zu produzieren, was ich danach bin. Aber es könnte ein Schritt auf dem Weg sein, den ich einfach nicht verstehe.
Ich habe überlegt, nur mit einem Johnson distribution from the suppdists package gehen, aber Johnson wird mir nicht den schönen bimodalen Buckel geben, die meine Daten haben.
mehr eine Statistik Frage als Programmierung ... –
weiß, dass ich die Statistiken. Ich möchte die Statistikmethode in einer bestimmten Sprache implementieren. Das ist Programmieren. –