Subversion Repositories wimsdev

Rev

Rev 7692 | 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=2\
 5 
slib_parms=2\
6 
  2,m\
 6 
  2,m\
7 
  2*m,N\
 7 
  2*m,N\
Line 20... Line 20...
20 
!default slib_mm=2
 20 
!default slib_mm=2
21 
!default slib_NN=$[2*$slib_mm]
 21 
!default slib_NN=$[2*$slib_mm]
22 
!default slib_K= R
 22 
!default slib_K= R
23 
 
 23 
 
24 
!if $slib_mm = 1
 24 
!if $slib_mm = 1
25 
 slib_out =1
 25 
  slib_out =1
26 
!else
 26 
!else
27 
 slib_cc = !values 1/$slib_mm for x = 1 to $[$slib_mm-1]
 27 
  slib_cc = !values 1/$slib_mm for x = 1 to $[$slib_mm-1]
28 
 
 28 
 
29 
 !readproc slib/stat/multinomial 1, $[$slib_NN-$slib_mm], $slib_cc
 29 
  !readproc slib/stat/multinomial 1, $[$slib_NN-$slib_mm], $slib_cc
30 
 slib_vv= $slib_out
 30 
  slib_vv= $slib_out
31 
 
 31 
 
32 
 !if $slib_K issametext Q
 32 
  !if $slib_K issametext Q
33 
  slib_out=!exec pari slib_a=Vec([$slib_vv]);vector($slib_mm, x,(slib_a[x]+1)/$slib_NN)
 33 
    slib_out=!exec pari slib_a=Vec([$slib_vv]);vector($slib_mm, x,(slib_a[x]+1)/$slib_NN)
34 
 !else
 34 
  !else
35 
  slib_out=!values (x+1)/$slib_NN for x in $slib_vv
 35 
    slib_out=!values (x+1)/$slib_NN for x in $slib_vv
36 
 !endif
 36 
  !endif
37 
!endif
 37 
!endif
38 
slib_out =!trim $slib_out
 38 
slib_out =!trim $slib_out