data("ha.sts")
# Show a (subset of a) single time series
plot(ha.sts[,7])
plot(ha.sts[year(ha.sts)==2006, 7])
# Map a single time point
plot(ha.sts[5*52+26,], type=observed~unit)
plot(ha.sts, type=observed~unit, tps=5*52+26) # same -> ?stsplot_space
# Restrict the data (and the map) to a subset of the districts
plot(ha.sts[,c("pank","lich")], type=observed~unit, labels=TRUE)
plot(ha.sts[,c("pank","lich"),drop=TRUE], type=observed~unit, labels=TRUE)
Run the code above in your browser using DataLab