Rev 4349 | Details | Compare with Previous | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
20 | reyssat | 1 | !if $wims_read_parm!=slib_header |
13586 | bpr | 2 | !goto proc |
20 | reyssat | 3 | !endif |
4161 | bpr | 4 | |
20 | reyssat | 5 | slib_parms=1\ |
6 | ,[data] or [data],[weight] |
||
4330 | bpr | 7 | slib_author=Adeline, Grelot |
4161 | bpr | 8 | |
20 | reyssat | 9 | slib_example=[1,2,5]\ |
10 | [1,2,3;2,4,5]\ |
||
11 | [1,3,3],[2,4,5] |
||
3265 | bpr | 12 | slib_require=pari |
20 | reyssat | 13 | !exit |
14 | |||
15 | :proc |
||
16 | slib_out= |
||
17 | |||
18 | !readproc slib/stat/dataproc $wims_read_parm |
||
19 | !if $slib_cnt<2 |
||
20 | slib_out=0,0 |
||
21 | !exit |
||
22 | !else |
||
23 | !if $slib_weight =$empty |
||
24 | slib_data1=!sum x for x in $slib_data |
||
25 | slib_mean=$[($slib_data1)/$slib_cnt] |
||
26 | slib_dev=!sum (x-($slib_mean))^2 for x in $slib_data |
||
27 | slib_out=$[sqrt($slib_dev/$slib_cnt)],$[sqrt($slib_dev/($slib_cnt-1))] |
||
28 | slib_out =!trim $slib_out |
||
29 | !exit |
||
30 | !else |
||
4349 | bpr | 31 | slib_out=!exec pari slib_mean=([$slib_data]*[$slib_weight]~)/$slib_tw;\ |
13586 | bpr | 32 | slib_dev=sqrt(sum(x=1,$slib_cnt,[$slib_weight][x]*([$slib_data][x]-slib_mean)^2));\ |
33 | print(slib_dev/sqrt($slib_tw)","slib_dev/sqrt($slib_tw-1)) |
||
34 | slib_out =!trim $slib_out |
||
20 | reyssat | 35 | !endif weight |
36 | !endif |