求助!关于proc/iml模块的FARMAFIT!

shiyiming

在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

可能是你的sas有问题,俺这边可以运行 <!-- s:twisted: --><img src="{SMILIES_PATH}/icon_twisted.gif" alt=":twisted:" title="Twisted Evil" /><!-- s:twisted: -->

shiyiming

不会吧, 如果真的是sas有问题那就难办了,

shiyiming

必须要更正一下,我用下面的做没问题.
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

对呀,我也用你上面的方法试过,确实没有问题,

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

会是什么原因呢?