surveillance (version 1.12.1)

plot.atwins: Plot results of a twins model fit

Description

Plot results of fitting a twins model using MCMC output. Plots similar to those in the Held et al. (2006) paper are generated

Usage

## S3 method for class 'atwins':
plot(x, which=c(1,4,6,7), ask=TRUE, ...)

Arguments

x

An object of class atwins.

which

a vector containing the different plot types to show [object Object],[object Object],[object Object],[object Object],[object Object],[object Object],[object Object]

ask

Boolean indicating whether to ask for a newline before showing the next plot.

...

Additional control for the plots, which are currently ignored.

Value

• This function does not return anything.

encoding

latin1

Details

For details see the plots in the paper. Basically MCMC output is visualized. This function is together with algo.twins still experimental.

References

Held, L., Hofmann, M., H�hle{Hoehle}, M. and Schmid V. (2006) A two-component model for counts of infectious diseases, Biostatistics, 7, pp. 422--437.

See Also

algo.twins

Examples

#Apparently, the algo.atwins can crash on some LINUX systems
#thus for now the example section is commented

#Load the data used in the Held et al. (2006) paper
data("hepatitisA")

#Fix seed - this is used for the MCMC samplers in twins
set.seed(123)

#Call algorithm and save result
otwins <- algo.twins(hepatitisA)

#This shows the entire output
plot(otwins,which=c(1,2),ask=FALSE)