gbm.valuemap: Decision Support Tool that generates (Marine) Protected Area...

In SimonDedman/gbm.auto: Automated Boosted Regression Tree Modelling and Mapping Suite

Decision Support Tool that generates (Marine) Protected Area options using species predicted abundance maps

Description

Scales response variable data, maps a user-defined explanatory variable to be avoided, e.g. fishing effort, combines them into a map showing areas to preferentially close. Bpa, the precautionary biomass required to protect the spawning stock, is used to calculate MPA size. MPA is then grown to add subsequent species starting from the most conservationally at-risk species, resulting in one MPA map per species, and a multicolour MPA map of all. All maps list the percentage of the avoid-variables total that is overlapped by the MPA in the map legend.

Usage

gbm.valuemap(
  dbase,
  loncolno = 1,
  latcolno = 2,
  goodcols,
  badcols,
  conservecol = NULL,
  plotthis = c("good", "bad", "both", "close"),
  maploops = c("Combo", "Biomass", "Effort", "Conservation"),
  savedir = tempdir(),
  savethis = TRUE,
  HRMSY = 0.15,
  goodweight = NULL,
  badweight = NULL,
  m = 1,
  alerts = TRUE,
  BnW = TRUE,
  shape = NULL,
  pngtype = c("cairo-png", "quartz", "Xlib"),
  byxport = NULL,
  ...
)

Arguments

dbase	Data.frame to load. Expects Lon, Lat & data columns: predicted abundances, fishing effort etc. E.g.: Abundance_Preds_All.csv from gbm.auto.
loncolno	Column number in dbase which has longitudes.
latcolno	Column number in dbase which has latitudes.
goodcols	Which column numbers are abundances (where higher = better)? List them in order of highest conservation importance first e.g. c(3,1,2,4).
badcols	Which col no.s are 'negative' e.g. fishing (where higher = worse)?
conservecol	Conservation column, from gbm.cons.
plotthis	Vector of variable types to plot. Delete any,or all w/ NULL.
maploops	Vector of sort loops to run.
savedir	Save outputs to a temporary directory (default) else change to current directory e.g. "/home/me/folder". Do not use getwd() here.
savethis	Export all data as csv?
HRMSY	Maximum percent of each goodcols stock which can be removed yearly, as decimal (0.15 = 15 pct). Must protect remainder: 1-HRMSY. Single number or vector. If vector, same order as goodcols. Required.
goodweight	Single/vector weighting multiple(s) for goodcols array.
badweight	Ditto for badcols array.
m	Multiplication factor for Bpa units. 1000 to convert tonnes to kilos, 0.001 kilos to tonnes. Assumedly the same for all goodcols.
alerts	Play sounds to mark progress steps.
BnW	Also produce greyscale images for print publications.
shape	Set coastline shapefile, else uses British Isles. Generate your own with gbm.basemap.
pngtype	File-type for png files, alternatively try "quartz" on Mac. Choose one.
byxport	Dummy param for package testing for CRAN, ignore.
...	Optional terms for gbm.map.

Value

Species abundance, abundance vs avoid variable, and MPA maps per species and sort type, in b&w if set. CSVs of all maps if set.

Author(s)

Simon Dedman, [e-mail]

SimonDedman/gbm.auto documentation built on Feb. 9, 2024, 5:12 p.m.