Я работаю со смертностью , прогнозирую и выполнил bootstrap, чтобы создать 2.000 полупараметрических c копий моих данных. Я использую пакет StMoMO - , где для данных о смертности доступна функция bootstrap. После запуска bootstrap -команды я получаю переменную со значением «bootParameters», которая представляет собой список 2000 репликаций. Каждая репликация содержит модель с 3 параметрами, в то время как я хочу сделать гистограмму по одному из коэффициентов.
Я использую следующий код:
# Creating 2000 semiparametric bootstrap samples
JA_lc_fitM_boot1 <- bootstrap(LCfit_JA_M, nBoot = 2000, type = "semiparametric")
# Creating the histogram over one of the parameters for value x = 30 (age)
JA_lc_fitM_boot1[["bootParameters"]][[1:2000]][["bx"]][[30]]
Но я получаю сообщение об ошибке, что
Ошибка рекурсивной индексации на уровне 3 ".
Однако, когда я просто смотрю на значение параметра для конкретной репликации c (например, 2 вместо 1: 2000), я получаю требуемое значение. Есть ли способ показать гистограмма более 2000 значений оценки параметра?
Заранее спасибо!
Андерс