作者tew (咖啡王子)
看板Statistics
标题Re: [程式] sas矩阵相乘
时间Thu Sep 29 15:39:51 2016
※ 引述《HeliosPlus (离风)》之铭言：
: TITLE请注名软体类别 如[程式] stata
: 发文前请把灰色的注解删除（可使用ctrl+y删除一整行）
: [软体程式类别]:
: SAS
: [程式问题]:
: 资料处理
: [软体熟悉度]:
: 新手
: [问题叙述]:
: 小弟我已经汇入了两分资料，A,B
: 想要把AB以矩阵方式相乘，请问该怎麽做呢？
: 网路上的教学都是手动建立矩阵，但我的资料量很大，一笔一笔打会想死…
: 谢谢大家帮忙
data first;
do i=1 to 100;
a=rannor(1);
b=rannor(2);
c=rannor(3);
output;
end;
drop i;
run;
data second;
do j=1 to 100;
x1=rannor(2+j);
x2=rannor(x1);
x3=rannor(3);
output;
end;
drop j;
run;
/*
建立两个档案
*/
proc iml;
use first var _all_;
/*读入第一个档案的所有变量数值*/
read all into a;
/*将之创建为矩阵A*/
use second var _all_;
/*读入第二个档案的所有变量数值*/
read all into b;
/*将之创建为矩阵b*/
x=a`*b;
/*矩阵X等於两矩阵相乘*/
print x;
quit;
好像也没有所谓的手动制作矩阵的说法
除非你是想求行列式的数值
--
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 59.72.30.154
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1475134794.A.681.html