if (requireNamespace("maptools")) {
# load disProg-object "ha" and convert to S4-class "sts"
data("ha")
shpfile <- system.file("shapes/berlin.shp",package="surveillance")
ha.sts <- disProg2sts(ha, map=maptools::readShapePoly(shpfile,IDvar="SNAME"))
} else {
data("ha.sts")
# is almost identical to the above except that German umlauts
# have been replaced in 'ha.sts@map@data$BEZIRK' for compatibility reasons
}
ha.sts
plot(ha.sts, type=observed ~ 1 | unit)
#Convert ts/mts object to sts
z <- ts(matrix(rpois(300,10), 100, 3), start = c(1961, 1), frequency = 12)
sts.z <- as(z,"sts")
plot(sts.z)
Run the code above in your browser using DataLab