R包'MSwM'结果不一致

时间:2013-12-20 03:49:32

标签: r statistics time-series markov

  

mod.mswm< - msmFit(lm(y~x),k = 2,p = 1,sw = rep(T,4),control = list(maxiter = 700,parallel = F))   摘要(mod.mswm)

     

我在这个r包中有一个不一致的结果..

     首次运行 - >制度1 =“估计= 0.05”
             regime 2 = "estimate = 0.90"
     第二次运行 - >制度1 =“估计= 0.90”
              regime 2 = "estimate = 0.05"
     

请帮帮我..

1 个答案:

答案 0 :(得分:0)

首先,我建议查看该函数msmFit的R代码,并特别注意是否使用了任何随机数生成器以及之前是否有set.seed。根据我对其他软件包的体验,结果不一致可能是由于这个原因。

其次,我看到两个政权在估计方面的两个结果之间似乎存在“转换” - 这对我来说是一个非常有趣的观察,我不确定是否存在任何潜在的隐藏的关系。