求助：关于iml解矩阵

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我是SAS/IML的新手。想解一个特定的矩阵。
要是单独解一个公司（比如公司A或者B），求出PV值的话，没有什么问题。可是如果要是想一起解两个公司的PV值的话怎么解不出来啊？大家能帮帮我吗，非常感谢。。。

这是我的编程 (由于字数限制，问题在3楼）：


1)求公司A的PV值 (没有问题）：

proc iml;

use lastAccumImpulse;read all where(ric='A' )into x;
x1=x[1,1];
x2=x[1,2];
x3=x[1,3];
x4=x[1,4];

use CovInnovation;
read all where(ric='A' )into y ;
y1=y[1:4,1];
y2=y[1:4,2];
y3=y[1:4,3];
y4=y[1:4,4];

z=T(x);

PV=x*y*z;

print x y z PV;
quit;

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：解矩阵 Innovation impulse Where print where

1# 大书页
2）求公司B的PV值 （没有问题）：

proc iml;

use lastAccumImpulse;

read all where(ric='B' )into x;

x1=x[1,1];

x2=x[1,2];

x3=x[1,3];

x4=x[1,4];

use CovInnovation;

read all where(ric='B' )into y ;

y1=y[1:4,1];

y2=y[1:4,2];

y3=y[1:4,3];

y4=y[1:4,4];

z=T(x);
PV=x*y*z;

print x y z PV;
quit;

2# 大书页
1）
3)求公司A 和B的PV值：

大家帮我看看怎么改才能一次性求出来他们的PV值（下面的是不正确的，大家帮我看看怎么改）：

proc
iml;

use lastAccumImpulse;

read all into x;

x1=x[1,1];

x2=x[1,2];

x3=x[1,3];

x4=x[1,4];

use CovInnovation;

read all into y ;

y1=y[1:4,1];

y2=y[1:4,2];

y3=y[1:4,3];

y4=y[1:4,4];

z=T(x);

PV=x*y*z;

print x y z PV;

quit;

各位大大，有熟悉IML的吗？帮帮忙。。。。。。
非常感谢。。。。。。。。。

我不熟悉sas的iml，但是给你本参考的，翻翻里面有没有吧~~
s Guide.pdf (4.55 MB)

2009-8-18 13:17:05 上传

哦对了，下这个不花钱的

非常感谢，哈哈，我去看看