ox报错

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位老师，ox做msvar跑程序出错的结果，该怎么办？
Ox version 3.40 (Windows) (C) J.A. Doornik, 1994-2004
MSVAR (c) H-M Krolzig, 1996-2005, package version 1.32a, object created on  6-11-2002
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarEm.ox (950): mle_msmh
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarEm.ox (951): mle_msmh
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarEm.ox (950): mle_msmh
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarEm.ox (951): mle_msmh
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarFilter.ox (36): filter
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarFilter.ox (36): filter
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarFilter.ox (36): filter
Warning: invert: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarEm.ox (541): Estimate
Warning: choleski: decomposition failed
C:\hmk\MSVAR\INCLUDE\hmk.ox (410): vechcholeskiblock
Warning: choleski: decomposition failed
C:\hmk\MSVAR\INCLUDE\hmk.ox (410): vechcholeskiblock
Warning: choleski: decomposition failed
C:\hmk\MSVAR\INCLUDE\hmk.ox (410): vechcholeskiblock
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarCom.ox (428): GetErrors
Warning: choleski: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarCom.ox (428): GetErrors
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarCom.ox (428): GetErrors
Warning: choleski: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarCom.ox (428): GetErrors
Warning: invertsym: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarCom.ox (428): GetErrors
Warning: choleski: decomposition failed
C:\hmk\MSVAR\INCLUDE\MsvarCom.ox (428): GetErrors


---------- Error occured! EM algorithm was stopped after  2 iterations ------------

*** Warning: Some transition probabilities are close to the border;
             Numerical stability endangered.

Runtime error: '[70 : 153][] in matrix[79][1]' invalid index range
Runtime error occurred in decomposeLambda(111), call trace:
C:\hmk\MSVAR\INCLUDE\MsvarCom.ox (111): decomposeLambda
C:\hmk\MSVAR\INCLUDE\MsvarEm.ox (286): fLogLik
C:\hmk\MSVAR\INCLUDE\MsvarEm.ox (363): Covar
C:\hmk\MSVAR\INCLUDE\MsvarEm.ox (381): GetStdErr
C:\hmk\MSVAR\INCLUDE\MsvarEm.ox (411): StdErr
C:\Program Files\Ox\package\msvar\ms3var1.ox (17): main

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：composition Transition iterations Abilities Stability 950 951 Windows version created

如果数据太大
就截成小部分上传
只要能运行说明错误即可

shuju.xls (264.5 KB)

2012-4-25 15:52:23 上传

数据

eoph老师：
    您好！下面是我的数据，是excel2000格式的。

helen_80223 发表于 2012-4-25 15:53
eoph老师：
    您好！下面是我的数据，是excel2000格式的。

你的数据使用msvar->SetModel(MSMH, 3);
的确有问题
若改为msvar->SetModel(MSIH, 3),则结果如下:
shuju_revised.rar (86.63 KB) 本附件包括：

• shuju_revised.xls

2012-4-25 18:43:04 上传

msvar->Load("shuju.xls");
decl p=1;  
msvar->Select(Y_VAR, {"bb",0,p,"ss",0,p,"ff",0,p,"cc",0,p,"gg",0,p,"ee",0,p,"eee",0,p});
msvar->SetSample(1, 1, 1568,1);      
msvar->SetModel(MSIH, 3);

---------- EM algorithm converged after 25 iterations ------------

EQ( 1) MSIH(3)-VAR(1) model of (bb,ss,ff,cc,gg,ee,eee)
       Estimation sample: 2 - 1568

no. obs. per eq. :       1567    in the system :      10969   
no. parameters   :        160    linear system :         84   
no. restrictions :         70
no. nuisance p.  :          6

log-likelihood   :  3882.2544    linear system :  1640.6422  

AIC criterion    :    -4.7508    linear system :    -1.9868
HQ  criterion    :    -4.5475    linear system :    -1.8800
SC  criterion    :    -4.2038    linear system :    -1.6996

LR linearity test:   4483.2245    Chi(70) =[0.0000] **  Chi(76)=[0.0000] **  DAVIES=[0.0000] **  


---------- matrix of transition probabilities ------

          Regime 1  Regime 2  Regime 3
Regime 1    0.8120    0.1795    0.0085
Regime 2    0.0725    0.8716    0.0559
Regime 3    0.0294    0.4330    0.5375


---------- regime properties ----------------------

              nObs     Prob.  Duration
Regime 1     416.9    0.2646      5.32
Regime 2     1019.1    0.6518      7.79
Regime 3     130.9    0.0836      2.16

---------- coefficients ----------------------------

                     bb         ss         ff         cc         gg         ee        eee
Const(Reg.1)   0.305202  -0.002169   0.002145   0.001591   0.060467   8.343099  -0.003348
Const(Reg.2)   0.561997  -0.003452   0.003867   0.003228   0.037343   7.691627   0.035148
Const(Reg.3)   0.011184  -0.003635   0.003890   0.003324  -0.006087   8.322107  -0.020444
bb_1           0.998721   0.000024  -0.000013  -0.000005  -0.000019   0.010895   0.000038
ss_1          -0.027467   0.998308   0.000183  -0.000518  -0.024502   0.143958  -0.001351
ff_1         -21.637451  -0.129114  -0.050704   0.005473  -3.021816 -33.296722  -0.939461
cc_1          47.507210   0.000440   0.036738  -0.059318   6.272398  18.950272  -0.127948
gg_1           0.001605   0.000296  -0.000089  -0.000064   0.998824   0.193847   0.001997
ee_1          -0.002364   0.000007  -0.000022  -0.000010  -0.000293   0.857617  -0.001096
eee_1         -0.003546  -0.000002   0.000006  -0.000028   0.001440  -0.018447   1.002541
  SE (Reg.1)   3.057883   0.048452   0.021299   0.018070   0.281792  10.431529   0.423327
  SE (Reg.2)   1.906766   0.011751   0.008685   0.009916   0.311347   5.070088   0.283263
  SE (Reg.3)   4.651493   0.011246   0.012597   0.014788   0.830067   8.690766   2.418986

---------- contemporaneous correlation -------------

Regime 1
           bb        ss        ff        cc        gg        ee       eee
bb     1.0000   -0.1612    0.5412    0.6497    0.1714    0.0273    0.3634
ss    -0.1612    1.0000   -0.1066   -0.0960   -0.0078    0.0925   -0.1438
ff     0.5412   -0.1066    1.0000    0.6375    0.0622    0.0221    0.3811
cc     0.6497   -0.0960    0.6375    1.0000    0.0627   -0.0082    0.3984
gg     0.1714   -0.0078    0.0622    0.0627    1.0000   -0.1066    0.0910
ee     0.0273    0.0925    0.0221   -0.0082   -0.1066    1.0000   -0.0414
eee    0.3634   -0.1438    0.3811    0.3984    0.0910   -0.0414    1.0000

Regime 2
           bb        ss        ff        cc        gg        ee       eee
bb     1.0000   -0.0765    0.3214    0.5205    0.1386   -0.0003    0.1808
ss    -0.0765    1.0000   -0.1500   -0.0840   -0.0508    0.0510   -0.0146
ff     0.3214   -0.1500    1.0000    0.3647    0.0303    0.0352    0.1060
cc     0.5205   -0.0840    0.3647    1.0000    0.1886    0.0198    0.1911
gg     0.1386   -0.0508    0.0303    0.1886    1.0000    0.0048    0.0875
ee    -0.0003    0.0510    0.0352    0.0198    0.0048    1.0000    0.0247
eee    0.1808   -0.0146    0.1060    0.1911    0.0875    0.0247    1.0000

Regime 3
           bb        ss        ff        cc        gg        ee       eee
bb     1.0000   -0.0523   -0.1865    0.6153    0.1845   -0.0066    0.1021
ss    -0.0523    1.0000    0.0402   -0.0135    0.1151   -0.0406    0.0128
ff    -0.1865    0.0402    1.0000   -0.0132   -0.0277   -0.0617    0.0211
cc     0.6153   -0.0135   -0.0132    1.0000    0.1933   -0.0871    0.0848
gg     0.1845    0.1151   -0.0277    0.1933    1.0000    0.2774    0.0840
ee    -0.0066   -0.0406   -0.0617   -0.0871    0.2774    1.0000    0.0540
eee    0.1021    0.0128    0.0211    0.0848    0.0840    0.0540    1.0000

helen_80223 发表于 2012-4-25 15:53
eoph老师：
    您好！下面是我的数据，是excel2000格式的。

如果你一定要用
msvar->SetModel(MSMH, 3);
试着把ff*100,cc*100

epoh老师：
    您好！
    我现在的情况是这样，我的程序跑出来了，也不抱错了，可是我不知道到底哪个样本是属于哪个机制下的。这个需要怎么弄呀？非常感谢！

epoh老师：
     您好！还有那个t检验值我的报告里面没有，这个怎么弄出来呢！

友情帮顶

helen_80223 发表于 2012-4-25 19:28
epoh老师：
     您好！还有那个t检验值我的报告里面没有，这个怎么弄出来呢！

先声明一下,上面结果还没跑完
还要一点时间time passed:  5:08.70

结果说明请参考
msvardoc.pdf
page 12 9.2 AnMS-VAR model of international business cycles

page 14 9.3 A Markov-switching vector equilibrium correction model

/////
---------- standard errors -------------------------

                      bb          ss          ff          cc          gg          ee         eee
Const(Reg.1)     0.48268    0.013058     0.25258     0.24168    0.059507     0.53314    0.071351
Const(Reg.2)     0.37601   0.0090159     0.18214     0.18519    0.050807     0.35157    0.056717
Const(Reg.3)     0.63801   0.0094241     0.25091     0.26563    0.098978     0.58407     0.22479
bb_1           0.0017937 6.4913e-005  0.00078443  0.00083777  0.00025382   0.0043895  0.00025897
ss_1             0.14396   0.0015729    0.085650    0.079615    0.016720     0.39270    0.021450
ff_1            0.052821  0.00039494    0.027908    0.027822   0.0066410     0.15314   0.0079682
cc_1            0.053792  0.00036071    0.025778    0.027334   0.0072340     0.14977   0.0079809
gg_1            0.019575  0.00012023   0.0086879   0.0095479   0.0029205    0.050919   0.0028644
ee_1           0.0046357 5.6192e-005   0.0023253   0.0024046  0.00063564    0.010986  0.00068029
eee_1          0.0077596 5.0978e-005   0.0036440   0.0040300   0.0012072    0.021217   0.0011721

---------- t - values ------------------------------

                    bb        ss        ff        cc        gg        ee       eee
Const(Reg.1)    0.6363   -0.1677    0.8400    0.6551    1.0165   15.6483   -0.0505
Const(Reg.2)    1.5008   -0.3843    2.1165    1.7420    0.7398   21.8801    0.6172
Const(Reg.3)    0.0237   -0.3857    1.5441    1.2465   -0.0654   14.2502   -0.0922
bb_1          556.7855    0.3655   -1.6064   -0.5620   -0.0836    2.4754    0.1486
ss_1           -0.1904  634.6933    0.2243   -0.6472   -1.4676    0.3682   -0.0593
ff_1           -4.0977   -3.2703   -1.8241    0.1945   -4.5526   -2.1718   -1.1782
cc_1            8.8260    0.0182    1.4256   -2.1742    8.6708    1.2686   -0.1653
gg_1            0.0733    2.4556   -1.0206   -0.6738  341.9984    3.8070    0.6948
ee_1           -0.5091    0.1354   -0.9308   -0.4242   -0.4551   78.0641   -1.6109
eee_1          -0.4591   -0.0417    0.1537   -0.6995    1.1947   -0.8692  855.3119


****        time passed:  5:08.70        ****

epoh老师：
     您的意思是需要等一段时间才能出来这个t值吗？非常感谢