Densité avec réflexion
L’estimation de la densité d’une distribution sous R avec la fonction density
pose des problèmes si ses valeurs sont bornées.
L’exemple suivant est un histogramme de 1000 tirages d’une loi uniforme entre 1 et 2.
L’estimation avec la fonction density
n’est pas nulle hors de l’intervalle \([1, 2]\).
Censurer l’estimation avec les arguments from
et to
revient simplement de supprimer les valeurs estimées hors de l’intervalle pour obtenir la courbe rouge.
L’estimation de la densité est donc incorrecte: elle ne somme pas à 1 et est sous-estimée près des bornes.
Le package GoFKernel (Pavia 2015) propose la fonction density.reflected
pour estimer correctement la densité précédente, en vert.
Suivre les liens pour accéder au tutoriel complet:
Référence
Pavia, Jose M. 2015. “Testing Goodness-of-Fit with the Kernel Density Estimator: GoFKernel.” Journal of Statistical Software 66 (Code Snippet 1). https://doi.org/10.18637/jss.v066.c01.