Plot Generation of the Observed and the Defined Outbreak States of a (Multivariate) Time Series
Plotting a (multivariate) disProg
object. The internal
function plot.disProg.one
is used as a helper function
to plot a univariate time series.
## S3 method for class 'disProg' plot(x, title = "", xaxis.years=TRUE, startyear = x$start[1], firstweek = x$start[2], as.one=TRUE, same.scale=TRUE, ...) plot.disProg.one(x, title = "", xaxis.years=TRUE, quarters=TRUE, startyear = x$start[1], firstweek = x$start[2], ylim=NULL, xlab="time", ylab="No. infected",type="hh",lty=c(1,1),col=c(1,1), outbreak.symbol = list(pch=3, col=3), legend.opts=list(x="top", legend=c("Infected", "Outbreak"), lty=NULL,pch=NULL,col=NULL), ...)
x |
object of class |
title |
plot title |
xaxis.years |
if |
quarters |
add quarters to the plot |
startyear |
year to begin the axis labeling (the year where
the oldest data come from). This arguments will be obsolete in |
firstweek |
number of the first week of January in the first year (just for axis labeling grounds) |
as.one |
if |
same.scale |
if |
ylim |
range of y axis |
xlab |
label of the x-axis |
ylab |
label of the y-axis |
type |
line type of the observed counts (should be |
lty |
line type of the observed counts |
col |
color of the observed count lines |
outbreak.symbol |
list with entries |
legend.opts |
a list containing the entries to be sent to the
An further arguments to the |
... |
arguments passed to |
a plot showing the number of infected and the defined alarm status for a time series created by simulation or given in data either in one single plot or in several plots for each individual time series.
M. Höhle with contributions by A. Riebler and C. Lang
# Plotting of simulated data disProgObj <- sim.pointSource(p = 0.99, r = 0.5, length = 208, A = 1, alpha = 1, beta = 0, phi = 0, frequency = 1, state = NULL, K = 5) # plot the simulated disease with the defined outbreaks plot(disProgObj) title <- "Number of Infected and Defined Outbreak Positions for Simulated Data" plot(disProgObj, title = title) plot(disProgObj, title = title, xaxis.years=TRUE, startyear = 1999, firstweek = 13) plot(disProgObj, title = title, xaxis.years=TRUE, startyear = 1999, firstweek = 14) # Plotting of measles data data(measles.weser) # one plot plot(measles.weser, title = "measles cases in the district Weser-Ems", xaxis.years=TRUE, startyear= 2001, firstweek=1) # plot cases for each "Kreis" plot(measles.weser, same.scale=TRUE, as.one=FALSE)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.