predictラスター ライブラリの関数を使用して、glm一連のラスター スタック予測子でモデルを予測しています。私のモデルではいくつかの因子レベルが欠落していたので、predict関数を実行するためにそれらの層を NA に変更しました。現在、出力マップに欠損値があり、この問題を解決する方法を探しています。他の連続変数を使用してこれらの地域を予測する方法はありますか?
これらは私のサンプルコードです:
m1 <- glm(loc ~ factor(var1) + var2 + var3 , data= data)
var1 <- raster(raster1)
var2 <- raster(raster2)
var3 <- raster(raster3)
id <- which(!(var1[] %in% data$var1))
var1[id] <- NA
predictors <- stack (var1, var2, var3)
prd_data <- predict(predictors, m1, type= "response", na.action="na.pass")