J'exécute une analyse MCP (domaine vital) en R sur mes données (appelé Ade) qui est un ensemble de points GPS pour les emplacements d'animaux.Le message d'erreur xy devrait être de classe SpatialPoints
Je courais le script suivant:
Ade <- read.csv("Ade.csv")
#Get the folder SP that is used to convert the data
library(sp)
#Convert
Ade <- matrix(runif(1564), ncol=2)
head (Ade)
Ade <- data.frame(Ade)
#Convert Ade to spatial points data frame by first ...
#Tell R that the coordinates are in the column X and column Y
coords <- SpatialPoints(Ade[, c("X", "Y")])
AdeSpat <- SpatialPointsDataFrame(coords, Ade)
#Specify the coordinate reference system of the data
proj4string(AdeSpat) <- CRS("+proj=XY +ellps=WGS84")
#View data
getClass(AdeSpat)
class(AdeSpat)
head (AdeSpat)
#Two steps in loading AdehabitatHR into R
install.packages("adehabitatHR")
library(adehabitatHR)
#Run MCP
mcp <- mcp(AdeSpat$coordinates, percent=95, unin = c ("km"), unout = c("km2"))
mais quand je lance ce que je reçois le message d'erreur suivant:
Error in mcp(AdeSpat$coordinates, percent = 95, unin = c("km"), unout =
c("km2")) : xy should be of class SpatialPoints
Quelqu'un peut-il me expliquer ce que cela signifie, et comment je peux résoudre le problème, j'ai cherché sur internet mais je ne trouve aucune réponse pour l'erreur "xy devrait être de classe SpatialPoints".
Quand je vérifie la trame de données, je reçois la sortie suivante: vous
> class(AdeSpat)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"
> head (AdeSpat)
coordinates ID ZONE X Y
1 (370984, 9752290) 1 36 M 370984 9752293
2 (370980, 9752290) 2 36 M 370980 9752288
3 (370980, 9752290) 3 36 M 370980 9752288
4 (370995, 9752290) 4 36 M 370995 9752287
5 (370995, 9752290) 5 36 M 370995 9752287
6 (371032, 9752310) 6 36 M 371032 9752312
Coordinate Reference System (CRS) arguments: +proj=XY +ellps=WGS84
Merci J