Description Usage Arguments Value Author(s) See Also Examples
update
updates a geostatistical linear model based
on the given model.
1 2 3 4 5 6 7 8 |
object |
An object produced by the |
mod |
A spatial dependence model object obtained
from one of the |
... |
Not implemented. |
Returns an object of the same class as
object
.
Joshua French
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 | # generate response
y = rnorm(10)
# generate coordinates
x1 = runif(10); x2 = runif(10)
# data frame for observed data
data = data.frame(y, x1, x2)
coords = cbind(x1, x2)
psill = 2 # partial sill
r = 4 # range parameter
evar = .1 # error variance
fvar = .1 # add finescale variance
# one can't generally distinguish between evar and fvar, but
# this is done for illustration purposes
cmod_std = cmod_std("exponential", psill = psill, r = r,
evar = evar, fvar = fvar)
cmod_std2 = cmod_std("exponential", psill = psill + 1,
r = r + .5, evar = evar + .01,
fvar = fvar)
# check geolm update for universal kriging
gear1 = geolm(y ~ x1 + x2, data = data, mod = cmod_std,
coordnames = c("x1", "x2"))
gear2 = geolm(y ~ x1 + x2, data = data, mod = cmod_std2,
coordnames = c("x1", "x2"))
gear2b = update(gear1, cmod_std2)
gear2$call = NULL
gear2b$call = NULL
identical(gear2, gear2b)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.