maotp_ds10

Groupe an3

Lecture du fichier

Télécharger le fichier "www.douillet.info/~douillet/maotp/maotp_ds10/dat_maotp_ds10.txt"

Déclarer le nom de la copie locale du fichier

> fich:= "F:/docs/Ensait/maotp/maotp_ds10/dat_maotp_ds10b.txt";

On lit tous les nombres

> lire:= table(): tmp:= 'tmp':
hdl:= fopen(fich, READ);
for i to 5000 while tmp <> [ ] do
tmp:= fscanf(hdl,"%a");
lire[i]:= op (tmp);
od:
fclose(fich); n:= i-2 ; i:='i':

On en fait une liste

> li2:= [seq(lire[k], k=1..n)]:

Paramètres de dispersion

> macro(moy= stats[describe, mean], var= stats[describe, variance],
ran= stats[describe, range], nmb= stats[describe, count] );

> m:= moy(li2);
s, s2:= sqrt(var(li2)), var(li2);
vals:= ran(li2);
ASSERT(n=nmb(li2));

Visualisation

> macro(histo = xhisto, talto=stats[transform, tallyinto] );

Détermination du nombre de barres

> histo(li2, area=1); pl1a:= %:

> histo(li2, area=1, numbars=100); pl1b:= %:

> nbx:= 40: histo(li2, area=1, numbars=nbx); pl1:= %:

> xima:= proc(pl); op(1,pl): convert(%,list): map(op,%): map2(op,2,%); (max@op)(%); end:

> xima(pl1a), xima(pl1b), xima(pl1); yx:= %[3];

> pl2:= plot([[m,y,y=0..yx], [m-s,y,y=0..yx], [m+s,y,y=0..yx]], color=red, linestyle=[24,8,8]):

> display(pl1, pl2);

Modèle lognormal

Paramètres loi lognormale

> M:='M': K:='K': m=M*sqrt(K), (s/m)^2 = K-1; solve({%}); assign(%);

> M, K;

> pl5:= plot(norlaw(ln(M), sqrt(ln(K)), ln(x))/x, x=0..rhs(vals), color=blue, linestyle=16):

> top:= exp(ln(M)+2.5*sqrt(ln(K)));

> cut:= proc(li, num);
stats[transform, split[num]](stats[transform, statsort](li)):
map(z-> Weight(ran(z), nmb(z)), %);
end;

> tal2:= cut(li2,20);

> pl3:= histo(tal2, area=1, color=yellow):

> display(pl1,pl2,pl5, view=[0..top, DEFAULT]); display(pl3,pl2,pl5, view=[0..top, DEFAULT]);

>

Loi "en masse"

> newcode:= unapply(Weight(z,z/m), z);

> li3:= map(newcode, sort(li2)):
nmb(li3);
m_:= moy(li3); v_:= var(li3): s_:= sqrt(%);

> pl21:=histo(li3, area=1, numbars=nbx): yx:= xima(%):

> pl22:= plot([[m_,y,y=0..yx], [m_-s_,y,y=0..yx], [m_+s_,y,y=0..yx]], color=red, linestyle=[24,8,8]):

> tal3:= cut(li3,20);

> pl23:= histo(tal3, area=1, color=yellow):

> jx:=nbx/2: [0, seq(exp(INorlaw(ln(K*M), sqrt(ln(K)), j/jx)), j=1..jx-1), rhs(vals)*1.1]:
les_interv:= seq(%[j]..%[j+1], j=1..nops(%)-1);
tal3th:= talto(li3, [les_interv]):