空間自己相関を検出するための座標投影に関する問題R

Question

緯度と経度の座標を含むデータセットと、それぞれが別々の列にあり、数値として格納されている属性情報があります。これらの座標は、地理座標系WGS 1984に基づいてジオコーディングされています。 "sp"パッケージを使用してバブルプロットで視覚化したいと考えているデータでは、大きな空間自己相関があることがわかります。ここでは、https://beckmw.wordpress.com/2013/01/07/breaking-the-rules-with-spatial-correlation/のように、オンラインで他の人の例をモデリングしています。私たちは、座標が "SP" の中にコマンドを使用しようとすると、しかし、我々は、エラーメッセージが出続ける：

コード例： COORDS < - data.frame（鉛の$ X、リード$ Y） 座標（COORDS） < - TRUE/FALSEが

を必要as.character（NA）欠損値 - C（ "lead6.X"、 "lead6.Y"）であれば（NCHAR（projargs）== 0）で

エラーが<をprojargs

ダイレクトコードは機密情報であり、インターネットにアクセスできない仮想環境でホストされているため、直接コードを読み込むことはできません。なぜ誰かがこれが起こっているかもしれないアイデアを持っていますか？私たちはproj4パッケージを調べましたが、投影システムを指定する方法を理解することはできません。誰かがRのパッケージや空間自己相関を視覚化する方法を知っていれば、それも非常に高く評価されます。

Answer 1

あなたのコードはちょっと "奇妙な"ものです。座標だけを含むデータセットを構築しようとしているようです。リード$ Zは、一般的な「変数」あなたは "あなたに適切なを与えるべき、

coordinates(data) <- c('X','Y')` 
proj4string(data) <- "+init=epsg:4326" 

、その後、検査するに対応して、

data <- data.frame(lead$X, lead$Y, lead$Z) 

：AFAIUは、あなたがこの行で何かを必要とするかもしれませんSpatialPointsDataframe "とWGS84の緯度経度座標（最初の行を削除することもできます）を使用して、すべての変数をspatialpointsデータフレームのデータに保持します。

HTH