R/alemdag.R

Defines functions alemdag

Documented in alemdag

#' Alemdag's (1978) tree competition index
#' 
#' @param dbh vector of DBH (diameter at breast height) measurements of competitor trees
#' @param dist vector of distances from focal tree to competitor trees
#' @param focal_dbh DBH of focal tree
#'
#' @return value of competition index for focal tree
#'
#' @details A spatially explicit competition index originally used in white 
#' spruce plantations. The value of this function increases when competitor 
#' trees are closer to the focal tree, or when competitor trees are larger. 
#' Alemdag's competition index is defined by the following equation:
#' \deqn{\sum^{n}_{j\neq{}i}(\pi[(l_{ij}d_{i})/(d_{i} + d_{j})]^2 (d_{j} / l_{ij}) / \sum(d_{j} / l_{ij}))}
#' where \eqn{l_{ij}} is the distance between focal tree \eqn{i} and competitor tree 
#' \eqn{j}, and \eqn{d_{i}} is the diameter (DBH) of focal tree \eqn{i}.
#' 
#' @references Alemdag I. S. (1978). Evaluation of some competition indexes for 
#' the prediction of diameter increment in planted white spruce. Forest 
#' Management Institute, Ottawa, Canada.
#'
#' @examples
#' data(bicuar)
#' nb <- nearNeighb(bicuar$x, bicuar$y, bicuar$stem_id, k = 4) 
#' lapply(nb, function(x) {
#'   nb <- merge(x, bicuar, by.x = "nb", by.y = "stem_id")
#'   focal_diam <- unique(bicuar[bicuar$stem_id == nb$focal,"diam"])
#'   alemdag(nb$diam, nb$nb_dist, focal_diam)
#'   })
#' 
#' @export
#'
alemdag <- function(dbh, dist, focal_dbh) {
	n <- length(dbh)
	eq1 <- dist * (dbh / (dbh + focal_dbh))
	eq2 <- dbh / dist
	
	n / sum(pi * eq1^2 * (eq2 / sum(eq2)))
}
johngodlee/compInd documentation built on Aug. 5, 2024, 8:44 a.m.