r - R 包“MSwM”的结果不一致

Question

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"

请帮我..

Answer 1

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

其次，就估计而言，我发现这两种制度的两种结果之间似乎存在“转换”——这（对我而言）是一个非常有趣的观察结果，我不确定是否存在任何潜在的隐藏关系。