Rev 7349 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
| Rev 7349 | 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 | ,[datax;datay] or [datax],[datay] |
6 | ,[datax;datay] or [datax],[datay] |
| 7 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
7 | slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou |
| Line 14... | Line 14... | ||
| 14 | :proc |
14 | :proc |
| 15 | slib_out= |
15 | slib_out= |
| 16 | !readproc slib/stat/dataproc $wims_read_parm |
16 | !readproc slib/stat/dataproc $wims_read_parm |
| 17 | 17 | ||
| 18 | !if $slib_weight=$empty or $slib_data=$empty |
18 | !if $slib_weight=$empty or $slib_data=$empty |
| 19 | !exit |
19 | !exit |
| 20 | !endif |
20 | !endif |
| 21 | slib_x=$slib_data |
21 | slib_x=$slib_data |
| 22 | slib_y=$slib_weight |
22 | slib_y=$slib_weight |
| 23 | slib_cnt=!itemcnt $slib_x |
23 | slib_cnt=!itemcnt $slib_x |
| 24 | slib_cnty=!itemcnt $slib_y |
24 | slib_cnty=!itemcnt $slib_y |
| 25 | !if $slib_cnt<>$slib_cnty |
25 | !if $slib_cnt<>$slib_cnty |
| 26 | !exit |
26 | !exit |
| 27 | !endif |
27 | !endif |
| 28 | 28 | ||
| 29 | slib_meanx=!sum x for x in $slib_x |
29 | slib_meanx=!sum x for x in $slib_x |
| 30 | slib_meanx=$[($slib_meanx)/$slib_cnt] |
30 | slib_meanx=$[($slib_meanx)/$slib_cnt] |
| 31 | slib_meany=!sum y for y in $slib_y |
31 | slib_meany=!sum y for y in $slib_y |
| 32 | slib_meany=$[($slib_meany)/$slib_cnt] |
32 | slib_meany=$[($slib_meany)/$slib_cnt] |
| 33 | 33 | ||
| 34 | slib_varx=!sum (x-($slib_meanx))^2 for x in $slib_x |
34 | slib_varx=!sum (x-($slib_meanx))^2 for x in $slib_x |
| 35 | slib_varx=$[$slib_varx] |
35 | slib_varx=$[$slib_varx] |
| 36 | 36 | ||
| 37 | slib_vary=!sum (y-($slib_meany))^2 for y in $slib_y |
37 | slib_vary=!sum (y-($slib_meany))^2 for y in $slib_y |
| 38 | slib_vary=$[$slib_vary] |
38 | slib_vary=$[$slib_vary] |
| 39 | 39 | ||
| 40 | slib_covar=!exec pari print(sum(slib_t=1,$slib_cnt,([$slib_x][slib_t]-($slib_meanx))*([$slib_y][slib_t]-($slib_meany)))); |
40 | slib_covar=!exec pari print(sum(slib_t=1,$slib_cnt,([$slib_x][slib_t]-($slib_meanx))*([$slib_y][slib_t]-($slib_meany)))); |
| 41 | 41 | ||
| 42 | slib_out=$[($slib_varx)/$slib_cnt],$[($slib_covar)/$slib_cnt];$[($slib_covar)/$slib_cnt],$[($slib_vary)/$slib_cnt] |
42 | slib_out=$[($slib_varx)/$slib_cnt],$[($slib_covar)/$slib_cnt];$[($slib_covar)/$slib_cnt],$[($slib_vary)/$slib_cnt] |
| 43 | slib_out=!append line $[($slib_varx)/($slib_cnt-1)],$[($slib_covar)/($slib_cnt-1)];$[($slib_covar)/($slib_cnt-1)],$[($slib_vary)/($slib_cnt-1)] to $slib_out |
43 | slib_out=!append line $[($slib_varx)/($slib_cnt-1)],$[($slib_covar)/($slib_cnt-1)];$[($slib_covar)/($slib_cnt-1)],$[($slib_vary)/($slib_cnt-1)] to $slib_out |
| 44 | - | ||
| 45 | - | ||