По сути, я пытаюсь провести регрессионный анализ в рамках исследования распределения видов.
#Regression model
backg_trainM2 <- SpatialPoints(backg_trainM2,
proj4string=CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"))
trainM2 <- bind(pres_trainM2, backg_trainM2)
backg_trainM2 <- backg[groupM2 != 1, ]
pb_trainM2 <- c(rep(1, nrow(pres_trainM2)), rep(0, nrow(backg_trainM2)))
Работает отлично, пока я не дохожу до pb_trainM2 - тогда выдает ошибку:
Error in rep(1, nrow(pres_trainM2)) : invalid 'times' argument
Что может быть причиной этого? Я видел ошибку в других вопросах здесь, но решения, которые они использовали, не работали для меня.
Извиняюсь, если это вопрос новичка, я все еще очень плохо знаком с R.
EDIT:
Это Prestrain, это пространственный точечный объект S4 (с 2802 координатами):
pres_trainM2
class : SpatialPoints
Особенности: 2802
степень: -9,360213, 26,46285, 36,07, 61,2114 (xmin, xmax, ymin, ymax)
коорд. ссылка : + proj = longlat + datum = WGS84 + ellps = WGS84 + towgs84 = 0,0,0
В то время как Backg_trainM2 - это объект с двумя столбцами, длиной 800 строк. Вот первые 20 строк:
> backg_trainM2
x y
[1,] 34.62083333 61.26250
[2,] -5.31250000 50.10417
[3,] 19.72083333 49.52917
[4,] -4.27083333 58.23750
[5,] -2.66250000 38.17917
[6,] -7.64583333 54.01250
[7,] 15.18750000 56.64583
[8,] -5.37916667 39.80417
[9,] 35.10416667 56.70417
[10,] 29.77083333 68.18750
[11,] 11.97083333 32.27083
[12,] 9.12083333 35.15417
[13,] 17.97916667 46.40417
[14,] 0.62916667 49.19583
[15,] 6.12916667 36.40417
[16,] 2.77916667 50.11250
[17,] -4.87083333 33.52917
[18,] 25.64583333 61.60417
[19,] 12.04583333 45.35417
[20,] 34.26250000 40.97083
Большое спасибо за ваше время и терпение.