1.帮你做了一个简单的例子,供参考
内含archTest, unit root test, modified R/S test,
figarch, fiegarch, and model checking.
请把你的数据和ibm.txt作一比较
#####
module(finmetrics)
ibm=matrix(scan(file='c:/ibm.txt'))
par(mfrow=c(2,1))
acf(ibm,24)
acf(ibm^2,24)
archTest(ibm,lag=24)
# UNIT ROOT ADF TEST
adft.out = unitroot(ibm, trend="nc", statistic="t", lags=1)
adft.out
# Modified R/S Long Memory Test
rosTest(abs(ibm))
# use percentage returns
ibm=ibm*100
#figarch
ibm.mod1=fgarch(ibm~1,~figarch(1,1))
summary(ibm.mod1)
# figarch model checking
autocorTest(residuals(ibm.mod1,standarize=T),lag=20)
archTest(residuals(ibm.mod1,standarize=T),lag=20)
#fiegarch
ibm.mod2=fgarch(ibm~1,~fiegarch(1,1),leverage=T)
summary(ibm.mod2)
# fiegarch model checking
autocorTest(residuals(ibm.mod2,standarize=T),lag=20)
archTest(residuals(ibm.mod2,standarize=T),lag=20)
plot(ibm.mod1)
[此贴子已经被作者于2007-6-12 7:36:59编辑过]