#Estimation des parametres a et lambda dans le modele Gamma. n=100 N=1000 a=15 #loi de Erlang lambda=10 X=rgamma(n*N,a,rate=lambda) X=matrix(X,N,n) #Echantillons dans matrice moy=apply(X,1,mean) #moyenne empirique de chaque echantillon varemp=apply(X,1,var) #variance empirique de chaque echantillon estima=moy^2/varemp #estimation du parametre a sur chaque echantillon estiml=moy/varemp #estimation du parametre lambda sur chaque echantillon boxplot(data.frame(estima)) #boxplot abline(h=a,col="red") #Vrai valeur x11() boxplot(data.frame(estiml)) #boxplot abline(h=lambda,col="red") #Vrai valeur