IDW() или Krige() Ошибка: размеры не совпадают, когда отсутствуют значения

Question

Функции idw() и krige() из gstat пакета продолжают сообщения об ошибках, когда либо ответа или переменная предсказатель содержит недостающие значения (NA), даже когда na.action установлен в na.omit:

require(gstat) 
data(meuse) 
coordinates(meuse) = ~x+y 
data(meuse.grid) 
gridded(meuse.grid) = ~x+y 

meuse2 <- as.data.frame(meuse) 
meuse2[1, 'zinc'] <- NA 
meuse2 <- SpatialPointsDataFrame(SpatialPoints(meuse), meuse2) 

# idw response var 
int <- idw(zinc ~ 1, meuse2, meuse.grid, na.action = na.omit) 
# Error: dimensions do not match: locations 310 and data 154 

# krige response var 
m <- vgm(.59, "Sph", 874, .04) 
int <- krige(zinc ~ 1, meuse2, meuse.grid, model = m, na.action = na.omit) 
# Error: dimensions do not match: locations 310 and data 154 

# krige predictor var 
meuse3 <- as.data.frame(meuse) 
meuse3[1, 'dist'] <- NA 
meuse3 <- SpatialPointsDataFrame(SpatialPoints(meuse), meuse3) 
int <- krige(zinc ~ dist, meuse3, meuse.grid, model = m, na.action = na.omit) 
# Error: dimensions do not match: locations 310 and data 154 

Это ошибка? Нужно ли нам сначала фильтровать наши данные и объединять результаты с исходными кадрами данных? Разве нет более простого решения? Почему есть опция na.action?

Answer 1

Аргумент na.action имеет недостающие значения в пределах newdata (не locations или data).

Это четко указан в ?idw/?krige/?predict.gstat

function determining what should be done with missing values in 'newdata'. The default is to predict 'NA'. Missing values in coordinates and predictors are both dealt with.

Там нет никакого способа, чтобы иметь дело с NA значений в пределах locations или data (и, следовательно, ошибок, которые в основном говорят, что есть еще два значения в местах, как данные (то есть. х и у координаты точки недостающих данных)

Вы можете получить его на работу, удалив местоположение с отсутствующим значением

int <- idw(zinc ~ 1, meuse2[!is.na(meuse2$zinc),],newdata= meuse.grid)