同求

DATA matrix;
INPUT a b c;
CARDS;
1 0.4 0.3
0.8 2 0.2
0.7 0.1 0.6
;
RUN;

proc iml;
use matrix;
read all into g;
close matrix;
P=g-i(nrow(g));
t_P=P`;
Q=P`*P;
eig_Q=eigval(Q);
print P Q t_P eig_Q;
quit;

谢谢！不过如果不能用IML模块的话该如何呢？我们的SAS版本里没有IML的授权，

ans1:
data p(drop=i);
set matrix;
array p(*) a b c;
do i=1 to dim(p);
   if i=_n_ then p(i)=p(i)-1;
end;
run;

ans2:
proc transpose data=p out=t_p;
run;

ans3:
data t1(rename=(a=a1 b=b1 c=c1));
set p p p;
run;

data t2;
if mod(_n_-1,3)=0 then set p;
set t1;
array p1(*) a b c;
array p2(*) a1 b1 c1;
p=0;
do i=1 to dim(p1);
   p=p+p1(i)*p2(i);
end;
run;


data Q;
array q(3);
do i=1 to 3;
   set t2(keep=p);
       q(i)=p;
end;
drop p i;
run;

用data步写矩阵太难受了，特征值不晓得怎么用data步写，
proc fcmp支持有限的矩阵运算，不过不是很好用，
实在不行用R吧，算出来结果再传回SAS

这是另一种方式来计算矩阵，但不包括特征值和特征向量.
如果数据不多， 可以直接在网上计算特征值和特征向量.