element_Lhood_poisson: # Likelihood of observing outbreaks of size y given...

Description Usage Arguments Details Value Author(s) Examples

View source: R/element_Lhood_poisson.R

Description

The function internally relies on pre-computing: 1) g: the probability of observing size y given z,rho; with z the true number of cases (reported and un-reported cases) and rho the reporting probability

Usage

1
element_Lhood_poisson(R, z, g, g0)

Arguments

R

is the reproduction number, i.e. the average number of secondary cases due to a single case. This can be any positive number.

z

is the true potential outbreak sizes. It is precomputed by proba_observation.

g:

the probability of observing size y given z,rho; with z the true number of cases (reported and un-reported cases) and rho the reporting probability. It is precomputed by proba_observation.

g0:

the probability of not observing an outbreak (i.e. y=0 or no cases) given z and rho. It is precomputed by proba_observation.

Details

2) g0: the probability of not observing an outbreak (i.e. y=0 or no cases) given z and rho.

g and g0 are obtained using the function proba_observation.

Value

The function returns the likelihood of the observations for a given R,rho.

Author(s)

Pierre Nouvellet (p.nouvellet@imperial.ac.uk)

Examples

1
2
x <- element_Lhood_poisson(.5,x$possible_size,x$p_y_z,x$p_0_z)
x

reconhub/branchr documentation built on May 27, 2019, 4:01 a.m.