Subversion Repositories wimsdev

Rev

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