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