方法如下:
1、比如说是2阶方阵,现在对其中的元素进行一些运算并把结果重新赋值到另一个矩阵中,
proc iml;
beta={2 3,3 4};
q12=beta[1,1]/2;
q22=beta[2,2]**2/3;
q={1 q12,q12 q22};
print beta q12 q22 q;
2、方法:
proc iml;
beta={2 3,3 4};
q=j(2,2,.);
q[1,2]=beta[1,1]/2;
q[2,2]=beta[2,2]**2/3;
q[1,1]=1;
q[2,1]=beta[1,1]/2;
print beta q;
%if you have a matrix A
idx=find(A==1.5); % find all 1.5
A(idx)=1; % set 1 to these indexes
%about NaN
idx=find(isnan(A)); % find all NaN value
A(idx)=1; % set 1 to these indexes