Get probabilities of abundances, as the multinomial probabilities, for each species at each time step as predicted by a TS model.

get_abund_probabilities(
  subsetted_dataset_item,
  fitted_lda,
  fitted_ts,
  max_sims = NULL
)

Arguments

subsetted_dataset_item

list with (at least) element $full: unaltered dataset.

fitted_lda

LDA fit to either a subsetted or a full dataset

fitted_ts

TS fit to that LDA

max_sims

defaults NULL. If supplied, will randomly draw max_sims draws from the posterior and return probabilities just for those. Useful for speeding things up.

Value

list of matrices of abundance probabilities for each species (columns) for each timestep (rows), with one list for each draw from the posterior (or max_sims, if supplied).

Details

Gets one matrix of probabiltiies for each draw from the posterior for estimates of changepoint locations and model parameters (in this case just the intercept).