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