/* sandm8.sas                     */
 /*  analysis of sandm7.r output   */
 /*                                */
 /* three methods, two param ests, */
 /*  conv codes, sigma-sq-a        */
  options ls=80 ;
  data a ;
    infile 'estsP4.dat' ;          * read data ;
    keep dseam dsear dsemr a2 m2 r2 siga ;
    input a1 a2 m1 m2 r1 r2 cm cr siga ;
 /*                                */
 /* these for comparing mse's      */
 /*                                */
    dseam = (a2-siga)**2 - (m2-siga)**2 ; 
    dsear = (a2-siga)**2 - (r2-siga)**2 ; 
    dsemr = (m2-siga)**2 - (r2-siga)**2 ;
 /* add labels for easier reading  */
 /*                                */ 
   label a1='Anova e' a2='Anova a'
         m1='MLE e'   m2='MLE a'
         r1='RML e'   r2='RML a'
         cm='codmle'  cr='codrml' ;
 run ;
 /* analyze to compare mse         */
 /* means & variance of estimators */
  proc means data=a mean stderr stddev ;
    class siga ;
    var a2 m2 r2 ;
    title 'Analysis of bias and variance' ;
 run ;
 /*  t-statistics for pairwise     */
 /*   comparison of mse's          */
 /*                                */
  proc means data=a mean stderr t ;
    class siga ;
    var dseam dsear dsemr ;
    title 'Analysis of Pairwise differences' ;
 run ;