GapSizeFDist: Forest Canopy Gap-size Frequency Distributions

Description Usage Arguments Value Details Author(s) References Examples

View source: R/GapSizeFDist.R

Description

This function quantifies the size frequency distribution of forest canopy gaps using Zeta distribution, which is a discrete power-law probability density.

Usage

1
GapSizeFDist(gaps_stats,...)

Arguments

gaps_stats

A data.frame containing basic statistics of forest gaps. Output of (GapStats) function.

...

Supplementary parameters for (plot).

Value

A log-log plot of gap-size Frequency Distributions and a vector containing the λ and minimum value of the likelihood. The parameter λ is the scaling exponent for the power-law Zeta distribution fitted to the data using maximum likelihood. See details section.

Details

For the Zeta distribution with parameter λ, the probability that gap size takes the integer value k is:

f(k) = k^-λ/ζ(λ)

where the denominator is the Riemann zeta function, and is undefined for λ = 1. The function calculates maximum likelihood estimates (MLE) of λ by minimizing a negative log-likelihood function (Asner et 2013).

Author(s)

Carlos Alberto Silva. This function was implemented based on the R script provided by Asner et al. (2013)

References

Asner, G.P., Kellner, J.R., Kennedy-Bowdoin, T., Knapp, D.E., Anderson, C. & Martin, R.E. (2013). Forest canopy gap distributions in the Southern Peruvian Amazon. PLoS One, 8, e60875.

White EP, Enquist BJ, Green JL (2008) On estimating the exponent of powerlaw frequency distributions. Ecology 89: 905<e2><80><93>912.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
#Loading raster library
library(raster)

# ALS-derived CHM over Adolpho Ducke Forest Reserve - Brazilian tropical forest
data(ALS_CHM_DUC)

# set height tresholds (e.g. 10 meters)
threshold<-10
size<-c(1,10^4) # m2

# Detecting forest gaps
gaps_duc<-getForestGaps(chm_layer=ALS_CHM_DUC, threshold=threshold, size=size)

# Computing basic statistis of forest gap
gaps_stats<-GapStats(gap_layer=gaps_duc, chm_layer=ALS_CHM_DUC)

# Gap-size Frequency Distributions
GapSizeFDist(gaps_stats=gaps_stats, col="forestgreen", pch=16, cex=1,
axes=FALSE,ylab="Gap Frequency",xlab=as.expression(bquote("Gap Size" ~ (m^2) )))
axis(1);axis(2)
grid(4,4)

ForestGapR documentation built on May 2, 2019, 1:26 a.m.