plot.incubestim {EpiLPS}

Plot the estimated incubation distribution

Description

This routine can be used to plot the estimated incubation distribution based on an object of class incubestim.

Usage

## S3 method for class 'incubestim'
plot(x, type = c("pdf", "cdf", "incubwin"), ...)

Arguments

x

An object of class incubestim.

type

The type of plot. Default is "pdf" for the plot of the density function. Setting it to "cdf" returns the cumulative distribution function and "incubwin" gives a bar plot showing the width of the incubation window.

...

Further arguments to be passed to plot.

Value

• The probability density function of the estimated incubation period with 95% credible envelope.

• The cumulative distribution function of the estimated incubation period with 95% credible envelope.

• A bar plot showing the width of the incubation windows.

Author(s)

Oswaldo Gressani [e-mail]

Examples

set.seed(123)
simdat <- incubsim(n = 30, tmax = 20) # Simulate incubation data
data <- simdat$Dobsincub              # Incubation bounds
incubfit <- estimIncub(x = data, niter = 500, tmax = 20, verbose = TRUE)
plot(incubfit)