Ich muss so genau wie möglich die Spitze der Kerndichteschätzung finden (Modalwert der kontinuierlichen Zufallsvariablen). Ich kann den ungefähren Wert finden:
%Vor% Aber bei der Berechnung von d$y ist die genaue Funktion bekannt. Wie kann ich den genauen Wert des Modus finden?
Hier sind zwei Funktionen für den Umgang mit Modi. Die dmode-Funktion findet den Modus mit der höchsten Spitze (dominate mode) und n.modes identifizieren die Anzahl der Modi.
%Vor% Wenn ich Ihre Frage verstehe, denke ich, dass Sie nur eine feinere Diskretisierung von x und y wollen. Um dies zu tun, können Sie den Wert von n in der Funktion density ändern (Standard ist n=512 ).
Vergleiche beispielsweise
%Vor%mit:
%Vor%Ich denke, Sie brauchen zwei Schritte, um zu archivieren, was Sie brauchen:
1) Finde den x-Achsenwert des KDE-Peaks
2) Erhalte den Desnity-Wert des Peaks
Also (wenn es Ihnen nichts ausmacht, ein Paket zu verwenden) würde eine Lösung, die das Paket hdrcde verwendet, so aussehen:
Bearbeiten: Sie könnten auch das
verwenden %Vor%wenn es genau genug für Sie ist!
Tags und Links r kernel-density