histogramme graphique de la distribution de Cauchy dans R

Question

Je suis en train de tracer un histogramme de la distribution de Cauchy dans R en utilisant le code suivant:

X = rcauchy(10^5) 
hist(X) 

et peu importe quelles sont les options que j'essayer dans la fonction hist(), je peux ne vois jamais plus de deux barres sur mon histogramme (essentiellement une pour les valeurs négatives et une pour les valeurs positives).

Cela fonctionne bien, cependant, quand j'utilise la distribution normale (ou autres).

Answer 1

Ceci résulte des propriétés de la distribution.

La plupart des valeurs sont relativement proches de zéro, mais de très grandes valeurs absolues sont beaucoup plus probables que pour la distribution normale. Il y a environ 1 valeurs% avec une valeur absolue supérieure à 50, et de 0,1% supérieur à 500.

Essayez tracer seulement une partie des valeurs:

hist(X[abs(X)<1]) 
hist(X[abs(X)<5]) 
hist(X[abs(X)<50]) 
hist(X) 

Vous pouvez également regarder la fonction de distribution cumulative:

plot(ecdf(X)) 

et vérifier les boxplot:

boxplot(X)