Subversion Repositories wimsdev

Rev

Rev 4161 | Go to most recent revision | Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
20 reyssat 1
!if $wims_read_parm!=slib_header
2
!goto proc
3
!endif
4
slib_title=Quantile of Binomial law
5
slib_parms=3\
6
,[data]\
7
5,n\
8
0.5,p
9
 
10
slib_author=Bernadette PERRIN-RIOU
11
slib_out=quantilei of binomial law of parameters n and p ((0<p<1)) at data
12
slib_comment=data can be a list of reals
13
slib_example=[0.2,0.4,0.6],5,0.5
14
!exit
15
 
16
:proc
17
 
18
!distribute  item $wims_read_parm into slib_data,slib_n,slib_p
19
!default slib_n= 5
20
!default slib_p=0.5
21
!if $slib_p>=1 or $slib_p<=0 
22
	slib_out=0
23
	!exit
24
!endif
25
slib_out=
26
slib_data=!declosing $slib_data
27
slib_p1=$[1-$slib_p]
28
slib_r=$[$slib_p/$slib_p1]
29
slib_a=$[$slib_p1^($slib_n)]
30
!!slib_bino=!values binomial($slib_n,x)*$slib_p^x*$slib_p1^($slib_n-x) for x =0 to $slib_n
31
slib_cnt= !itemcnt $slib_data
32
slib_pp=$[$slib_p/$slib_p1]
33
!for slib_i =1 to $slib_cnt
34
	slib_di=!item $slib_i of $slib_data
35
	slib_P=0
36
	slib_a=$[$slib_p1^($slib_n)]
37
	!for slib_k=0 to $slib_n
38
		!!slib_a=!item $[$slib_k+1] of $slib_bino
39
		slib_P=$[$slib_P+$slib_a]
40
		!if $slib_P>= $slib_di
41
			slib_out=!append item  $slib_k to $slib_out
42
			!break
43
		!endif
44
		slib_a=$[$slib_a*$slib_pp*($slib_n-$slib_k)/($slib_k+1)]
45
	!next slib_k
46
!next slib_i
47
slib_out=!trim $slib_out