Ich habe einen einfachen Kernel Density Estimator in Java entwickelt, der auf ein paar Dutzend Punkten (vielleicht bis zu einhundert oder so) und einer Gaußschen Kernelfunktion basiert. Die Implementierung gibt mir zu jedem Zeitpunkt das PDF und CDF meiner Wahrscheinlichkeitsverteilung.
Ich möchte jetzt eine einfache Stichprobenmethode für diese KDE implementieren . Eine naheliegende Wahl wäre natürlich, aus den Punkten des KDE zu ziehen, aber ich möchte in der Lage sein, Punkte abzurufen, die sich geringfügig von denen im KDE unterscheiden.
Ich habe bisher keine Stichprobentechnik gefunden, die ich leicht implementieren könnte, um dieses Problem zu lösen (ohne auf externe Bibliotheken für die numerische Integration oder komplexe Berechnungen angewiesen zu sein). Irgendwelche Ratschläge? Ich habe keine besonders hohen Anforderungen an Präzision oder Effizienz. Mein Hauptanliegen ist es, eine Stichprobenfunktion zu haben, die funktioniert und einfach implementiert werden kann. Vielen Dank!
rnorm(n, sample(dx$x, n, prob = dx$y, replace = TRUE), dx$bw)
Wodx
wird diedensity
Funktion ausgegeben ? Das Argumentprob
muss angegeben werden, da Sie sonst eine einheitliche Stichprobe erstellen.Antworten:
Wie von Procrastinator erwähnt, gibt es eine einfache Möglichkeit, eine Stichprobe mit einem Kernel-Dichteschätzer zu erstellen:
quelle