Hi all, I got these set of codes to compute mcp which worked well in drawing the polygons but I don't have any idea about how the get the areas of the various polygons calculated for each study site. Here is what I did using version 2.15.2; data for location 1 attached with mail.
library(ggplot2) library(plyr) ##Location 1 mcpyankari<-read.csv("mcpyankari.csv") head(mcpyankari) find_hully <- function(mcpyankari) mcpyankari[chull(mcpyankari$long, mcpyankari$lat), ] hullsy <- ddply(mcpyankari, "mcpyankari$Individual", find_hully) plot <- ggplot(data = hullsy, aes(x = long, y = lat, colour=Individual, fill = Individual)) + geom_point() + geom_polygon(data = hullsy, alpha= 0.5) + labs(x = "Longitude", y = "Latitude") plot ##Location 2 mcpkainji<-read.csv("mcpkainji.csv") head(mcpkainji) find_hullk <- function(mcpkainji) mcpkainji[chull(mcpkainji$long, mcpkainji$lat), ] hullsk <- ddply(mcpkainji, "mcpkainji$Individual", find_hullk) plot <- ggplot(data = hullsk, aes(x = long, y = lat, colour=Individual, fill = Individual)) + geom_point() + geom_polygon(data = hullsk, alpha = 0.5) + labs(x = "Longitude", y = "Latitude") plot Thanks Chima WebRep Overall rating This site has no rating (not enough votes) WebRep Overall rating This site has no rating (not enough votes) WebRep Overall rating This site has no rating (not enough votes) WebRep Overall rating This site has no rating (not enough votes) WebRep Overall rating This site has no rating (not enough votes) WebRep Overall rating This site has no rating (not enough votes) WebRep Overall rating This site has no rating (not enough votes)
mcpyankari.csv
Description: MS-Excel spreadsheet
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo