这的确是一般的时间序列
data(zeroyld)
zeroyld
short.run long.run
1 2.183 1.575
2 2.246 1.545
3 2.308 1.762
.....
.....
480 8.103 6.842
481 8.069 6.531
482 8.069 6.431
zhangtao兄,记得添加
library pgraph;
这个程序要执行很久,
最好修改一下参数.
初步看看执行的结果(当然准确度有差异).
gn=300; @ number of gridpoints for gamma @ ==>30
bn = 300; @ number of gridpoints for beta @ ==>30
boot=5000; @ number of bootstrap replications @ ==>500
作者除了提供Gauss code
也有提供 Rcode,Matlab code
|