我正在通过R中的BRugs
包在OpenBUGS中估算一个相当简单的模型。该模型只需要一个参数即可估算,因此计算很快,一切看起来都很好。但是,该模型是较大的仿真研究的一部分,在该研究中,模型被估算了数千次,并且运行需要几个小时。对于并行化,我将foreach-function
与%dopar%
运算符配合使用。
在仿真过程中,我有时会收到一条错误消息,指出“ BUGSHE〜1.EXE已停止工作”。但是,从R控制台和cpu工作负载来看,模拟仍在继续。但是,在某个时间点(计算时间为18小时之后的最后一次),R将停止foreach
循环,这是因为在过程中某处产生了NA。当我执行几十个而不是数千个模拟时,整个代码运行得非常好。
有人知道这个错误消息是什么意思吗?到目前为止,Google和SO上的搜索栏都没有帮助。
或者有任何提示可以提示我如何调查模拟错误的地方?在进行并行计算时,我很难找出任何错误...
很抱歉,我现在无法提供可复制的示例,也许我可以在接下来的几天中想到一个示例。