MACRO
aranova resp subs effs
# responses stacked in resp, group ids are in subs, *repeat* F-statistics
# will be put into effs
# copyright 2005 Robert L. Raymond University of St Thomas. Free use for
# teaching & noncommercial research - retain this notice in this and derivative works
MCONSTANT I repeat j len replic sxbars vxbars num sumssqd den
MCOLUMN resp subs effs inp rsubs rxbars rssqd
NOECHO
NOTE How many treatments are you comparing ?
SET inp;
   FILE "TERMINAL";
   NOBS 1.
LET I = inp(1)
NOTE How many F-statistics do you want ?
BRIEF 0
SET inp;
   FILE "TERMINAL";
   NOBS 1.
LET repeat = inp(1)
COUNT subs len
DIVIDE len I replic
DO j = 1:repeat
   SAMPLE len subs rsubs
   STATS resp;
      BY rsubs;
      MEAN rxbars;
      VARIANCE rssqd.
   STDEV rxbars sxbars
   MULTIPLY sxbars sxbars vxbars
   MULTIPLY vxbars replic num
   SUM rssqd sumssqd
   DIVIDE sumssqd I den
   LET effs(j) = num/den
ENDDO
BRIEF 2
ECHO
ENDMACRO