estHeight {TapeS}

Estimate tree height by DBH according to BWI3

Description

Function calculates expected tree height given diameter in breast height and species code

Usage

estHeight(d13, sp, qtl = NULL)

Arguments

d13

numeric vector of diameter in breast height [cm]

sp

TapeS species code, see also tprSpeciesCode

qtl

desired quantile of height, either NULL (LS-regression) or one of 17, 50, 83 (quantile regression)

Details

Function evaluates the expected value of a Pettersen-Height Function based on diameter in breast height and tree species code. The Pettersen-Function (h = 1.3 + (a + \frac{b}{dbh})^{-3}) was fitted on NFI 3 (BWI 3) data, using the main stand only.

d13 and sp should be of equal length or one of it can be > 1 if the other is of length 1. Then, the shorter object will be extended to match the length of the longer object. See examples.

The quantile option return tree height at quantiles 17, 50 or 83. If qtl is NULL, the result of a nonlinear least-squares regression is provided.

Value

a vector of same length as d13 or sp, with tree height in [m].

Examples

sp <- 1
d13 <- 30
estHeight(d13, sp)

sp <- 1
d13 <- seq(15, 50, 5)
estHeight(d13, sp)

sp <- 1:36
d13 <- 30
estHeight(d13, sp)