```
# NOT RUN {
# separable model: spatial and temporal sill will be ignored
# and kept constant at 1-nugget respectively. A joint sill is used.
separableModel <- vgmST("separable",
space=vgm(0.9,"Exp", 147, 0.1),
time =vgm(0.9,"Exp", 3.5, 0.1),
sill=40)
# product sum model: spatial and temporal nugget will be ignored and kept
# constant at 0. Only a joint nugget is used.
prodSumModel <- vgmST("productSum",
space=vgm(39, "Sph", 343, 0),
time= vgm(36, "Exp", 3, 0),
k=15)
# 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 fitted. The spatial,
# temporal and jont nuggets are set to 0.
simpleSumMetricModel <- vgmST("simpleSumMetric",
space=vgm(20,"Lin", 150, 0),
time =vgm(20,"Lin", 10, 0),
joint=vgm(20,"Exp", 150, 0),
nugget=1, stAni=15)
# metric model
metricModel <- vgmST("metric",
joint=vgm(60, "Exp", 150, 10),
stAni=60)
# }
```

Run the code above in your browser using DataCamp Workspace