Map of Disease Incidence
For each period (row) or for the overall period of the
observed
matrix of the "sts"
object, a map
showing the counts by region is produced.
It is possible to redirect the output into files, e.g., to generate an
animated GIF.
stsplot_spacetime(x, type, legend = NULL, opts.col = NULL, labels = TRUE, wait.ms = 250, cex.lab = 0.7, verbose = FALSE, dev.printer = NULL, ...)
x |
an object of class |
type |
a formula (see |
legend |
An object of type
If |
opts.col |
A list containing the two elements
|
labels |
|
wait.ms |
Number of milliseconds to wait between each plot |
cex.lab |
|
verbose |
|
dev.printer |
Either |
... |
Extra arguments sent to the plot function. |
The animate.sts
method provides a
re-implementation and supersedes this function!
Michael Höhle
Other stsplot
types, and animate.sts
for
the new implementation.
data("ha.sts") print(ha.sts) ## map of total counts by district plot(ha.sts, type=observed ~ 1 | unit) ## only show a sub-period total for two selected districts plot(ha.sts[1:20,1:2], type=observed ~ 1 | unit) ## Not run: # space-time animation plot(aggregate(ha.sts,nfreq=13), type= observed ~ 1 | unit * time) #print the frames to a png device #and do the animation without extra sleeping between frames imgname <- file.path(tempdir(), "berlin") plot(aggregate(ha.sts,nfreq=13), type = observed ~ 1 | unit * time, wait.ms=0, dev.printer=list(name=imgname)) #Use ImageMagick (you might have to adjust the path to 'convert') system(paste0("convert -delay 50 ", imgname, "*.png ", imgname, "-animated.gif")) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.