自動化されたバリオグラム モデリングに automap パッケージを 2 日前に使い始めて、いじってみました (手動で行うのは面倒で、あまり多くの利点はありません... ただし、automap が Matérn 関数を本当に気に入っているように見えるのは奇妙だと思います)バリオグラムのモデリング用)。私の頭に浮かんだのは、データセットに負の値が存在することでしたが、それらはこのパッケージによって生成されたマップには表示されませんでした。
これは私が貼り合わせたばかりの図です。1 は、autoFitVariogram と autoKrige を使用した後の出力です。負の値はないようです。グラフィック出力が気に入らなかったので、autoFitVariogram (automap から) (Nr. 2) からのバリオグラムを使用する gstat で通常のクリギングを実行し、負の降水量の領域を明確に示します。autoKrige から生成されたマップを QGIS にエクスポートして再確認したところ、多くの負の値 (Nr. 3; 黒く塗りつぶされた領域) がありました。
だから私の質問：
1.) autoKrige によって生成されるクリギング予測マップは、「if var.pred < 0 then 0」などのスキームに従うマップを生成するようです。これはそのようなマップを作成するのに役立つことがわかりますが、これは実際のデータでどのように意味があるのでしょうか? データセットをさらに分析するために使用したいのですが、負の降雨値が少し表示されます... 非現実的です。この「if var.pred < 0 then 0」スキームをデータにも適用することは、物理的に妥当である場合、数学的に適切でしょうか?
2.) automap パッケージの独自の引数として一般化最小二乗法 (GLS) を使用して、さらにバリオグラム パラメーターを決定できます。ただし、gstat ( fit.variogram.reml )からの REML フィットを使用したいと考えています。自分の REML 推定 (修正したバリオグラムも生成) に置き換えてみvariogram$var_modelましたが、autoKrige でクリギングを実行すると、元のバリオグラム モデルが再び使用されます。パッケージで提供されている以外の方法でモデルを適合させる方法はありますか? (autoKrige マニュアルには引数 'data_variogram' がありますが、これは別のことをしているようです)。
ありがとう！