この問題は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 を削除すると、グリッドの総数は存在/不在ポイントの数に等しくなります。今何をすべきかわからない：
私が試したことの1つはこれでした：
prec1<-values(prec)[!is.na(values(prec))]
これにより、ラスターから NA が削除されましたが、モデルを実行した後、いくつかの非常に面白いグラフが表示されますが、これは完全に間違っていると確信しています。予測子がラスターで応答がベクターの場合にゲームを実行する方法に関するアドバイス。ラスターと同じ投影法で応答の XY 座標があります。ご協力ありがとうございました。