yydy322a.sas

高老师的经典多元统计分析教材源代码非常经典使用

/*  yydy322a.sas  */

 data  d322;
   input  type x1-x4;
   cards;
 1      65     35     25     60
 1      75     50     20     55
 1      60     45     35     65
 1      75     40     40     70
 1      70     30     30     50
 1      55     40     35     65
 1      60     45     30     60
 1      65     40     25     60
 1      60     50     30     70
 1      55     55     35     75
 2      55     55     40     65
 2      50     60     45     70
 2      45     45     35     75
 2      50     50     50     70
 2      55     50     30     75
 2      60     40     45     60
 2      65     55     45     75
 2      50     60     35     80
 2      40     45     30     65
 2      45     50     45     70
 ;

  proc iml;
   n=10;m=10; p=4;
   use d322(obs=10);
   xx={x1 x2 x3 x4};
   read all  var xx into x; print x;
   ln={[10] 1} ;
   x0=(ln*x)/n; print x0;
   mx=i(n)-j(n,n,1)/n;
   a1=x`*mx*x; print a1;

   use d322(firstobs=11);
   read all  var xx into y; print y;
   lm={[10] 1} ;
   y0=(lm*y)/m; print y0;
   my=i(m)-j(m,m,1)/m;
   a2=y`*my*y; print a2;

   a=a1+a2; xy=x0-y0;
   ai=inv(a); print a ai;
   dd=xy*ai*xy`; d2=(m+n-2)*dd;
   t2=n*m*d2/(n+m) ;
   f=(n+m-1-p)*t2/((n+m-2)*p);
   print  d2 t2 f;
   pp=1-probf(f,p,m+n-p-1);
   print pp;
 run;