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)