sample_sequences: Sample Sequence Observation Model

In epiforecasts/forecast.vocs: Forecast Case and Sequence Notifications using Variant of Concern Strain Dynamics

Sample Sequence Observation Model

Description

Sample Sequence Observation Model

Usage

sample_sequences(frac_voc, seq_total, phi)

Arguments

frac_voc	A numeric vector of expected proportions positive for the variant of concern.
seq_total	An integer vector of total sequences available.
phi	The overdispersion of the sampling process. If not supplied then no overdispersion is used (i.e a binomial observation model vs a beta binomial observation model).

Value

A vector of observed sequences positive for the variant of concern.

See Also

Functions to generate simulated data generate_obs()

Examples

# dummy sequence data
frac_voc <- seq(0, 1, by = 0.1)
seq_total <- seq(10, length.out = length(frac_voc), by = 100)

# binomial observation model
sample_sequences(frac_voc, seq_total)

# beta binomial observation model
sample_sequences(frac_voc, seq_total, 0.5)

epiforecasts/forecast.vocs documentation built on May 14, 2023, 2 p.m.