Hi Folks:
Another question for you...
I have run the following code to produce a directionally kriged map. Although the code seems to work, I am trying to understand the warning I am getting. Can you shed some light on this for me?
Here is the code:
```{r, warning = TRUE, fig.height = 8, fig.width = 5.5}
dbgrid4 <-
kriging(db_noout.db, dbgrid2, data.model.4dir, neigh.unique, radix = "dir")
plot(
dbgrid4,
pos.legend = 7,
cex = 0.7,
title = paste(property, "Kriging With Directionality"),
xlab = "X (UTM)"
)
plot(
dbgrid4,
name.contour = "dir.P_KH_md.estim",
nlevels = 11,
col = ("black") ,
add = TRUE
)
plot(
db_noout.db,
add = TRUE,
pch = 20,
name.post = property,
col = "blue",
size = 2
)
```
Here is warning:
is.na() applied to non-(list or vector) of type 'S4'is.na() applied to non-(list or vector) of type 'S4'is.na() applied to non-(list or vector) of type 'S4'is.na() applied to non-(list or vector) of type 'S4'NAs introduced by coercionNAs introduced by coercionis.na() applied to non-(list or vector) of type 'S4'is.na() applied to non-(list or vector) of type 'S4'is.na() applied to non-(list or vector) of type 'S4'is.na() applied to non-(list or vector) of type 'S4'NAs introduced by coercionNAs introduced by coercionNAs introduced by coercionNAs introduced by coercionNAs introduced by coercionNAs introduced by coercionNAs introduced by coercionNAs introduced by coercion
I sent Didier the db structure (there is one one outlier which is suppressed)
I also sent Didier the resulting figure which IS CORRECT! The question is to understand the warning.