2017-12-14 21 views
0

위스콘신지도를 얻고 공간 데이터 포인트에 대한 Clark-Evans 테스트를 시도했습니다.R : spatstat에서 오류가 있습니다. 바깥쪽에 누워있는 것으로 거부되었습니다.

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으로 완료해야했지만 이제 st_transform 기능을 가진 최신의 간단한 기능 패키지 sf을 사용하는 것이 더 편리 할 수 ​​있습니다.

사이드 노트 : 왜이 데이터에 대해 클라크 에반스 테스트를 수행 하시겠습니까? 가장 가까운 이웃 거리는 동일한 관측점 창에있는 동일한 포인트 수와 을 가진 균일 한 포아송 과정과이 점 패턴에 대해 유사하게 보이는지 확인합니다. 따라서 참조 프로세스는 모든 위스콘신에 대해 데이터와 동일한 수의 포인트를 균등하게 분배하며 데이터에는 위스콘신의 남서부에만 포인트가 있습니다.

+0

사실 내 모든 데이터는 위스콘신의 남서부에 있습니다. 내 데이터 중에서 클러스터링 패턴이 존재하는지 확인하고 싶기 때문에 Clark Evans Test를하고 싶습니다. 데이터에서 얻은 핫스팟이 통계적으로 중요한지 확인하려고합니다. 나는 많은 문학을 확인했지만 그것을 테스트 할 수있는 직접적인 방법을 찾지 못했습니다. 그래서 Clark Evans Test를 사용해야합니다 ... 더 나은 테스트를 알고 있습니까? – Leonard

+0

Clark-Evans 테스트는 전체 관측 창에서 균질 패턴을 관찰 한 것으로 가정합니다. 이는 분명히 여기에 해당하지 않습니다. 분석을 진행하기 전에 실제 관찰 창이 무엇인지 알아야합니다. –