Anonymous user
Permutation test: Difference between revisions
m
→{{header|REXX}}: added more whitespace, changed some comments.
m (→{{header|REXX}}: added more whitespace, changed some comments.) |
|||
Line 940:
This REXX program is modeled after the C version, with some generalizations and optimization.
<lang rexx>/*REXX program does a permutation test on N + M subjects (volunteers):*/
n=9
data=85 88 75 66 25 29 83 39 97 68 41 10 49 16 65 32 92 28 98
w=words(data); m=w-n
say 'volunteer population given treatment:' right(n,length(w))
say ' control population given a placebo:' right(m,length(w))
say
say 'treatment population efficacy % (percentages):' subword(data,1,n)
say ' control population placebo % (percentages):' subword(data,n+1)
say
do v= 0 for w ; #.v=word(data,v+1) ; end
treat=0; do i= 0 to n-1 ; treat=treat+#.i ; end
total=1; do j=19 to m+1 by -1 ; total=total*j ; end
do k= 9 to 1 by -1
gt=
le=total-gt
say "<= " format(100*le/total,,3)'%'
say " > " format(100*gt/total,,3)'%' gt
exit /*stick a fork in it, we're done.*/
/*──────────────────────────────────PICKER subroutine───────────────────*/
if rest==0 then return eff>treat
if it>rest then q=
else q=0
itP=it-1
return
'''output''' using the default input:
<pre>
volunteer population given treatment: 9
|