Je cherche à réaliser une ENFA avec le package adehabitatHS. Pour cela, j'ai suivi l'exemple donné par l'aide de R, mais j'obtiens une erreur au moment de faire l'enfa.
Voilà mon code :
Code : Tout sélectionner
locs <- read.table("Stephanitis nashi ENFA.csv", header = TRUE, sep = ";")
bio1 <- read.asciigrid("bio1.asc", as.image = FALSE)
bio2 <- read.asciigrid("bio2.asc", as.image = FALSE)
xy.sp = SpatialPoints(locs)
join <- cbind(bio1,bio2)
xy.sp <- as(xy.sp, "SpatialPointsDataFrame")
join2 <- as(join, "SpatialPixelsDataFrame")
bio1b <- as(bio1, "SpatialPixelsDataFrame")
tab <- slot(bio1b, "data")
tab <- na.omit(tab)
pr1 <- count.points(xy.sp, bio1b)
pr <- slot(pr1, "data")[,1]
pc <- dudi.pca(tab, scannf = FALSE, nf = 1)
enfatest <- enfa(pc, pr, scannf = FALSE)
J'ai un warning lorsque je réalise le Count.points() qui me donne :
Code : Tout sélectionner
Warning messages:
1: In count.points(xy.sp, join2) :
several columns in the SpatialPointsDataFrame, no id considered
2: In points2grid(points, tolerance, round) :
grid has empty column/rows in dimension 2
Puis l'erreur finale :
Code : Tout sélectionner
Error in eigen(Se) : infinite or missing values in 'x'
Mon tableau locs est constitué des longitudes et latitudes de chaque individu dont voilà un extrait
decimallatitude decimallongitude
35.66006 139.68521
35.66006 139.68521
35.66006 139.68521
35.66006 139.68521
35.95455 139.0381
35.95455 139.0381
35.86161 139.64488
35.86161 139.64488
35.86161 139.64488
35.86161 139.64488
35.86161 139.64488
35.86161 139.64488
35.86161 139.64488
35.86161 139.64488
Si vous pouviez m'aider à avancer je vous en remercie d'avance !