先mark一下  最近正在学MRS-GARCH模型  现在似然函数编出来的，却不知道该如何做参数估计，还好有楼主的帖子做个参考，实在感谢

你好，又要打扰楼主了，你发的程序我现在看了1/4多吧，但是又几个疑惑，是在等不急，就先问你一下吧：
首先是帖子第二楼的“start gibbs toop中，求了三个似然函数（09行——42行），ARCH效应中本期的条件方差是用前几期（这程序里是滞后阶数是3）的残差平方表示的，但是在程序中却是先把残差按波动的大小分为3个状态，而在3个状态中分别求条件方差h,那没在时间上，显然t期的h不是用t-1,t-2,t-3期的残差平方表示了。所以这里求h的方法让我无法理解，楼主能不能把原理讲一下，谢谢。
另一个问题仍是在GIBBS TOOP中，第49行有一个norm_rnd函数，这个是不是从正态分布中抽取随机数？因为我的matlab中现实没有这个函数，所以我就改成了normrnd（），但是又显示
Error using normrnd (line 28)
Requires at least two input arguments.

Error in code_ms_3s (line 147)
   beta_can = beta + normrnd(beta_sig);
我不知道问题出在哪里，还是希望楼主能指点一下，谢谢

49行可不可以改成
beta_can=beta+[normrnd(0,(0.003/4)^2);normrnd(0,(0.013/4)^2);normrnd(0,0.018/4)^2)]

不好意思，过去几个月在工作。很少来论坛交流，也有7个月的时间做客户营销，没时间运行和修改matlab程序；
你问的问你很深入，也比较难以解决。毕竟MR-GARCH不是简单的。Griddy gibbs 抽样技术，也许你们自己有书，我从论坛找到国外的一个作者的贝叶斯金融统计建模，里面也提到过这个。不过很多我还没理解。
关于你提到的 normrnd()函数，这个是出自 lesage 的 JPL7 计量经济学工具箱里，是他本人自己写的一个多元正态分布随机函数，你可以下载JPL7，将这里的代码放进去，就可以正常运行。
如果你改用系统的多元正态分布随机数函数，因为输入参数和输出，也就是我们说的函数调用不匹配问题，导致出错。

明天要考试。所以不深入聊下去。我很高兴你这样深入的看这个程序和提这么多问题。有时间的时候多交流。
我的Q：280201722 潇湘品茗 30岁，大学毕业5年了。

楼主，我觉得“start gibbs toop”中的似然函数（第9——42行）可以改成这样，会和模型更符合一些
beta_by_s=zeros(n,1);
for t=1:n
    if S(t)==1
       beta_by_s(t)=beta(1,1);
    else
       if S(t)==2
          beta_by_s(t)=beta(2,1);
       else
          beta_by_s(t)=beta(3,1);
       end
    end
end
err_by_s=y-beta_by_s;
h_by_s=zeros(n,1);
   if S(1)==1
      h_by_s(1)=a_1(1)/(1-a_1(2)-a_1(3)-a_1(4));
   else
      if S(2)==2
         h_by_s(1)=a_2(1)/(1-a_2(2)-a_2(3)-a_2(4));
      else
         h_by_s(1)=a_3(1)/(1-a_3(2)-a_3(3)-a_3(4));
      end
   end

   if S(2)==1
      h_by_s(2)=a_1(1)+a_1(2)*err_by_s(1)^2;
   else
      if S(2)==2
         h_by_s(2)=a_2(1)+a_2(2)*err_by_s(1)^2;
      else
         h_by_s(3)=a_3(1)+a_3(2)*err_by_s(1)^2;
      end
   end

   if S(3)==1
      h_by_s(3)=a_1(1)+a_1(2)*err_by_s(2)^2+a_1(3)*err_by_s(1)^2;
   else
      if S(3)==2
         h_by_s(3)=a_2(1)+a_2(2)*err_by_s(2)^2+a_2(3)*err_by_s(1)^2;
      else
         h_by_s(3)=a_3(1)+a_3(2)*err_by_s(2)^2+a_3(3)*err_by_s(1)^2;
      end
   end

  for t=4:n
      if S(t)==1
         h_by_s(t)=a_1(1)+a_1(2)*err_by_s(t-1)^2+a_1(3)*err_by_s(t-2)+a_1(4)*err_by_s(t-3);
      else
        if S(t)==2
           h_by_s(t)=a_2(1)+a_2(2)*err_by_s(t-1)^2+a_2(3)*err_by_s(t-2)+a_2(4)*err_by_s(t-3);
        else
           h_by_s(t)=a_3(1)+a_3(2)*err_by_s(t-1)^2+a_3(3)*err_by_s(t-2)+a_3(4)*err_by_s(t-3);
        end
      end
  end

  lterm_by_s=(err_by_s.^2)./h_by_s;
  h_s1=h_by_s(S==1);
  h_s2=h_by_s(S==2);
  h_s3=h_by_s(S==3);
  lterm_s1=lterm_by_s(S==1);
  lterm_s2=lterm_by_s(S==2);
  lterm_s3=lterm_by_s(S==3);
  llike1=-0.5*sum(log(h_s1))-0.5*sum(lterm_s1);
  llike2=-0.5*sum(log(h_s2))-0.5*sum(lterm_s2);
  llike3=-0.5*sum(log(h_s3))-0.5*sum(lterm_s3);

我对MRS-GARCH了解有限，上面的程序是我自己的看法，不知道对不对。

请问48.beta_sig = [(0.003/4)^2 0 0; 0 (0.013/4)^2 0; 0 0 (0.018/4)^2 ];是如何设定出来的？

楼主能否请教你markov的问题呢？？？？不胜感激

这个好像是单元的，有木有人，有二元或者是多元的MS-GARCH程序。欢迎多讨论，一起多交流。我的扣扣：529820052

thanks a lot