問題タブ [geostatistics]

R パッケージRandomFieldsを使用する場合、さまざまなコンポーネントを追加して、データの空間分布を記述するモデルを指定します。
例: データ内の近い点は、遠い点よりも互いによく似ているため、指数バリオグラムに従います。しかし、まだかなりのノイズがあり、バリオグラムにナゲット効果があります。あなたは次のように書くでしょう：

ナゲットは別途追加。
x座標の線形関数としてデータ値がさらに増加すると想像してください。+ RMtrend()トレンドを指定して、上記のモデルに追加することは間違いありません。
そして、これが私が立ち往生するところです。私は本当に-ページを取得しませんhelp(RMtrend)-特に例は私を立ち往生させます(この-関数は何をしR.pますか?ヘルプページは、-言語をまったく知らない私にとって完全に意味C不明です)。

とにかく、RMtrend引数は 1 つだけですmean。これは式ではなく、数値ベクトルまたは別のRMmodel. RMmodel線形トレンドまたは多項式トレンドを定義する必要がなかったかのように私には思えますRMtrend。RMtrendまた、式を与えるようにだまそうとしました

後でエラーが発生し、それRFformulaが誤用されたとのことです。RMtrend()また、シミュレーション グリッドのグリッド ノードの数の長さを持つベクトルを、それぞれの「平均」(= モデル化された傾向) を指定して与えることはできないようです。
たとえば、私のグリッドには 62250 ポイントがあるはずです。そして、トレンドなしでシミュレートすると、実際に 249 行 250 列の行列が得られますが、62250 の数値エントリを持つベクトルを にRMtrend() RFsimulate()渡しても、エラーが発生します (または R の致命的なエラーが発生します)。

RMtrendでは、どのようにして線形または多項式の傾向を与えることができるのでしょうか?

2012 年から、8187 の郵便番号 (および列としてリストされている約 10 の変数) に対して 8187 行のドイツのデータ フレームを取得しましたが、座標はありません。さらに、8203行の別のシェープファイルの座標を取得しました(ほとんど同じ郵便番号も含まれています)。

初期データ フレームの 8178 ケースに割り当てるには、8203 ケースの正しい座標が必要です。

問題: 必要な正しい割り当ての違いは、8178 ではなく、16 のケースが欠落している (8203 - 8187 = 16)、それ以上です。2012 年のいくつかの町 (郵便番号付き) は、最新のシェープファイルにはリストされていませんが、その逆も同様です。

(I) おそらく最も簡単な解決策は、2012 年の座標を取得することです (投影されていない: CRS("+init=epsg:4326"))。--> この目的のためのオープン ソース プラットフォームを知っている人はいますか? そして、郵便番号は正確に 8187 ですか?

(II) または: から別の年のデータ セットに座標を割り当てた経験がある人はいますか? - または、境界線と座標がわずかに変更されているため (特に、2012 年からデータをマップしてポリゴンで視覚化する必要がある場合)、新しいデータ セットの古い "and" にリストされていない町があるため、これを回避する必要がありますか?

この問題にアプローチする (そしてできれば解決する) 方法について、専門家のアドバイスをいただければ幸いです。

編集 - MWE:

したがって、

--> df2 にリストされておらず、次の場合:

--> df1 には記載されていません。

（I）と（II）に関するアイデアはありますか？