我的使用也遇到问题了,向高手求助。我用的arfima软件包和ox edit。因为我要求d值,所以我把预测5的程序直接搬过来了,程序如下:#include <oxstd.h> #import <packages/arfima/arfima>
main()
{
decl dy = loadmat("orginal.in7");
decl dly = diff0(log(dy), 1)[1:][];
print(dly);
Arfima::EstimateGPH(dly, sizerc(dly) / 2, TRUE);
Arfima::EstimateGSP(dly, sizerc(dly) / 2, TRUE);
}
我的数据是指数形式的,我发现通过decl dy = loadmat("orginal.in7"); decl dly = diff0(log(dy), 1)[1:][];已经把数据处理成先Log再差分的形式,这点通过我加了个 print(dly);可以看出来。
可是,算出来的gph的d结果跟我用其他方法算出来的不一样。
结果:---- Log Periodogram Regression ----
d parameter 0.812188 (0.120608) [0.0000]
No of observations: 95; no of periodogram points: 47
---- Gaussian semiparametric estimate ----
d parameter 0.600154 (0.0729325) [0.0000]
No of observations: 95; no of periodogram points: 47
我不明白0.812188 和0.600154就是结果吗?结果应该在0.5以下才和我用其他方法做的一致啊?
请高人指点。