R rcpp错误：如何禁用估计值的输出

Question

我现在正在使用rcppbugs来优化生态系统模型的参数。情况是我的观察数据是一个非常长的时间序列，它至少需要10,000次迭代。由于run.model()每次都会输出估算值（此处为 y .hat ），因此会出现一个大问题，例如：

List of 6
 $ factor1 : num [1:2] 0.921 0.921
 $ factor2 : num [1:2] 2 2
 $ factor3 : num [1:2] 0.954 0.954
 $ tau.y   : num [1:2] 0.17 0.17
 $ y.hat   : num [1:2, 1:70176] 0.139 0.139 0.138 0.138 0.136 ...
 $ y.lik   : NULL
 - attr(*, "acceptance.ratio")= num 0

显然每次都很难保存，但同时我需要记录所有factor1，factor2和factor3值。有人可以帮忙吗？感谢。

Answer 1

如果通过R message或warning进行打印，则可以分别使用suppressMessages和suppressWarnings。如果使用cat，print生成输出或从C / C ++ / Fortran函数发送到“std out”，则可以使用capture.output。

Answer 2

y.hat变量由rcppbugs C ++源代码中的getHistory()函数输出。因此，修改此功能足以解决上述问题。