Аргумент функции среднего квадрата ошибки (MSE) "poly"
Я пытаюсь оценить начальное значение пропускной способности для ядра, сглаживающего мои данные, вычисляя среднеквадратичные ошибки (функция "mse2d"), но застрял с ошибкой в коде реализации многоугольника в аргументе функции.
Я следую примеру в Bivand, Pebesma и Gómez-Rubio (2013, стр. 186) с набором данных spatstat Redwood, который работает отлично:
library(spatstat)
library(splancs)
data(redwood)
spred<-as(redwood, "SpatialPoints")
#this is the function that computes the Mean Squared Error of the points
#distribution: mse2d(pts,poly,nsmse, range)
mserwq <- mse2d(as.points(coordinates(spred)), as.points(list(x = c(0,1, 1, 0), y = c(0, 0, 1, 1))), 100, 0.15)
bwq <- mserwq$h[which.min(mserwq$mse)]
bwq
Хорошо. Я воспроизвел этот код, используя мой набор данных. Вот мой код с подвыборкой точек, которая возвращает ту же ошибку:
m23.Xs<-c(17349,13212,11551,16659,9461,12062,12802,9638,9835,9803)
m23.Ys<-c(576,13600,6372,11763,11081,5462,15802,11667,11552,11121)
# matrix of coordinates
loc23<-matrix(c(m23.Xs,m23.Ys),nrow=10,byrow=FALSE)
MSEm23<-mse2d(as.points(coodinates(loc.m23),as.points(list(x=c(0,20000,20000,0),y=c(0,0,20000,20000))),100,0.15))
Error in storage.mode(poly) <- "double" : argument "poly" missing, without pattern
На самом деле, я перевел сообщение об ошибке на английский, не знаю точно ли это.
Как видите, мой код воспроизвел структуру в книге, но я продолжаю получать ошибки с этим аргументом "поли". Я пытался научиться конструировать полигоны в пакете splancs, но не смог понять, как и не нашел документации по этому поводу. После поиска в сети я также попытался поместить объекты и матрицы "owin" вместо списка в аргумент "poly", но ни одна из них не сработала.
Буду признателен за любые комментарии или помощь. Спасибо всем заранее.
Лейла