Distribuţia Binomială: Modelare Statistică, Optimizare Numerică, cu ...

Distribuţia Binomială: Modelare Statistică, Optimizare Numerică, cu Aplicaţii în Bioinformatică şi Biochimie
i := 0;
repeat
for j := 0 to 5 do t[i][j] := v_s[j]-1;
i := i + 1;
j := 0;
v_s[j] := v_s[j] + 1;
while(v_s[j]>=3)do begin
v_s[j+1] := v_s[j+1] + 1;
v_s[j] := v_s[j] - 3;
j := j + 1;
end;
until(i>=729);
end;
Intervalele de încredere OAB
Intervalele de încredere OAB se obţin din matricea de probabilitate pentru n şi α fixaţi şi
pentru X şi Y variind în domeniul de posibilităţi (0..n). Algoritmul constă în identificarea limitelor
intervalului de încredere respectând monotonia de variaţie a limitelor acestuia şi urmărind cea mai
mică abatere a sumei probabilităţilor exterioare intervalului ([XS,n] şi [0,XI]) faţă de valoarea
impusă (α).
÷ Formule matematice:
Xi
CIOAB(X,m) = (Xi,Xs) | ∑PB
( m,
X,
i)
+ ∑PB
( m,
X,
i)
− α = min.
,
i=
0
49(157)
m
i=
Xs
Xi(X-1) ≤ Xi(X); Xs(X-1) ≤ Xs(X)
÷ Algoritmi de calcul:
const
alpha=5.0; lmax=1000; fname='rez';gname='cie';
type
arr_data=array[0..lmax] of integer;
row_data=array[0..lmax] of extended;
prowdata=^row_data;
mat_data=array[0..lmax] of prowdata;
var
p_m:mat_data;err:row_data;cil,ciu:arr_data;
var
m,i,j:integer;
begin
mat_init(p_m);
write('m=');readln(m);
mat_pro(p_m,m);
cie_init(cil,ciu,err,m);
ci_comp(p_m,cil,ciu,err,m);
cie_afi(cil,ciu,err,m);
mat_afi(p_m,m);
mat_done(p_m);
end.
procedure cie_init(var cil,ciu:arr_data;var
err:row_data;m:integer);
var
i:integer;