2017-12-14 10 views
0

私はウィスコンシンの地図を取得しようとし、空間データポイントでClark-Evansテストを試みました。spatstatのR:点が外にあると拒否されました

Warning message: 244 points were rejected as lying outside the specified window

私は次のコードを使用:しかし、私はウィスコンシン州のマップをインポートし、窓に変換しようとすると、私はエラー会っ

wisc.map <- map("state", "wisconsin", fill=T) 
w <- owin(poly=data.frame(x=rev(wisc.map$x), y=rev(wisc.map$y))) 

を私は走った後、私はエラーを得ました spat = as.ppp(data, w)

それはプロットは以下のように示したすべての私のデータポイントはウィスコンシン州にあることを明確に示します。 enter image description here

は、私はここにいくつかのサンプルデータを置く:

latitude longitude 
1 43.30716 -89.43137 
2 43.30388 -89.42823 
3 43.31845 -89.44372 
4 43.33766 -89.44634 
5 43.50815 -89.52010 
6 43.34179 -89.01988 

何が起こっているため任意のアイデア?


P.S

おかげで私は私のデータのうち、クラスタリングパターンが存在するかどうかを確認したいのでクラーク・エヴァンスのテストをしたいです。私は、データから得られたホットスポットが統計的に有意であることを確認したい。私は多くの文献をチェックしましたが、それを直接テストする方法は見つからなかったので、Clark Evans Testを使用しなければなりません...クラスタリングパターンやホットスポットの重要性をチェックするためのテストは誰でも知っていますか?

答えて

1

これはエラーではありませんが、出力として明示的に警告しています。それはあなたがpppに供給したポイントのうち244が、観測ウィンドウであると主張している多角形領域の外にあることを示しています。それらは端に非常に近く、ウィスコンシンに近いポリゴンは詳細ではないかもしれません。

注:緯度と経度のデータにspatstat関数を使用する意味がないのでspatstatパッケージは現在、球面幾何学を処理しません。プロジェクションは、spパッケージのspTransformで行われていましたが、より現代的なシンプルな機能パッケージsfを使用する方が便利です(st_transform)。

サイドノート:なぜこのデータについてClark Evansテストを行いますか?最も近い隣接距離が、この点パターンと、同一の観測ウィンドウ上の同じ期待数点とを有する均一なポアソン過程と類似しているかをチェックする。したがって、参照プロセスでは、ウィスコンシン全域に均等にデータを分配し、データにはウィスコンシンの南西部にのみポイントがあります。

+0

実際、私のすべてのデータはウィスコンシンの南西部にあります。私は自分のデータの中にクラスタリングパターンが存在するかどうかをチェックしたいので、Clark Evans Testをしたい。私は、データから得られたホットスポットが統計的に有意であることを確認したい。私は多くの文献をチェックしましたが、直接テストする方法は見つからなかったので、Clark Evans Testを使用しなければなりません...それに対処するためのより良いテストがありますか? – Leonard

+0

Clark-Evansテストでは、観測ウィンドウ全体に均一なパターンが観察されていることを前提としています。分析を進める前に、実際の観測ウィンドウが何であるかを知る必要があります。 –

関連する問題