为什么beta输出来的不是数组呢??提示错误 (execution) Matrix has not been set to a value.
proc iml ;
A={0.483 13.2682 0 4.3646 5.1070,
0.4035 13.4909 39.0131 3.6151 5.5005,
0.8979 25.7776 9.0513 4.8920 7.5342,
0.5927 16.0245 13.2935 4.4529 6.5913
};
r=nrow(A);
c=ncol(A);
B=A;
do i=1 to c;
do p=1 to r;
if i=3
then B[p,i]=(A[<>,i]-A[p,i])/(A[<>,i]-A[><,i]);
else
B[p,i]=(A[p,i]-A[><,i])/(A[<>,i]-A[><,i]);
end;
end;
do i=1 to c;
deta[i]=std(B[1,i],B[2,i],B[3,i],B[4,i]);
end;
print deta;
quit;