Skip to contents

Get trajectories from a Bayesian model result

Usage

trajectoriesFromBayesian(
  bayesianResults,
  N0 = NULL,
  cPars = demographyDefaults(),
  returnSamples = TRUE,
  doSummary = TRUE,
  ...
)

Arguments

bayesianResults

A result from estimateBayesianRates

N0
cPars
returnSamples
doSummary
...

Examples

surv_data <- bboudata::bbousurv_a %>% filter(Year > 2010)
recruit_data <- bboudata::bbourecruit_a %>% filter(Year > 2010)
bbouInformative <- estimateBayesianRates(surv_data, recruit_data,
                                         return_mcmc = TRUE)

trajB <- trajectoriesFromBayesian(bbouInformative)
str(trajB, max.level = 1)
#> List of 5
#>  $ summary     :'data.frame':	60 obs. of  8 variables:
#>  $ samples     :'data.frame':	198000 obs. of  8 variables:
#>  $ surv_data   :'data.frame':	84 obs. of  8 variables:
#>  $ recruit_data:'data.frame':	6 obs. of  8 variables:
#>  $ popInfo     :'data.frame':	3000 obs. of  4 variables:
plotTrajectories(trajB)
#> Warning: Removed 210 rows containing missing values or values outside the scale range
#> (`geom_line()`).
#> Warning: Removed 6 rows containing missing values or values outside the scale range
#> (`geom_ribbon()`).
#> Warning: Removed 6 rows containing missing values or values outside the scale range
#> (`geom_line()`).