Subversion Repositories wimsdev

Rev

Rev 7692 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 7692 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=3\
 5 
slib_parms=3\
6 
 1,m\
 6 
 1,m\
7 
  ,n\
 7 
  ,n\
Line 20... Line 20...
20 
slib_q=!declosing $slib_q
 20 
slib_q=!declosing $slib_q
21 
slib_t=!itemcnt $slib_q
 21 
slib_t=!itemcnt $slib_q
22 
 
 22 
 
23 
slib_s=!sum x for x in $slib_q
 23 
slib_s=!sum x for x in $slib_q
24 
!if $slib_s>1
 24 
!if $slib_s>1
25 
!exit
 25 
  !exit
26 
!endif
 26 
!endif
27 
!if $slib_s>0 and $slib_s<1
 27 
!if $slib_s>0 and $slib_s<1
28 
slib_q=$slib_q, $[1-$slib_s]
 28 
  slib_q=$slib_q, $[1-$slib_s]
29 
!advance slib_t
 29 
  !advance slib_t
30 
!endif
 30 
!endif
31 
 
 31 
 
32 
slib_mult=
 32 
slib_mult=
33 
!for slib_r=1 to $slib_M
 33 
!for slib_r=1 to $slib_M
34 
slib_s=1
 34 
  slib_s=1
35 
slib_N=$slib_nn
 35 
  slib_N=$slib_nn
36 
 
 36 
 
37 
slib_V=
 37 
  slib_V=
38 
!for slib_i=1 to  $[$slib_t-1]
 38 
  !for slib_i=1 to  $[$slib_t-1]
39 
slib_j=!item $slib_i of $slib_q
 39 
    slib_j=!item $slib_i of $slib_q
40 
!readproc slib/stat/binomial  1,$slib_N, $[$slib_j/$slib_s]
 40 
    !readproc slib/stat/binomial  1,$slib_N, $[$slib_j/$slib_s]
41 
 
 41 
 
42 
slib_V=!append item $slib_out to $slib_V
 42 
    slib_V=!append item $slib_out to $slib_V
43 
!distribute item $[$slib_s-$slib_j], $[$slib_N-$slib_out] into slib_s,slib_N
 43 
    !distribute item $[$slib_s-$slib_j], $[$slib_N-$slib_out] into slib_s,slib_N
44 
 
 44 
 
45 
!next slib_i
 45 
  !next slib_i
46 
slib_V=!append item $slib_N to $slib_V
 46 
  slib_V=!append item $slib_N to $slib_V
47 
 
 47 
 
48 
slib_mult=!append line $slib_V to $slib_mult
 48 
  slib_mult=!append line $slib_V to $slib_mult
49 
!next slib_r
 49 
!next slib_r
50 
 
 50 
 
51 
slib_out=!nonempty line $slib_mult
 51 
slib_out=!nonempty line $slib_mult
52 
slib_out=!trim $slib_out
 52 
slib_out=!trim $slib_out