我正在使用predict光栅库中的函数来预测我的glm模型在一组光栅堆栈预测器上。我的模型中缺少一些因子水平，因此我将这些层设置为 NA 以运行该predict函数。现在我的输出映射缺少值，我正在寻找一种方法来解决这个问题。有没有办法使用其他连续变量来预测这些区域？
这些是我的示例代码：
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")