epoh老师，帮我看看吧~实在搞不定

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

这个就是我的序列数据，怎么都搞不成时间序列格式，拜托你看看吧

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏4 回帖

关键词：epoh SOSO 序列数据 时间序列 sos 老师

module(finmetrics)
data=importData(file="pitaya.xls", type="EXCEL")
pitaya=cbind(data$ndfr,data$ahpr)
pitaya.bekk = mgarch(pitaya~-1, ~bekk(1,1))
summary(pitaya.bekk)
Call:
mgarch(formula.mean = pitaya ~ -1, formula.var =  ~ bekk(1, 1))

Mean Equation: structure(.Data = pitaya ~ -1
, class = "formula"
)

Conditional Variance Equation: structure(.Data =  ~ bekk(1, 1)
, class = "formula"
)

Conditional Distribution:  gaussian

--------------------------------------------------------------

Estimated Coefficients:
--------------------------------------------------------------
                   Value Std.Error   t value   Pr(>|t|)
       A(1, 1)  0.085875  0.008769    9.7933 0.000e+000
       A(2, 1)  0.028472  0.084970    0.3351 7.376e-001
       A(2, 2)  0.169469  0.098120    1.7272 8.450e-002
ARCH(1; 1, 1)  0.441704  0.029476   14.9854 0.000e+000
ARCH(1; 2, 1) -0.684249  0.248628   -2.7521 6.048e-003
ARCH(1; 1, 2)  0.009955  0.004621    2.1545 3.148e-002
ARCH(1; 2, 2)  0.234791  0.029056    8.0806 2.220e-015
GARCH(1; 1, 1)  0.830065  0.022975   36.1292 0.000e+000
GARCH(1; 2, 1)  0.325959  0.160278    2.0337 4.229e-002
GARCH(1; 1, 2) -0.002823  0.002767   -1.0200 3.080e-001
GARCH(1; 2, 2)  0.962251  0.008904  108.0645 0.000e+000

--------------------------------------------------------------

AIC(11) = 3260.095
BIC(11) = 3312.305

Normality Test:
--------------------------------------------------------------
   Jarque-Bera    P-value Shapiro-Wilk P-value
v1      182.89 0.000e+000       0.9813 0.02085
v2       28.76 5.682e-007       0.9879 0.74634

Ljung-Box test for standardized residuals:
--------------------------------------------------------------
   Statistic P-value Chi^2-d.f.
v1     16.01  0.1906         12
v2     23.44  0.0242         12

Ljung-Box test for squared standardized residuals:
--------------------------------------------------------------
   Statistic P-value Chi^2-d.f.
v1     6.632 0.88094         12
v2    23.236 0.02579         12

Lagrange multiplier test:
--------------------------------------------------------------
    Lag 1    Lag 2   Lag 3   Lag 4  Lag 5   Lag 6   Lag 7  Lag 8   Lag 9  Lag 10 Lag 11   Lag 12        C
v1 -1.035 -0.01092 -0.9605  0.5189 0.0635 -0.8196  0.0453 0.9842 -0.5974  0.7353 1.1013  0.02773 -0.06196
v2  1.274 -0.22540  3.3285 -0.1522 1.3225 -0.4175 -2.1950 0.2379 -1.6314 -0.6673 0.2407 -0.14026  0.85735

     TR^2 P-value F-stat P-value
v1  6.203 0.90551 0.5681 0.94316
v2 22.598 0.03134 2.1112 0.08006

epoh 发表于 2013-2-20 13:28
module(finmetrics)
data=importData(file="pitaya.xls", type="EXCEL")
pitaya=cbind(data$ndfr,data$a ...

谢谢您！！！！！！！！！！！！！~~~~~~~~~~~~~

epoh 发表于 2013-2-20 13:28
module(finmetrics)
data=importData(file="pitaya.xls", type="EXCEL")
pitaya=cbind(data$ndfr,data$a ...

再咨询您一下，我这个生成的BEKK没有C矩阵，是不是做wald检验的时候只要定义4*11的矩阵就行了[1,5]=[2,6]=[3,9]=[4,11]=1检验双波动溢出，然后令[1,5]=[2,9]=1检验AtoB的，令[1,6]=[2,10]=1检验BtoA的？

火龙果汽水 发表于 2013-2-20 14:23
再咨询您一下，我这个生成的BEKK没有C矩阵，是不是做wald检验的时候只要定义4*11的矩阵就行了[1,5]=[2,6] ...

应该是[1,5]=[2,6]=[3,9]=[4,10]=1  ??

epoh 发表于 2013-2-20 14:58
应该是[1,5]=[2,6]=[3,9]=[4,10]=1  ??

数据文件给您传上了，就是那个叫做数据的

火龙果汽水 发表于 2013-2-23 14:31
数据文件给您传上了，就是那个叫做数据的

module(finmetrics)
data=importData(file="pitaya1.xls", type="EXCEL")
pitaya1=cbind(data$v1,data$v2)
pitaya1.bekk = mgarch(pitaya1~-1, ~bekk(1,1))
summary(pitaya1.bekk)

Call:
mgarch(formula.mean = pitaya1 ~ -1, formula.var =  ~ bekk(1, 1))

Mean Equation: structure(.Data = pitaya1 ~ -1
, class = "formula"
)

Conditional Variance Equation: structure(.Data =  ~ bekk(1, 1)
, class = "formula"
)

Conditional Distribution:  gaussian

--------------------------------------------------------------

Estimated Coefficients:
--------------------------------------------------------------
                   Value Std.Error   t value   Pr(>|t|)
       A(1, 1)  0.001321  0.005064   0.26095 7.942e-001
       A(2, 1) -0.001187  0.011400  -0.10416 9.171e-001
       A(2, 2)  0.003367  0.005387   0.62506 5.322e-001
ARCH(1; 1, 1)  0.293512  0.082547   3.55569 4.125e-004
ARCH(1; 2, 1) -0.001479  0.083964  -0.01761 9.860e-001
ARCH(1; 1, 2) -0.280921  0.123263  -2.27904 2.308e-002
ARCH(1; 2, 2)  0.020918  0.118460   0.17658 8.599e-001
GARCH(1; 1, 1)  0.924886  0.070220  13.17124 0.000e+000
GARCH(1; 2, 1)  0.055811  0.061759   0.90369 3.666e-001
GARCH(1; 1, 2)  0.094547  0.127677   0.74052 4.593e-001
GARCH(1; 2, 2)  0.902398  0.108385   8.32590 8.882e-016

--------------------------------------------------------------

AIC(11) = -6136.881
BIC(11) = -6090.324

Normality Test:
--------------------------------------------------------------
   Jarque-Bera    P-value Shapiro-Wilk  P-value
v1       13.33 1.272e-003       0.9766 0.006411
v2       26.14 2.112e-006       0.9794 0.044397

Ljung-Box test for standardized residuals:
--------------------------------------------------------------
   Statistic P-value Chi^2-d.f.
v1    15.397  0.2204         12
v2     8.404  0.7528         12

Ljung-Box test for squared standardized residuals:
--------------------------------------------------------------
   Statistic P-value Chi^2-d.f.
v1    10.700  0.5548         12
v2     4.404  0.9750         12

Lagrange multiplier test:
--------------------------------------------------------------
    Lag 1  Lag 2   Lag 3   Lag 4  Lag 5   Lag 6  Lag 7   Lag 8  Lag 9  Lag 10  Lag 11    Lag 12       C
v1 0.5538 1.7195 -0.7045 -1.0199 0.2528  0.4327 0.4468 0.72121 1.5802 -1.0795 -0.3441 -0.007436  0.5234
v2 0.2656 0.9985  0.2316  0.4073 1.0403 -0.6377 1.0230 0.07003 0.0611 -0.9343 -0.3995 -0.367546 -0.8732

     TR^2 P-value F-stat P-value
v1 10.223  0.5964 0.9489  0.6026
v2  5.676  0.9315 0.5220  0.9648

epoh 发表于 2013-2-23 14:57
module(finmetrics)
data=importData(file="pitaya1.xls", type="EXCEL")
pitaya1=cbind(data$v1,da ...

我的程序代码没有问题啊，为什么您能做出来呢

火龙果汽水 发表于 2013-2-23 14:58
我的程序代码没有问题啊，为什么您能做出来呢

如果你是持续用同一个command window
请注意不要有相同的对象名称
譬如第一个程序是用data,pitaya.bekk,…
那么第二个程序就用data1, ,pitaya1.bekk