/* sandm97.sas                     */
 /*  analysis of sandm77.r output   */
 /*                                */
 /* three methods, two param ests, */
 /*  conv codes, sigma-sq-a        */
  options ls=80 ;
  data a ;
    infile 'estsP47.dat' ;          * read data ;
    keep sqere sqera method tensiga rep ;
    input a1 a2 m1 m2 r1 r2 cm cr siga ;
    rep = mod(_n_-1,100) ;
    tensiga = put(10*siga,3.) ;
 /* make 3 obs from each -- bs & i  */
    method = 'anova' ;
    sqere = (a1-1)**2 ;
    sqera = (a2-siga)**2 ;
    output ;
    method= 'ml' ;
    sqere = (m1-1)**2 ;
    sqera = (m2-siga)**2 ;
    output ;
    method = 'reml' ;
    sqere = (r1-1)**2 ;
    sqera = (r2-siga)**2 ;
    output ;
 run ;
 /* analyze to compare mse         */
 /*  means & stdev of estimators   */
  proc glm data=a ;
    class method tensiga rep ;
    model sqera = method tensiga tensiga*rep ;
    random tensiga*rep ;
    lsmeans method tensiga ;
    title 'Analysis of MSE' ;
 run ;