Suite

Ajouter des « marques » au fichier de points dans spatstat ?


En utilisant l'aide de ce site, j'ai déjà obtenu le code suivant à partir duquel je peux effectuer certaines analyses. Cependant, j'aimerais maintenant effectuer des analyses sur une certaine variable (colonne z) dans le fichier de formes de points. Pour ce faire, je dois établir des « marques » sur la colonne z (la 7e colonne). Cette colonne fournit soit un « 4 » soit un « 10 », définissant deux types de point.

Quelqu'un peut-il aider?

##Exécutez les bibliothèques : library ("spatstat") library ("maptools") library ("sp") library ("rgdal") ## définissez le répertoire de travail setwd ("Z:/CP/spatstat/data/") # #lire dans les fichiers de formes Points <- readOGR(getwd(), "MyPointShapefile") StudyArea <- readOGR(getwd(), "MyPolygonShapefile") ##créer la fenêtre à partir du polygone shp w <- as.owin(StudyArea) ## forcer les points vers la fenêtre pts.ppp <- as.ppp(coordinates(Points), w) plot(pts.ppp)

Si les marques sont créées à partir d'un autre attribut, vous pouvez probablement simplement utiliser une instruction if-else. Vous pouvez simplement créer une nouvelle colonne en créant une variable, c'est-à-dire.

Points$NewMarks <- ifelse(X, Y, Z)

Où:

X = votre argument (instruction if), Y = votre instruction then (pour vous : 4) et Z = instruction Else (pour vous : 10)