The following contains the R code used to run the analysis for kittiwake.Gannets
setwd("~/Work/Project/CEH/Displacement/SimplifiedModel/CurrentModel/")
source("Simplified-Displacement-Model-Version17.R")
pars.full <- (read.csv("Parameters/parameters-gps.csv"))[,3:9]
meta <- list(wf = read.csv("Meta/wfinfo.csv"), colony = read.csv("Meta/colonies.csv"), colonysize = read.csv("Meta/colonysize.csv"), wfstats = read.csv("Meta/wfstats.csv"))
levels(meta$colony$colony) <- c("ForthIslands", "Boddam", "FAR", "Fowlsheugh", "ForthIslands", "StAbbsHead", "TPL")
levels(meta$wfstats$SPA)[3] <- "ForthIslands"
Species-specific part :
dat <- read.csv("/tmpshare/adam/projects/seabirds/Displacement/out-sixth/kittiwake.csv")
levels(dat$Site)[3] <- "ForthIslands"
out.kittiwake <- vsdm.sensitivity(dat, species = "Kittiwake", meta = meta, pars.full = pars.full, ppop = 1/3)
write.csv(getsummary.zones(out.kittiwake), file="Output/sensitivity-kittiwake.csv", quote=FALSE, row.names=FALSE)