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.
Eric Marcon
Eric Marcon
Chercheur en écologie

Mes centres d’intérêts en recherche incluent l’écologie, l’économie et la programmaton avec R.