Rev 7209 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 7209 | Rev 13586 | ||
---|---|---|---|
Line 1... | Line 1... | ||
1 | !if $wims_read_parm!=slib_header |
1 | !if $wims_read_parm!=slib_header |
2 | !goto proc |
2 | !goto proc |
3 | !endif |
3 | !endif |
4 | 4 | ||
5 | slib_parms=1\ |
5 | slib_parms=1\ |
6 | ,[data] or [data],[weight] |
6 | ,[data] or [data],[weight] |
7 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
7 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
Line 21... | Line 21... | ||
21 | !readproc slib/stat/dataproc $wims_read_parm |
21 | !readproc slib/stat/dataproc $wims_read_parm |
22 | slib_cnt=!itemcnt $slib_data |
22 | slib_cnt=!itemcnt $slib_data |
23 | !if $slib_cnt=0 |
23 | !if $slib_cnt=0 |
24 | slib_out=0 |
24 | slib_out=0 |
25 | !else |
25 | !else |
26 | !if $slib_weight=$empty |
26 | !if $slib_weight=$empty |
27 | slib_data1=!sum x for x in $slib_data |
27 | slib_data1=!sum x for x in $slib_data |
28 | slib_mean=$[($slib_data1)/$slib_cnt] |
28 | slib_mean=$[($slib_data1)/$slib_cnt] |
29 | slib_dev=!sum (x-($slib_mean))^2 for x in $slib_data |
29 | slib_dev=!sum (x-($slib_mean))^2 for x in $slib_data |
30 | slib_out=$[$slib_dev/$slib_cnt],$[$slib_dev/($slib_cnt-1)] |
30 | slib_out=$[$slib_dev/$slib_cnt],$[$slib_dev/($slib_cnt-1)] |
31 | slib_out =!trim $slib_out |
31 | slib_out =!trim $slib_out |
32 | !exit |
32 | !exit |
33 | !else |
33 | !else |
34 | slib_out=!exec pari { slib_mean=([$slib_data]*[$slib_weight]~)/$slib_tw;\ |
34 | slib_out=!exec pari { slib_mean=([$slib_data]*[$slib_weight]~)/$slib_tw;\ |
35 | slib_dev=sum(x=1,$slib_cnt,[$slib_weight][x]*[$slib_data][x]^2)-$slib_tw*(slib_mean)^2;\ |
35 | slib_dev=sum(x=1,$slib_cnt,[$slib_weight][x]*[$slib_data][x]^2)-$slib_tw*(slib_mean)^2;\ |
36 | print([1.*slib_dev/$slib_tw,1.*slib_dev/($slib_tw-1)])} |
36 | print([1.*slib_dev/$slib_tw,1.*slib_dev/($slib_tw-1)])} |
37 | slib_out =!trim $slib_out |
37 | slib_out =!trim $slib_out |
38 | !endif weight |
38 | !endif weight |
39 | 39 | ||
40 | slib_out =!trim $slib_out |
40 | slib_out =!trim $slib_out |
41 | !endif |
41 | !endif |
42 | slib_out=!trim $slib_out |
42 | slib_out=!trim $slib_out |