基本上,我正在尝试进行回归分析,作为物种分布研究的一部分。
#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还是很陌生。
编辑:
这是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是2列对象对象,长度为800行。这是第1 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
非常感谢您到目前为止的所有时间和耐心。