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 |