SAS中文论坛

标题: 求助!关于proc/iml模块的FARMAFIT! [打印本页]

作者: shiyiming 时间: 2004-4-14 21:25
标题: 求助!关于proc/iml模块的FARMAFIT!
在work里有一数据集rv1,包含变量logrv,现要对该序列进行arfima模拟:

proc iml;
use rv1;
read all var{logrv};
call farmafit(d,ar,ma,sigma,logrv) p=1 q=1;
put d ar ma sigma;

错误提示:
ERROR: Cannot compute the function value.
ERROR: Read Access Violation In Task ( IML ]
Exception occurred at (67E691D3)
就是call 函数执行不了,不知何故,

请高人指点!

作者: shiyiming 时间: 2004-4-15 23:49
可能是你的sas有问题,俺这边可以运行 <img src="{SMILIES_PATH}/icon_twisted.gif" alt=":twisted:" title="Twisted Evil" />

作者: shiyiming 时间: 2004-4-16 15:23
标题: 不会吧!
不会吧, 如果真的是sas有问题那就难办了,

作者: shiyiming 时间: 2004-4-17 01:54
必须要更正一下,我用下面的做没问题.
proc iml;
d= 0.3;
phi = 0.5;
theta= -0.1;
call farmasim(yt, d, phi, theta);
call farmafit(d, ar, ma, sigma, yt) p=1 q=1;
print d ar ma sigma;
quit;

但是换成数据,用下面的程序做就出现和你一样的问题,读入的数据就是上面生成的 yt
data tep;
input yt;
cards;
-3.854829488
-0.781470481
0.76046816
2.646407414
0.747332333
...
-0.11705902
;run;
proc iml;
use tep;
read all var{yt} into yt;
call farmafit(d,ar,ma,sigma,yt) p=1 q=1;
print d ar ma sigma;
quit;
这个问题确实很怪,但绝对不会是程序的问题.

作者: shiyiming 时间: 2004-4-18 12:24
对呀,我也用你上面的方法试过,确实没有问题,

可是用read语句来读入数据后,再用call farmafit就不行了,

会是什么原因呢?

欢迎光临 SAS中文论坛 (http://mysas.net/forum/)