問題タブ [spdep]

spdepR のライブラリを使用して Queen Contiguity Neighbors をプロットした次の結果に困惑しています。

ワールド形状ファイル wrld_simp を使用して重心点 (空間座標) を取得する

女王の隣人を取得するために、私はそうします

次に、形状ファイルをプロットします

およびクイーン隣接ネイバー:

それをプロットすると (または単に を実行するとsummary(queen_nb))、多くの近隣が接続されていないことがわかります (たとえば、チリは孤立しています)。

ヒントはありますか？

ブラジルの大西洋岸森林の森林伐採を調査するプロジェクトのために、空間自動コレログラムを実行しようとしています。

しかし、なぜこの問題に直面しているのかについては混乱しています。

問題

コードの最初の部分を実行すると、次のエラーが表示されます

エラー: ncol(x) == 2 は TRUE ではありません

私のコードは

そして、このコードを実行したいと思います

私のデータは

見出し付きのベクター データ セット

k-means クラスタリングを使用して、自己組織化マップからコードブックをクラスタリングしたいと考えています。ただし、データの「空間的」性質を考えると、隣接するノードのみが一緒にクラスター化されるように、クラスター化を制限したいと考えています。skaterいろいろ調べた結果、spdepパッケージの機能を使ってみることにしました。

これが私が行ってきたことの例です。

上記のコードが本当にぎこちなく、洗練されていないことは承知しています (ご容赦ください)。使用してみmstree(nb2listw(nbs))[,1:2]ましたが、すべてのリンクがリストされていませんでした。これが何をしているのかよくわからないので、エッジのマトリックスを手動で作成しました。

次に、この行列をスケーター関数に使用しようとしました

しかし、次のエラーメッセージが表示されます。 Error in colMeans(data[id, , drop = FALSE]) : error in evaluating the argument 'x' in selecting a method for function 'colMeans': Error in data[id, , drop = FALSE] : subscript out of bounds

ただし、mstreeエッジを使用すると、エラー メッセージは表示されませんが、結果はまったく意味がありません。

このエラー メッセージに関するヘルプ (または、私がやりたい空間的に制約されたクラスタリングの実行方法に関する別の提案) は大歓迎です。

この質問で受け入れられた回答の議論で、user3050574 は次のように述べています。

「... mat2listw は、現在バイナリ形式の行列から行の標準化された重み行列を作成することを理解しています。」

これは私がこの種のことわざについて読んだ唯一の場所です. 「spdep」Rドキュメントでは、

「この関数は、正方形の空間ウェイト マトリックス、オプションでスパース マトリックスをウェイト リスト オブジェクトに変換します...」

この変換には行の標準化が含まれますか?

適用したい正確な重みとして各要素を持つ重み行列があります。したがって、関数 mat2listw が行の標準化された重み行列を生成するかどうかを確認することは、私にとって非常に重要です。

Rで使用spdepし、Jarque-Bera検定を実行して、空間誤差モデルの最尤推定から得られる残差の非正規性を確認した後、一般化モーメント法 ( ) を使用しました。ここで、このステップにとどまるか、 SHACモデルの推定に進むかを決定するために、現在の不均一分散の問題からの残差をテストしたいと思います。GMerrorsarGMerrosar

モデルspdepから異分散性をチェックできるのはどのテストですか? GMerrorsarこの場合、Breusch-Pagan ( bptest.sarlm) 検定は適用できないことがわかりました。