On considère la famille de lois $\mathcal{U}([a,b])$ dont le paramètre $\theta=(a,b)\in\R^{2}$ est inconnu. On souhaiterait obtenir une estimation de l'étendue $b-a$ de l'intervalle $[a,b]$. Pour cela, on considère un échantillon $(X_{1},\dots,X_{n})$ mutuellement indépendant et identiquement distribué selon la loi $\mathcal{U}([a,b])$ dont on possède une observation $(X_{1}(\omega),\dots,X_{n}(\omega))$. On pose alors :
$$\ds E_{n}(\omega)=\sup\{X_{1}(\omega),\dots,X_{n}(\omega)\}-\inf\{X_{1}(\omega),\dots,X_{n}(\omega)\}$$Ce réel est appelé étendue empirique. On va alors justifier que la suite de variables aléatoires $(E_{n})_{n\geqslant1}$ est un estimateur de $b-a$.
Simulation
On choisit les réels $a$ et $b$ au hasard dans, respectivement, $[0,1[$ et $[1,2[$. Ainsi, $b-a>0$. On détermine ensuite un 2000-échantillon de la loi $\mathcal{U}([a,b])$ puis on calcule l'étendue empirique (biaisée), l'étendue empirique corrigée (sans biais) et on affiche ensuite l'étendue réelle $b-a$.
n=2000 ; a=rand() ; b=1+rand() Xn=grand(1,n,"unf",a,b) Ee=max(Xn)-min(Xn) ; Eesb=Ee*(n+1)/(n-1) ; Er=b-a disp("Etendue estimee : ") ; disp(Ee) disp("Etendue estimee sans biais : ") ; disp(Eesb) disp("Etendue reelle : ") ; disp(Er)