1

虽然这个问题已经在stackoverflow上讨论过，但我不太清楚如何解决我的具体情况

我有一个光栅文件，我将其读入：

prec<-raster("R:/rsrch/model/prec.tif")
summary(prec)
            prec
Min.         3.0
1st Qu.     52.0
Median     104.5
3rd Qu.    173.0
Max.      1850.0
NA's     19201.0

我的回答是在场和不在场（分别编码为 1 和 0）（共 9674 分）。我gam使用包构建了一个模型mgcv，如下所示：

mdl<-gam(species ~ s(prec), family=binomial)

当我运行它时，我收到以下错误：

 Error in model.frame.default(formula = species ~ 1 + prec,drop.unused.levels = TRUE): 
 invalid type (S4) for variable 'prec'

我认为这个错误是由我的栅格中的 NA 引起的。如果删除 NA，则网格的总数等于存在-不存在点的数量。我不确定现在该怎么做：

我尝试过的一件事是：

prec1<-values(prec)[!is.na(values(prec))]

这从我的光栅中删除了 NA，但是在我运行模型之后，它带有一些非常有趣的图表，我完全确定这是错误的。当预测变量是光栅并且响应是向量时如何运行游戏的任何建议。我在与栅格相同的投影中具有响应的 XY 坐标。谢谢您的帮助。

4

0 回答 0