R rcppbugs: как отключить вывод оценочных значений

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 

Очевидно, что трудно сохранить y.hat каждый раз, но пока мне нужно запишите все значения factor1, factor2 и factor3. Может ли кто-нибудь помочь? Благодарю.

Answer 1

Переменная y.hat выводится функцией getHistory() в исходном коде rcppbugs C++. Таким образом, достаточно изменить эту функцию для решения вышеуказанной проблемы.

Answer 2

Если печать осуществляется через R's message или warning, вы можете использовать suppressMessages и suppressWarnings, соответственно. Если вывод производится с использованием cat, print или отправляется на «std out» из функции C/C++/Fortran, вы можете использовать capture.output.