Subversion Repositories wimsdev

Rev

Rev 4330 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4330 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 16... Line 16...
16 
 
 16 
 
17 
!readproc slib/stat/dataproc $wims_read_parm
 17 
!readproc slib/stat/dataproc $wims_read_parm
18 
slib_x=$slib_x
 18 
slib_x=$slib_x
19 
slib_y=$slib_y
 19 
slib_y=$slib_y
20 
!if $slib_y=$empty or $slib_x=$empty
 20 
!if $slib_y=$empty or $slib_x=$empty
21 
 !exit
 21 
  !exit
22 
!endif
 22 
!endif
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/$slib_cnt]
 35 
!!slib_varx=$[$slib_varx/$slib_cnt]
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/$slib_cnt]
 38 
!!slib_vary=$[$slib_vary/$slib_cnt]
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=$[sqrt(($slib_varx)/($slib_vary))],$[$slib_covar/sqrt($slib_varx*($slib_vary))];$[$slib_covar/sqrt(($slib_varx)*($slib_vary))],$[sqrt(($slib_vary)/($slib_varx))]
 42 
slib_out=$[sqrt(($slib_varx)/($slib_vary))],$[$slib_covar/sqrt($slib_varx*($slib_vary))];$[$slib_covar/sqrt(($slib_varx)*($slib_vary))],$[sqrt(($slib_vary)/($slib_varx))]
43 
 
 - 
 
44 
 
 -