Download Package `gstat`
Transcript
70
vv
time= vgm(36, "Exp",
sill=41, nugget=17)
3, 0),
# sum metric model: spatial, temporal and joint nugget will be estimated
sumMetricModel <- vgmST("sumMetric",
space=vgm( 6.9, "Lin", 200, 3.0),
time =vgm(10.3, "Lin", 15, 3.6),
joint=vgm(37.2, "Exp", 84,11.7),
stAni=77.7)
# simplified sumMetric model, only a overall nugget is
# temporal and jont nuggets are set to 0.
simpleSumMetricModel <- vgmST("simpleSumMetric",
space=vgm(20,"Lin", 150,
time =vgm(20,"Lin", 10,
joint=vgm(20,"Exp", 150,
nugget=1, stAni=15)
fitted. The spatial,
0),
0),
0),
# metric model
metricModel <- vgmST("metric",
joint=vgm(60, "Exp", 150, 10),
stAni=60)
vv
Precomputed variogram for PM10 in data set air
Description
Precomputed variogram for PM10 in data set air
Usage
data(vv)
Format
data set structure is explained in variogramST.
Examples
## Not run:
# obtained by:
data(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)