Para crear un gráfico de densidad tipo núcleo tendrás que estimar primero la densidad. Para ello puedes hacer uso de la función density
y luego pasar el objeto a la función plot
.
El argumento kernel
de la función density
usa por defecto un kernel gaussiano (kernel = "gaussian"
), pero existen otros tipos de kernel disponibles, como "rectangular"
, "triangular"
, "epanechnikov"
, "biweight"
, "cosine"
y "optcosine"
. La seleción dependerá de tus datos, pero el valor por defecto es el más recomendado.
Kernel rectangular
Kernel triangular
Kernel epanechnikov
Kernel biweight
Kernel cosine
El argumento bw
de la función density
permite cambiar el parámetro ventana (bandwidth). Puedes pasar un valor, una cadena de una regla de selección o una función. El valor por defecto es "nrd0"
(o bw.nrd0(.)
), que implementa el procedimiento conocido como rule-of-thumb. Otras opciones son:
Variación de la rule-of-thumb de Scott (1992)
"nrd"
o bw.nrd(.)
Validación cruzada insesgada
"ucv"
o bw.ucv(.)
Validación cruzada sesgada
"bcv"
o bw.bcv(.)
Métodos de Sheather & Jones (1991)
"SJ"
o bw.SJ(.)
El bandwidth se debe elegir con precaución. Una ventana demasiado pequeña creará una curva sobreajustada mientras que una ventana demasiado grande creará una curva demasiado suave.
También te puede interesar