SDLfilter (version 2.0.1)

boot_area: Cumulative analysis of collective areas by bootstrapping

Description

Function to calculate collective areas (merged x% Utilisation Distributions or UDs) of n individuals by bootstrapping.

Usage

boot_area(
  data,
  cell.size = NA,
  R = 1000,
  percent = 50,
  quantiles = c(0.25, 0.5, 0.75)
)

Arguments

data

A matrix or list of RasterLayer objects. Each row of the matrix or each RasterLayer object contains a utilisation distribution (or other statistics that sums to 1 - e.g. proportion of time spent). The grid size and geographical extent must be consistent across each row of the matrix or each RasterLayer object. The function assumes that each column of the matrix is associated with a unique geographical location or that each RasterLayer has exactly the same geographical extent and resolution.

cell.size

A numeric value specifying the grid cell size of the input data in metres.

R

An integer specifying the number of iterations.

percent

An integer specifying the percent volume of each UD to be considered in the analysis.

quantiles

A vector or a number to specify the quantiles to be calculated in the summary of the results.

Value

A list containing two data frames - raw results and summary (mean, sd, sem and quantiles at each sample size).

Details

This function calculates collective areas (e.g. 50% UDs) of 1 to n individuals by bootstrapping.

References

Shimada et al. (in prep)

Examples

# NOT RUN {
#1 Utilisation distributions of flatback turtles.
data(curtis)

#2 Calculate collective areas from 1000 random permutation
area <- boot_area(curtis, R = 1000, percent = 50)

#3 Find the minimum sample size required to estimate the general distribution.
a <- asymptote(area)

#4 Plot the mean collective area and rational function fit relative to the sample sizes.
ggplot(data = area$summary)+
  geom_point(aes(x = N, y = mu/1e+6), alpha = 0.5) + 
  geom_path(data = a$results, aes(x = x, y = ys/1e+6)) +        
  labs(x = "N", y = expression(Area~(km^2)))
# }