Question

我有一个行矩阵（向量）A和另一个方阵B.如何使用proc iml或其他方法将矩阵B的每一行与SAS中的行矩阵A相乘？

让我们说

a = {1 2 3}
b =
{2 3 4
1 5 3
5 9 10}

我的输出c将是：

{2 6 12
1 10 9
5 18 30}

谢谢！

Answer 1

在IML中使用逐元素乘法运算符#：

proc iml;
a = {1 2 3};
b = {2 3 4, 
     1 5 3,
     5 9 10};

c = a#b;
print c;
quit;

Answer 2

当然有一个非IML解决方案，或二十，虽然IML，因为Dom注释可能是最简单的。这是两个。

首先，将它们放到一个数据集上，其中a数据集位于每一行（带有一些其他变量名称） - 见下文。然后，要么只是进行数学运算（使用数组），要么使用PROC MEANS或类似数据集作为权重使用a数据集。

data a;
  input w_x w_y w_z;
  datalines;
1 2 3
;;;;
run;

data b;
  input x y z;
  id=_n_;
  datalines;
2 3 4
1 5 3
5 9 10
;;;;
run;

data b_a;
  if _n_=1 then set a;
  set b;
  *you could just multiply things here if you wanted;
run;


proc means data=b_a;
class id;
types id;
var x/weight=w_x;
var y/weight=w_y;
var z/weight=w_z;
output out=want sum=;
run;

如何在sas中将每一行与另一个矩阵元素的每一行相乘？

2 个答案: