Les deux révisions précédentesRévision précédenteProchaine révision | Révision précédenteProchaine révisionLes deux révisions suivantes |
math:2:estimation_intervalle_confiance [2015/03/03 11:27] – [Estimation par intervalle de confiance de l'espérance d'une loi] Alain Guichet | math:2:estimation_intervalle_confiance [2016/02/15 21:55] – Alain Guichet |
---|
^ **[[:math:2:index#chapitre_19|Estimation > ]]** | [[:math:2:introduction_estimation|Intro]] | [[:math:2:estimation_parametre_bernoulli|Ex 1]] | [[:math:2:estimation_amplitude|Ex 2]] | [[:math:2:problematique_estimation|Problématique]] | [[:math:2:estimation_ponctuelle|Estimation ponctuelle]] | [[:math:2:estimation_intervalle_confiance|Estimation intervalle]] | | ^ **[[:math:2:index#estimation|Estimation > ]]** | [[:math:2:introduction_estimation|Intro]] | [[:math:2:estimation_parametre_bernoulli|Ex 1]] | [[:math:2:estimation_amplitude|Ex 2]] | [[:math:2:problematique_estimation|Problématique]] | [[:math:2:estimation_ponctuelle|Estimation ponctuelle]] | [[:math:2:estimation_intervalle_confiance|Estimation intervalle]] | |
| |
| |
On teste le théorème à l'aide du programme suivant : | On teste le théorème à l'aide du programme suivant : |
| |
rand("seed",getdate("s")) // initialisation du hasard | <code scilab> |
p=grand(1,1,"unf",0,1) | rand("seed",getdate("s")) // initialisation du hasard |
n=1000 ; alpha=0.05 ; t=1.96 ; m=10000 | p=grand(1,1,"unf",0,1) |
IBT=0 ; TLC1=0 ; TLC2=0 | n=1000 ; alpha=0.05 ; t=1.96 ; m=10000 |
for k=1:m | IBT=0 ; TLC1=0 ; TLC2=0 |
Xn=mean(grand(1,n,"bin",1,p)) | for k=1:m |
if abs(Xn-p)<t/sqrt(alpha*n)/2 then IBT=IBT+1 | Xn=mean(grand(1,n,"bin",1,p)) |
end | if abs(Xn-p)<t/sqrt(alpha*n)/2 then IBT=IBT+1 |
if abs(Xn-p)<t*sqrt(Xn*(1-Xn)/n) then TLC1=TLC1+1 | end |
end | if abs(Xn-p)<t*sqrt(Xn*(1-Xn)/n) then TLC1=TLC1+1 |
if abs(Xn-p)<t/sqrt(n)/2 then TLC2=TLC2+1 | end |
end | if abs(Xn-p)<t/sqrt(n)/2 then TLC2=TLC2+1 |
end | end |
disp("Proportion d''intervalles de type IBT contenant p : ") ; disp(100*IBT/m) | end |
disp("Proportion d''intervalles de type TLC1 contenant p : ") ; disp(100*TLC1/m) | disp("Proportion d''intervalles de type IBT contenant p : ") ; disp(100*IBT/m) |
disp("Proportion d''intervalles de type TLC2 contenant p : ") ; disp(100*TLC2/m) | disp("Proportion d''intervalles de type TLC1 contenant p : ") ; disp(100*TLC1/m) |
disp("p = ") ; disp(p) | disp("Proportion d''intervalles de type TLC2 contenant p : ") ; disp(100*TLC2/m) |
| disp("p = ") ; disp(p) |
| </code> |
| |
Tableau récapitulatif de résultats pour $n=1000$ et $\alpha=0.05$ (10000 répétitions) : | Tableau récapitulatif de résultats pour $n=1000$ et $\alpha=0.05$ (10000 répétitions) : |
| |
| |
^ **[[:math:2:index#chapitre_19|Estimation > ]]** | [[:math:2:introduction_estimation|Intro]] | [[:math:2:estimation_parametre_bernoulli|Ex 1]] | [[:math:2:estimation_amplitude|Ex 2]] | [[:math:2:problematique_estimation|Problématique]] | [[:math:2:estimation_ponctuelle|Estimation ponctuelle]] | [[:math:2:estimation_intervalle_confiance|Estimation intervalle]] | | ^ **[[:math:2:index#estimation|Estimation > ]]** | [[:math:2:introduction_estimation|Intro]] | [[:math:2:estimation_parametre_bernoulli|Ex 1]] | [[:math:2:estimation_amplitude|Ex 2]] | [[:math:2:problematique_estimation|Problématique]] | [[:math:2:estimation_ponctuelle|Estimation ponctuelle]] | [[:math:2:estimation_intervalle_confiance|Estimation intervalle]] | |