bonsoir, <div>  ci-joint un programme un peu plus simple obtenir l'intervalle de fluctuation exact et centré concernant la loi binomial au seuil de 0.95, (ici en Xcas) traduisible très simplement en Scilab , voir sur calculatrice, avec le tableau de conversion des langages joint aussi ( ressource n°120 sur le site planète maths de l'académie de Grenoble)</div>

<div><br></div><div>Cordialement,</div><div>     Georges Dubouloz</div><div><br><br><div class="gmail_quote">Le 30 janvier 2013 07:12, Roland Thiers <span dir="ltr"><<a href="mailto:roland.thiers@gmail.com" target="_blank">roland.thiers@gmail.com</a>></span> a écrit :<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Bonjour,<br>
je ne suis pas du tout expert pour programmer (en Scilab ou dans un autre langage).<br>
J'aimerais bien savoir si le petit programme ci-dessous (relatif à l'intervalle de fluctuation) peut être amélioré.<br>
En particulier pour déterminer a et b j'ai fait une boucle , ce qui est je crois déconseillé, y-a-t-il un autre moyen ?<br>
merci à celle ou celui qui aura la patience de regarder ça !<br>
Cordialement, Roland Thiers<br>
<br>
// fonction donnant, pour un échantillon de taille n, l'intervalle<br>
// de fluctuation au seuil de 95 %  de la  fréquence de réalisation de la<br>
// variable aléatoire  X suivant la loi binomiale B(n,p).<br>
function IntFluctu(n,p)<br>
    binc = zeros(1,n+1);// liste de n+1 zéros dans la variable binc<br>
    binc(1)=loi_binomiale(n,p,0) // proba de 0 réalisation.<br>
    for k = 1:n<br>
binc(k+1)=loi_binomiale(n,p,k)<u></u>+binc(k);<br>
    end<br>
    // détermination de a, plus petit entier tel que p(X<=a)>0,025<br>
a=0<br>
while binc(a+1)<=0.025 do a=a+1<br>
end<br>
    // détermination de b, plus petit entier tel que p(X<=b)>=0,975<br>
b=0<br>
while binc(b+1)<0.975 do b=b+1<br>
end<br>
// pour afficher la loi binomiale B(n,p) cumulée<br>
//    u=0:n;<br>
//afficher(["k ",string(u) ; "P(X<=k)",string(binc)])<br>
  afficher("l''intervalle de fluctuation au seuil de 95% est ["+string(a/n)+" ; "+string(b/n)+"]")<br>
endfunction<br>
______________________________<u></u>_________________<br>
enseignement mailing list<br>
<a href="mailto:enseignement@lists.scilab.org" target="_blank">enseignement@lists.scilab.org</a><br>
<a href="http://lists.scilab.org/mailman/listinfo/enseignement" target="_blank">http://lists.scilab.org/<u></u>mailman/listinfo/enseignement</a><br>
</blockquote></div><br></div>