https://github.com/cran/gstat
Tip revision: 0d54c223814f3b1cd7515490af55d9e099a134f4 authored by Edzer J. Pebesma on 13 April 2007, 14:45:18 UTC
version 0.9-38
version 0.9-38
Tip revision: 0d54c22
na.action.R
library(gstat)
data(meuse)
data(meuse.grid)
set.seed(13131) # reproduce results
# select 10 random rows;
# create two missing values in the coordinates:
m = meuse.grid[sample(nrow(meuse.grid), 10), ]
m[c(2,8), "x"] = NA
## this is not allowed anymore:
try(krige(log(zinc)~1,~x+y,meuse,m, na.action = na.pass))
try(krige(log(zinc)~1,~x+y,meuse,m, na.action = na.omit))
try(krige(log(zinc)~1,~x+y,meuse,m, na.action = na.exclude))
try(krige(log(zinc)~1,~x+y,meuse,m, na.action = na.fail))
# select 10 random rows;
# create two missing values in the regressor variable:
m = meuse.grid[sample(nrow(meuse.grid), 10), ]
m[c(3,7), "dist"] = NA
krige(log(zinc)~dist,~x+y,meuse,m, na.action = na.pass)
krige(log(zinc)~dist,~x+y,meuse,m, na.action = na.omit)
krige(log(zinc)~dist,~x+y,meuse,m, na.action = na.exclude)
try(krige(log(zinc)~dist,~x+y,meuse,m, na.action = na.fail))