R (spatstat) : 내장 함수가 원하는 한계로 평가되지 않음

Question

그래서 내가 만난 문제는 아주 기본적인 것이지만 해결할 수없는 것 같습니다. 나는 spatstat의 F, G, J, K 및 L 함수를 사용하여 작업했으며, 다른 점 프로세스의 F 함수와 동일한 플롯에서 한 포인트 프로세스의 F 함수를 플롯하려고합니다. 이것은 내 포인트 프로세스 중 하나가 함수를 r = 20으로 평가하고 다른 하나는 r = 15 ('r'은 x 축의 독립 변수 임)를 제외하고는 쉽습니다. 둘 다 같은 r 값, 바람직하게는 20으로 평가되기를 바랍니다. 어떻게해야합니까? 본질적으로 이것은 문제입니다. 어떤 이유로,

mydata.Fest <- Fest(mydata) 

다음 코드 행 :

plot(mydata.Fest) 

이

이 같은 'R'로 평가하지 않습니다 (

먼저 나는 함수의 값을 포함하는 데이터 프레임을 생성 이 코드 줄 독립 변수) 값 :

plot(mydata.Fest, xlim=c(0,20), ylim=c(0,0.8)) 

두 플롯

이 (xlim =의 C에서 0,20 이동) 및 ylim = c (0,0.8)이지만, 어떤 이유에서 두 번째 것은 x = 15에서 함수를 평가하는 것을 멈추는 것입니다. 두 변수가 같은 객체로부터 구획 된 것을 고려할 때, 왜 같은 것을 생성하지 않는지에 대한 단서가 없습니다 음모.

나는 F, J, K 함수로 작업 할 때 한계 (특히 독립 변수)를 조정하면 R이 예기치 않게 중단되는 것을 발견했습니다. 이 문제가 있거나 해결책을 알고 있다면 알려주십시오.

감사합니다.

Answer 1

FAQ입니다. spatstat.org FAQ 페이지에서 논의됩니다.

Fest 및 Kest와 같은 spatstat의 요약 함수에는 독립 변수에 '권장 범위'가 있으며, 이는 함수가 계산 된 값의 범위보다 짧을 수 있습니다.권장 범위는 기능 평가가 통계적으로 신뢰할 수있는 범위입니다. 기본적으로이 함수는 권장되는 한도로만 표시됩니다. (이것은 공간 통계의 표준 사례입니다. 우리가 이것을하지 않으면, 함수 값이 폭발하기 때문에 대부분의 플롯이 이상하게 보일 것입니다. 신뢰할 수있는 정보는 왼쪽 하단의 작은 부분을 차지하도록 축소됩니다. 플롯.)

r의 유효 범위와 권장 범위를 검사하려면 함수 객체 (이름을 입력하기 만하면 됨)를 인쇄하고 마지막 몇 줄의 출력을 봅니다.

플롯 제한을 제어하려면 xlim 및 ylim을 사용하십시오. 패키지 저자 입력, 앤디에 대한

Answer 2

음 조금 이상합니다. 미안하지만 대답은 아니지만 다른 사람들이 설명했던 행동의 재현 가능한 예를 남겨 두어 조사하는 것이 가장 좋을 것이라고 생각했습니다.

내가 과거에 한 일은 일관된 포인트 패턴을 평가할 수있는 일관된 벡터 벡터를 만드는 것입니다.

library(spatstat) 
data(lansing) 

mydata1 <- lansing[lansing$marks == "blackoak",] 
mydata2 <- lansing[lansing$marks == "hickory",] 

my_r <- (0:40*.002) 

mydata1.Fest2 <- Fest(mydata1, r = my_r) 
mydata2.Fest2 <- Fest(mydata2, r = my_r) 

plot(mydata1.Fest2) 
plot(mydata2.Fest2, add = TRUE) 

mydata1.Fest2$r 
mydata2.Fest2$r 

볼 수 있듯이 거리 저장소는 동일하지만 기능은 동일한 세트에서 평가되지 않았습니다. 나는 이것이 아마도 Fest에 고유 한 동작이며, Kest을 사용할 때 예상되는 해결책을 얻을 수 있었는지 확인했습니다. 아래 예.

mydata1.Kest <- Kest(mydata1) 
my_k <- mydata1.Kest$r 
mydata2.Kest <- Kest(mydata2, r = my_k) 

plot(mydata1.Kest) 
plot(mydata2.Kest, add = TRUE) 

나는 Fest이 예상되지 않는 행동을 가정,하지만 난 기술에 익숙한 아니에요 그것은 리플리의 K와 다른, 그래서 확실히 오인 될 수있는 방법이 몇 가지 이유로 예상되는 동작 나에게 알려지지 않았다. 출력 장치에서 축 제한을 설정할 때 설명했던 것과 비슷한 예기치 않은 충돌이 발생했지만 재현 할 수 없었습니다. 어쩌면 당신은 무언가에 있거나 우리 둘다 미쳤습니다!

Answer 3

문제에 관해 Adrian Baddeley (spatstat '의 작성자)에게 연락하여이 예제가 add = TRUE 인 경우 plot.fv의 버그를 보여줍니다. 이 문제는 spatstat 버전 1.28-1에서 수정 될 예정입니다.