求助：关于MSVAR程序包

epoh 发表于 2012-4-16 15:32
Q1:_n_x, _n_z,只有在typmod=4 typmod=5 才用得上
      (typmod==4);        
        /* ===   MS- ...

老师你好：
我现在用我的论文数据做MSVAR时候，如果选择了_typmod=3 那么一旦运行就出现下面的错误提示，如果选_typmod=1 则可以运行。我想请问下老师是为什么？我想选择3做，有什么办法呢？附件是我的数据和修改的程序。请老师帮忙看看谢谢！
C:\Gauss\MSVAR\MSlib\MSVAR_VecMat.prg(167) : error G0048 : Matrix singular
Currently active call: MSVAR_Ergodic [167] C:\Gauss\MSVAR\MSlib\MSVAR_VecMat.prg
Stack trace:
   MSVAR_Ergodic called from C:\Gauss\MSVAR\MSlib\MSVAR_VecMat.prg, line 43
   MSVAR_Vec_Mat called from C:\Gauss\MSVAR\MSlib\MSVAR_FiltHmm.prg, line 41
   MSVAR_Filt called from C:\Gauss\MSVAR\MSlib\MSVAR_MaxHmm.prg, line 40
   MSVAR_likfcn called from c:\gauss\src\optutil.src, line 445
   _optmum called from c:\gauss\src\optmum.src, line 264
   optmum called from C:\Gauss\MSVAR\MSlib\MSVAR_MaxHmm.prg, line 64
   MSVAR_MaxHmm called from C:\gauss\msvar\MSlib\MSVAR.prg, line 336
   MSVAR_Estimate called from C:\gauss\msvar\MSlib\MSVAR.prg, line 214

ferguson06 发表于 2012-5-4 10:36
老师你好：
我现在用我的论文数据做MSVAR时候，如果选择了_typmod=3 那么一旦运行就出现下面的错误提示 ...

这是你的数据所造成
你的数据应该要取log
or 100*log

可以参考KROTO.xls,GNP82.xls,WBC.xls

epoh 发表于 2012-5-4 14:47
这是你的数据所造成
你的数据应该要取log
or 100*log

好的 谢谢老师 我去试试

epoh 发表于 2012-4-15 21:13
Date:  4/15/12   Initial time: 21:05:09

Start and ending period

谢谢你们的讨论，我也学到了很多东西，不过我对结果还是很糊涂，向大家请教下，比如你们的例子
==============Regime      1.00000000=============

Beta
      -0.41293422      1.09166441

Delta

      0.24226903     -0.35544298
      -0.20933353      0.64755363

Sigma

      0.65384647      0.00000000
       0.00000000      0.11022738

==============Regime      2.00000000=============

Beta
      -0.19425370      0.46819728

Delta

      0.24226903     -0.35544298
      -0.20933353      0.64755363

Sigma

      0.07106441      0.00000000
       0.00000000      0.08369711
------------------------------------------------------------------
这里是_typmode=3, 所以公式是y(t)= mu(S(t)) + [ yt-1,..yt-p]*Delta(t) +u(t)  = 1.*Beta_S(t)+z(t).Delta(t)+ u(t)

m=2, k=2, p=1 只是均值变，因为k=2,所以每个regime里有2个均值，一共4个beta变量，我不明白的是为什么在结果里delta参数在每个regime里有4个？我以为是2个，分别对应y1t-1和y2t-1，在不同的regime，这两个值是不变的，但这里为什么是4个，这个结果应该怎么看, 对应程序运行结果这里delta1,2分别是什么，还是我对msvar模型理解错了
regime1
y1t=-0.41293422 +y1t-1*Delta1+u(t)
y2t=1.09166441+y2t-1*Delta2+u(t)

Regime2
y1t= -0.19425370  +y1t-1*Delta1+u(t)
y2t=0.46819728+y2t-1*Delta2+u(t)
谢谢

为什么我直接运行会得到下面的结果？
(0) : error G0008 : '#lineson' : Syntax error
(0) : error G0008 : '#include C:\Gauss\MSVAR\MSLIB\MSVAR_Call.prg' : Syntax error

我建了个群255607063 大家可以讨论下啊

ding