如下要求的一个程序编写

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

具体内容见附件

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：程序编写 程序编写

proc iml;
x=rannorm(j(10,1,0));
lagx={0}//x[1:nrow(x)-1];
T=((lagx<=0 & x>0)+(lagx>0 & x<=0))[+];
print T;

x0不知道怎么设置，设置为0

ziyenano 发表于 2012-9-6 10:59
proc iml;
x=rannorm(j(10,1,0));
lagx={0}//x[1:nrow(x)-1];

x=rannorm(j(10,1,0));不对吧，似乎应该为x=normal(j(10,1,0));吧，请检验一下，你给的这个语句我这里无法运行。

harlon1976 发表于 2012-9-6 14:43
x=rannorm(j(10,1,0));不对吧，似乎应该为x=normal(j(10,1,0));吧，请检验一下，你给的这个语句我这里无 ...

normal,rannorm都是标准正态分布的随机数；
语句没有问题，可以运行啊

ziyenano 发表于 2012-9-6 14:52
normal,rannorm都是标准正态分布的随机数；
语句没有问题，可以运行啊

在IML中，我这里运行不起来，在DATA步中rannorm是可以的，我再问你一个问题，在IML中不用循环语句如何产生一个初值为0的随机游走模型,且样本为T个，即x(t)=x(t-1)+e(t),x(0)=0,e(t)可以是标准正态分布，也可以是其它的分布如卡方分布和t分布等。谢谢！

应该是可以的吧，这边版本是9.2，你报什么错，截个图；
初值为0随机游走模型，A=normal(j(T,1,0))，实现A元素的累加，最后加上初值为0，得到序列；
不过iml好像没有累加函数，或者生成一个上三角全为1的方阵，A*这个方阵，相当于累加功能

ziyenano 发表于 2012-9-6 16:17
应该是可以的吧，这边版本是9.2，你报什么错，截个图；
初值为0随机游走模型，A=normal(j(T,1,0))，实现A元 ...

这个思路是可行的，回头我试试，我用的是9.1.3版本的，运行后有日志窗口显示：
x=rannorm(j(10,1,0));
ERROR: Invocation of unresolved module RANNORM.

可否告诉我你的QQ，我以后多向您请教SAS程序。