特定のラグ（R単位）で空間自己相関を使用する問題

Question

数日から、私はRの空間自己相関を含む新しい困難な空間分析に苦労しています。具体的には、aで設定されたポイント間の自己相関を検証することに興味があります。グリッドは50m（多かれ少なかれ）です。私の目的は、これらの点（私がデータを収集した場所）間の自己相関をテストし、自己相関がそれらの間の距離を増やすことが減少するかどうかを検証することです（これは予想されます）。私の考えは、各点（50m、100m、150mなど）の周りに特定のメーターの異なる半径を生成し、モランのI自己相関指数をテストすることです。最後に、ggplotを使用して各特定の距離結果でMIを表示したいと思います（MI出力を取得すると簡単に取得できます）。

私の開始データフレームには、収集されたデータのポイントID、その特定ポイントで測定された値（z）、経度（x）と緯度（y）のcoloumn、および緯度（y）のデータが表示されます次のように：

#install libraries 
library(sp) 
library(spdep) 
library(splm) 
library(ape) 

ID<- c(1,2,3,4,5,6) 
x<-c(20.99984,20.99889, 20.99806,20.99800,20.99700,20.99732) 
y<-c(52.21511,52.21489,52.21464,52.21410,52.21327,52.21278) 
z<-c(1.16,0.54,0.89,0.60,1.27,1.45) 

data <- data.frame(ID,x,y,z) 

私は実際に私が興味を持ってるものを示しており、オンラインで多くのことを読んで、このチュートリアル https://mgimond.github.io/Spatial/spatial-autocorrelation-in-r.html#morans-i-as-a-function-of-a-distance-band が見つかりました：しかし、それは本当に私の座標から始まる、実際の最初から動作していません、私は問題があると思うし、Rの適切な形式でそれらを変換する方法を知らない。これは私が得るエラーメッセージです：

data <- data.frame(dataPOL$Long , dataPOL$Lat, dataPOL$Human_presence) 
coordinates(data) <- c('x','y')` 
proj4string(data) <- "+init=epsg:4326" 
S.dist <- dnearneigh(coordinates, 0, 50) #radius of 50 meters 

Error in dnearneigh(coordinates, 0, 50) : Data non-numeric 

Answer 1

私はどんな答えを受信しませんでしたが、私は解決策を見つけてしまった：私は（私の場合、モランIで）最も使用されるパッケージはRで空間的自己相関で動作することを発見した はspdepと猿です。

私は両方を試みました：sp​​depはまだ動作しませんでしたが、猿はしました。ここで私は私の特定のケースのために続くチュートリアルは次のとおりです。

https://stats.idre.ucla.edu/r/faq/how-can-i-calculate-morans-i-in-r/

前モラン指数を計算、あなたは距離行列を生成する必要があり、私はパッケージ「フィールド」から「rdist.earth」でそれをやりました。 この関数は、座標に基づいてデータポイントの各セット間の距離を測定します。この関数は、世界がフラットではないことを認識し、大きな円の距離として知られているものを計算します。私は特定のケースで距離をKmで指定しました。 （popdists> 0 & popdists < = 0.06）60メートルの＃半径（ を覚えている -

ライブラリ（APE）

pop.dists.1 <：

はモランIを計算するために、私はこれを実行しましたそのフィールドパッケージはキロまたはマイルで動作します） Moran.I（mydataframe $ myzvariable、pop.dists.1）

これは、私は、この特定の半径で得た出力です：

pop.dists.1 < - （popdists> 0 & popdists < = 0.06）＃1 60メートル Moran.I（dataPOL $ Human_presence、pop.dists.1） $観察[1] 0.3841241 #Moran度：間-1から1、ここで、60 m以内ポイント は$は 予想

を自己相関している[1] -0。009615385

$は、SD [1] 0.08767598

の$ p.value [1] 7.094019e-06

私は、私が興味を持っていた距離の公式を繰り返す：それは本当によくして上げる働きを距離、Moran Iのインデックスはおおよそ0です（これは私が予想したものです）。

私はz変数の空間自己相関の傾向に従うために、常にggplotを使って単一の出力をプロットします。

これは必要な場合に役立ちます。