Rev 20 | Rev 4161 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
20 | reyssat | 1 | !if $wims_read_parm!=slib_header |
2 | !goto proc |
||
3 | !endif |
||
4 | slib_title=Binomial law |
||
5 | slib_parms=3\ |
||
6 | ,[data]\ |
||
7 | 5,n\ |
||
8 | 0.5,p |
||
9 | |||
10 | slib_author=Bernadette PERRIN-RIOU |
||
11 | slib_out=values of binomial law of parameters n and p ((0<p<1)) on data |
||
12 | slib_comment=data can be a list of integers |
||
13 | slib_example=[2,4,5],20,0.7 |
||
3265 | bpr | 14 | slib_require=pari |
20 | reyssat | 15 | !exit |
16 | |||
17 | :proc |
||
18 | |||
19 | !distribute item $wims_read_parm into slib_data,slib_n,slib_p |
||
20 | |||
21 | !default slib_n= 5 |
||
22 | !default slib_p=0.5 |
||
23 | |||
24 | !if $slib_p>=1 or $slib_p<=0 |
||
25 | slib_out=0 |
||
26 | !exit |
||
27 | !endif |
||
28 | slib_out= |
||
29 | slib_data=!declosing $slib_data |
||
30 | slib_cnt= !itemcnt $slib_data |
||
31 | slib_data1=!values rint(x) for x in $slib_data |
||
32 | |||
33 | slib_n=$[rint($slib_n)] |
||
34 | slib_p1=$[1-$slib_p ] |
||
35 | slib_out=!exec pari {slib_data=[$slib_data1];\ |
||
36 | print(vector($slib_cnt,i,binomial($slib_n,slib_data[i])\ |
||
37 | *$slib_p^(slib_data[i])*$slib_p1^($slib_n-slib_data[i])));} |
||
38 | |||
39 | slib_out=!trim $slib_out |
||
40 | !for slib_i=1 to $slib_cnt |
||
41 | slib_tmp=!item $slib_i of $slib_data |
||
42 | !if $[rint($slib_tmp)-$slib_tmp]<>0 |
||
43 | slib_out=!replace internal item number $slib_i by 0 in $slib_out |
||
44 | !endif |
||
45 | !next slib_i |
||
46 | slib_out=!trim $slib_out |