Точечные паттерны в Spatstat

У меня возникли проблемы с созданием шаблона точек в пакете R Spatstat.

> HI06mfav <- read.table("MarHI06mfav.txt", header = TRUE)
>  attach(HI06mfav)
The following object(s) are masked _by_ '.GlobalEnv':

x, y
The following object(s) are masked from 'HI06mfav (position 3)':

x, y
> HI06mfav.P <- ppp(x, y, c(-1, 11), c(-1, 11))
Warning message:
In ppp(x, y, c(-1, 11), c(-1, 11)) :
10 points were rejected as lying outside the specified window

> HI06mfav
      x        y
1  4.100800 6.526913
2  3.710310 6.189953
3  3.030026 6.696919
4  2.747816 6.417674
5  2.955485 6.125787
6  3.299942 5.922913
7  2.932207 5.549772
8  2.681645 5.641663
9  1.858555 5.686039
10 2.069101 5.427891
11 2.387927 5.227459
12 2.530470 5.383951
13 2.783648 5.019703
14 2.219491 4.767310
15 6.591515 8.160451
16 6.600408 8.567767
17 3.098551 3.561141
18 2.634305 3.571341

Я использовал этот же код раньше, и все работало нормально, поэтому я не понимаю, почему он не работает сейчас. Я знаю, что мне может понадобиться альтернатива для attach(), но я не понимаю ошибку после того, как попытаюсь создать файл ppp.

Он даже делает то же самое с примером кода, взятого прямо из руководства Spatstat:

> w <- owin(c(-1,1), c(-1,1))
> w <- as.mask(w)
> X <- raster.x(w)
> Y <- raster.y(w)
> M <- (X^2 + Y^2 <= 1)
> pp <- ppp(x, y, c(-1,1), c(-1,1), mask=M)
Warning message:
In ppp(x, y, c(-1, 1), c(-1, 1), mask = M) :
10 points were rejected as lying outside the specified window
> pp
planar point pattern: 0 points 
window: binary image mask
100 x 100 pixel array (ny, nx)
enclosing rectangle: [-1, 1] x [-1, 1] units  

Я возвращаюсь к руководству по спатстату, чтобы попытаться понять это, но любая помощь будет очень признательна.


person user1549252    schedule 02.08.2013    source источник
comment
Ваш код spatstat работает нормально, поэтому он должен обращаться к столбцам прикрепленного фрейма данных. library(spatstat);x <- c(4.1008,3.71031,3.030026,2.747816,2.955485,3.299942,2.932207,2.681645,1.858555,2.069101,2.387927,2.53047,2.783648,2.219491,6.591515,6.600408,3.098551,2.634305);y <- c(6.526913,6.189953,6.696919,6.417674,6.125787,5.922913,5.549772,5.641663,5.686039,5.427891,5.227459,5.383951,5.019703,4.76731,8.160451,8.567767,3.561141,3.571341);ppp(x, y, c(-1, 11), c(-1, 11)).   -  person Andy W    schedule 03.08.2013
comment
Нет проблем, просто есть несколько точек, выходящих за пределы окна наблюдения, и все.   -  person Simon O'Hanlon    schedule 03.08.2013
comment
Что ж, @SimonO101, если вы посмотрите на результаты, в окно попадает 0 баллов! Сообщение об ошибке о маскировании столбцов в HI06mfav.P довольно простое, просто напечатайте x и y, чтобы увидеть, что они на самом деле, или используйте HI06mfav.P$x и т. д. при формировании шаблона точек.   -  person Andy W    schedule 03.08.2013
comment
Просто сделайте ppp(HI06mfav$x, HI06mfav$y, c(-1, 11), c(-1, 11))   -  person Simon O'Hanlon    schedule 03.08.2013