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