Qu'est-ce qu'une région à plus haute densité (HDR)?


23

Dans l'inférence statistique , le problème 9.6b, une "région de plus haute densité (HDR)" est mentionné. Cependant, je n'ai pas trouvé la définition de ce terme dans le livre.

Un terme similaire est la plus haute densité postérieure (HPD). Mais cela ne rentre pas dans ce contexte, car 9.6b ne mentionne rien sur un a priori. Et dans la solution suggérée , il dit seulement que "évidemment est un HDR".c(y)

Ou le HDR est-il une région contenant le (s) mode (s) d'un pdf?

Qu'est-ce qu'une région à plus haute densité (HDR)?


Oui. La page amazon est le livre (page d'achat). Le pdf est la solution aux problèmes du livre.
user3813057

Réponses:


33

Je recommande l'article de Rob Hyndman de 1996 «Calcul et représentation graphique des régions à plus haute densité» dans The American Statistician . Voici la définition du HDR, tirée de cet article:

Soit soit la fonction de densité d'une variable aléatoire . Alors le HDR est le sous-ensemble de l'espace échantillon de tel que où est la plus grande constante telle que F(X)X100(1-α)%R(Fα)X

R(Fα)={X:F(X)Fα},
Fα
P(XR(Fα))1-α.

La figure 1 de cet article illustre la différence entre le 75% HDR (donc ) et diverses autres régions de probabilité de 75% pour un mélange de deux normales ( est le -ème quantile, la moyenne et le écart type de la densité):α=0,25cqqμσ

HDR

L'idée dans une dimension est de prendre une ligne horizontale et de la déplacer vers le haut (jusqu'à ) jusqu'à ce que la zone au-dessus et en dessous de la densité soit . Alors le HDR est la projection sur l' axe de cette zone.y=Fα1-αRαX

Bien sûr, tout cela fonctionne avec n'importe quelle densité, qu'elle soit postérieure bayésienne ou autre.

Voici un lien vers le code R, qui est le hdrcdepackage (et vers l'article sur JSTOR).


4

Une densité postérieure [intervalle] la plus élevée est fondamentalement l'intervalle le plus court sur une densité postérieure pour un certain niveau de confiance donné. Une région de densité la plus élevée est probablement la même idée appliquée à toute densité arbitraire, donc pas nécessairement une distribution postérieure.

1-αq1-α/2+cqα/2-c

F()unebF(une)=F(b)

En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.