1

我有两列（c1，c2）的 ax 矩阵。我想修复第一列（c1），将 10 列添加到 X 矩阵中，每列都有值 C2+m、C2+m...C2+m，m 是一个随机整数。最后矩阵将是：

C1, C2+m, C2+m, C2+m...C2+m;

代码：

proc iml;

use nonpar;

read all var{treat response} into x;

do i=1 to 10;

call randseed(123);

call randgen(u, "Uniform");

Max = 300; Min = 68;

m = min + floor( (1+Max-Min)*u );

x = x[,1]||x[,2]+m;

end;

quit;

有人可以帮我解决这个问题。。谢谢

4

1 回答 1

1

有几件事可以引导您朝着正确的方向前进。

首先，预先创建完整的目标矩阵；不要经常连接。因此，一旦您将数据集读入x，请创建另一个x_new与行数相同x但有 11 列的数据集。 j会为你做这件事。

其次，您可以一次生成所有随机数，但您必须首先定义要填充的矩阵的大小，再次使用j. 这是假设您希望 10 列中的每一列和每一行都有一个新的随机整数；如果您只想要每一行或总共只有一个“m”，您需要以不同的方式执行此操作，但您需要澄清这一点。如果您只想要一行 10 m，那么您可以先执行此操作（生成u具有 10 列 1 行的 a），然后使用矩阵乘法将其扩展为 x 的全部行数。

这是一个使用 SASHELP.CLASS 的简化示例，展示了这两个概念的作用。

proc iml;
  use sashelp.class;
  read all var {age weight} into x;
  x_new = j(nrow(x),11);  *making the new matrix with lots of columns;
  x_new[,1] = x[,1];      *copy in column 1;
  call randseed(123); 
  u = j(nrow(x),10);      *make the to be filled random matrix;
  call randgen(u,'Uniform',68,300); *the min/max parameters can go here;
  u = floor(u+0.5);  *as Rick noted in comments, needed to get the max value properly;
  x_new[,2:11] = u[,1:10] + x[,2]; *populate x_new here;
  print x_new;
quit;
于 2014-09-24T16:30:55.037 回答