Я использую пакет R
StMoMo
для данных о смертности, но у меня возникли некоторые проблемы с пониманием процедуры начальной загрузки.Как я понял, функция boostrap
позволяет создавать параметры начальной загрузки, чтобы учесть их неопределенность (поскольку вычислить неопределенность модели классическим способом непросто).После получения nBoot
загрузочных выборок необходимо использовать функцию simulate
для выборки из улучшенных параметров и прогноза.В списке аргументов этой функции есть nsim
, то есть количество образцовых путей, которые будут сгенерированы из параметров Boostrapped, но в виньетке пакета, Villegas et al.не использовал аргумент nsim, что означает, что только 1 образец был сгенерирован из каждого параметра начальной загрузки.Итак, правильно ли увеличивать количество выборок, чтобы иметь лучший доверительный интервал?Ниже приведены два разных результата, фиксирующих nBoot
до 100 и варьирующиеся nsim
от 1 до 100.
100 с начальной загрузкой nsim
= 1
100начальная загрузка с nsim
= 100
Как мы видим во втором рис. доверительный интервал кажется лучше, чем в первом (хорошо, цель должна состоять в том, чтобы построить доверительный интервал, включающий красную точку),Мой вопрос может быть немного запутанным, но в любом случае, каково правильное количество симуляций из n
образцов начальной загрузки?Спасибо