Subversion Repositories wimsdev

Rev

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