plots: Abundance and Sector Effects Plots

In ABbiodiversity/cure4insect: Custom Reporting for Intactness and Sector Effects

Abundance and Sector Effects Plots

Description

Plots the sector effects for a single or a group of species.

Usage

plot_sector(x, ...)
## S3 method for class 'c4iraw'
plot_sector(x,
    type=c("unit", "regional", "underhf"), main, ylab, subset=NULL, ...)
## S3 method for class 'c4idf'
plot_sector(x,
    type=c("unit", "regional", "underhf"), main, ylab, subset=NULL, ...)

plot_intactness(x, ...)
## S3 method for class 'c4idf'
plot_intactness(x, type=c("SI", "SI2"), col, ...)

plot_abundance(species, type, plot=TRUE, paspen=0, ...)

Arguments

x	an object of class 'c4iraw' (from calculate_results) or a data frame (class 'c4idf') with flattened (1-row) results per species.
species	character, species ID (see get_all_species).
type	type of the plot, see Details.
main	title for the plot, if single species results are displayed the default is to use the species ID.
ylab	character, optional label for the y axis.
subset	subset of sectors to be plotted, can be any suitable index.
plot	logical, if a plot is to be drawn.
paspen	numeric in [0, 1], the probability of aspen occurrence (proxy for climatic suitability for treed vegetation).
col	color.
...	other possible arguments passed to underlying functions, e.g. ylim, xlab, ylab, col (for sector colors), or method (one of "kde", "fft", or "hist").

Details

"unit" type sector effects are based on regional current and reference abundances, and the regional sector effects are standardized by footprint area.

"regional" sector effects includes native and disturbed habitats when comparing regional abundance under current and reference conditions.

"underhf" (under human footprint) type sector effects consider only the abundance that us 'under the footprint', meaning that the current designation is disturbed.

The single species sector effect plots are different kinds of bar plots. The multi-species plot represents violin (carrot, vase) plots based on kernel density, fast Fourier transform, or binning (histogram).

Intactness plots are either one sided (0-100%, "SI"), or two-sided (0-200%, "SI2") differentiating increased (>100%) and decreaser (<100%) species.

Abundance plots depend on the type argument: "veg_coef" type abundance plots show relative abundances across various land cover (incl. disturbance) classes, "soil_coef" shows relative abundance by soil and disturbance types. The "veg_lin" and "soil_lin" types show average relative abundances compared to 10% vegetated (soft) and non-vegetated (hard) linear disturbance.

Value

Called for the side effect of drawing a plot. Returns the plotted data invisibly.

Author(s)

Peter Solymos <[e-mail]>

See Also

calculate_results

Examples

## Not run: 
## *res*ults from calculate_results, all province, all species
fn <- paste0("http://science.abmi.ca/reports/",
    getOption("cure4insect")$version, "/misc/raw_all.rda")
con <- url(fn)
load(con)
close(con)

plot_sector(res[["CanadaWarbler"]], "unit")
plot_sector(res[["CanadaWarbler"]], "regional")
plot_sector(res[["CanadaWarbler"]], "underhf")

z <- do.call(rbind, lapply(res, flatten))
class(z) <- c("c4idf", class(z))
plot_sector(z, "unit") # all species
plot_sector(z[1:100,], "regional") # use a subset
plot_sector(z, "underhf", method="hist") # binned version

plot_intactness(z, "SI")
plot_intactness(z, "SI2", method="hist")

## land cover associations
load_common_data()
plot_abundance("Achillea.millefolium", "veg_coef")
plot_abundance("Achillea.millefolium", "soil_coef")
plot_abundance("Achillea.millefolium", "veg_lin")
plot_abundance("Achillea.millefolium", "soil_lin")

## R markdown file with worked examples
file.show(system.file("doc/example-species-report.Rmd", package="cure4insect"))

## End(Not run)

ABbiodiversity/cure4insect documentation built on Oct. 11, 2022, 9:09 p.m.