## Not run:
# # obtained by:
# library(spacetime)
# library(gstat)
# data(air)
#
# if (!exists("rural"))
# rural = STFDF(stations, dates, data.frame(PM10 = as.vector(air)))
# rr = rural[,"2005::2010"]
# unsel = which(apply(as(rr, "xts"), 2, function(x) all(is.na(x))))
# r5to10 = rr[-unsel,]
# vv = variogram(PM10~1, r5to10, width=20, cutoff = 200, tlags=0:5)
# ## End(Not run)
Run the code above in your browser using DataLab