2010-06-11 4 views
8

Ich habe Daten dieser Form:Statistik Frage: Kernel Smoothing in R

x y 
1 0.19 
2 0.26 
3 0.40 
4 0.58 
5 0.59 
6 1.24 
7 0.68 
8 0.60 
9 1.12 
10 0.80 
11 1.20 
12 1.17 
13 0.39 

Ich Plotten derzeit einen Kernel-geglätteten Dichteschätzung der x gegen y mit diesem Code:

smoothed = ksmooth(d$resi, d$score, bandwidth = 6) 
    plot(smoothed) 

Ich möchte einfach eine grafische Darstellung der x gegen geglätteten (y) Werte, die ist ##

Die documentation for ksmooth schlägt jedoch vor, dass dies nicht die beste Kernel-Glättung Paket verfügbar ist :

Diese Funktion ist rein für die Kompatibilität mit S umgesetzt, obwohl es bei weitem nicht so langsam wie die S Funktion ist. Bessere Kernelglätter sind in anderen Paketen verfügbar.

Welche anderen Kernelglätter sind besser und wo können diese Glätter gefunden werden?

Antwort

11

Wenn Sie "einfach eine Darstellung des x gegen geglätteten (y)" wollen, dann empfehle ich loess in Paket stats - es ist einfach, schnell und effektiv. Wenn Sie stattdessen eine Regression basierend auf der Kernel-Glättung wünschen, können Sie locpoly im Paket KernSmooth oder npreg im Paket np ausprobieren.

+0

"Bessere Kernelglätter sind in anderen Paketen verfügbar" ... in welchem ​​Sinne ist Locpoly besser? –