Erreur ENFA package adehabitatHS

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Margaux Julien
Messages : 2
Enregistré le : 18 Mai 2018, 07:39

Erreur ENFA package adehabitatHS

Messagepar Margaux Julien » 07 Juin 2018, 10:24

Bonjour,

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 !

Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité

cron