getOutputTables: Summarize results of Bayesian demographic model in tables

View source: R/getOutputTables.R

getOutputTablesR Documentation

Summarize results of Bayesian demographic model in tables

Description

Produces summary tables for Bayesian caribou population model results. Optionally calculates Kolmogorov–Smirnov distances between Bayesian model results and national model results.

Usage

getOutputTables(
  caribouBayesDemogMod,
  startYear = min(caribouBayesDemogMod$inData$disturbanceIn$Year),
  endYear = max(caribouBayesDemogMod$inData$disturbanceIn$Year),
  paramTable = data.frame(param = "observed"),
  exData = NULL,
  simNational = NULL,
  getKSDists = FALSE
)

Arguments

caribouBayesDemogMod

caribou Bayesian demographic model results produced by calling caribouBayesianPM()

startYear, endYear

year defining the beginning of the observation period and the end of the projection period.

paramTable

data.frame. Optional. Scenario parameters see simulateObservations()

exData

data.frame. Optional. Output of simulateObservations() that records the true population metrics of the population that observations were simulated from.

simNational

National simulation results, produced by calling getSimsNational()

getKSDists

logical. Should Kolmogorov–Smirnov distances be calculated?

Value

a list of tables:

  • rr.summary.all: Mean parameter values for each year and standard deviation, upper and lower credible intervals projected by the Bayesian model, as well as scenario input parameters.

  • sim.all: Mean parameter values and upper and lower credible intervals from the National model for each year, as well as scenario input parameters.

  • obs.all: Observed parameter values with column "type" identifying if it is the "true" value of the simulated population or the "observed" value simulated based on the collaring program parameters.

  • kdDists: Kolmogorov–Smirnov distances and p-values, NA if getKSDists is false

See Also

Caribou demography functions: caribouBayesianPM(), caribouPopGrowth(), compositionBiasCorrection(), demographicCoefficients(), demographicProjectionApp(), demographicRates(), getPriors(), getScenarioDefaults(), getSimsNational(), plotRes(), popGrowthTableJohnsonECCC, runScnSet(), simulateObservations()

Examples

scns <- getScenarioDefaults(projYears = 10, obsYears = 10,
                            obsAnthroSlope = 1, projAnthroSlope = 5,
                            collarCount = 20, cowMult = 5)

simO <- simulateObservations(scns)

out <- caribouBayesianPM(survData = simO$simSurvObs, ageRatio = simO$ageRatioOut,
                          disturbance = simO$simDisturbance,
                          Nchains = 1, Niter = 100, Nburn = 10,
                          Nthin = 2)

getOutputTables(out, exData = simO$exData, paramTable = simO$paramTable,
                simNational = getSimsNational(), getKSDists = FALSE)

LandSciTech/caribouMetrics documentation built on Feb. 25, 2025, 9:34 a.m.