E_HDxoR_HmDm_Ht.f {TapeS}

Find Height of diameter under bark via uniroot

Description

Functional equivalent to E_HDx_HmDm_HT.f, finding the height of a given diameter *without* bark, i.e. double bark thickness needs to be added on top of given diameter to find appropriate height.

Usage

E_HDxoR_HmDm_Ht.f(DxoR, Hm, Dm, mHt, sHt = 0, par.lme, Rfn = NULL, ...)

HxoR_root.f(Hx, DxoR, Hm, Dm, mHt, sHt, par.lme, Rfn, ...)

Arguments

DxoR

Scalar. Diameter under bark for which to return height.

Hm

Numeric vector of stem heights (m) along which diameter measurements were taken for calibration. Can be of length 1. Must be of same length as Dm.

Dm

Numeric vector of diameter measurements (cm) taken for calibration. Can be of length 1. Must be of same length as Hm.

mHt

Scalar. Tree height (m).

sHt

Scalar. Standard deviation of stem height. Can be 0 if height was measured without error.

par.lme

List of taper model parameters obtained by TapeR_FIT_LME.f, enhanced by the attribute 'spp', which refers to the tree species used for double bark thickness

Rfn

setting for residuals error matrix, defaults to "sig2", see details.

...

not currently used

Hx

height at which taper curve is evaluated

Details

finds height of given diameter via uniroot.

Value

A scalar. Estimated height (m) given a diameter without bark.

Functions

• HxoR_root.f(): function to be searched

Examples

tmp <- tprTrees()
spp <- spp(tmp)
Hm <- Hm(tmp)
Dm <- Dm(tmp)
H <- Ht(tmp)
data("SKPar")
sppSK <- BaMap(spp, 1) # tree species for taper curve
## diameter in 5m height
TapeR::E_DHx_HmDm_HT.f(c(5, 10), Hm, Dm, mHt=H, sHt = 0, par.lme = SKPar[[sppSK]])$DHx
(D5m <- TapeR::E_DHx_HmDm_HT.f(c(5, 10), Hm, Dm, mHt=H, sHt = 0, par.lme = SKPar[[sppSK]])$DHx)
## bark thickness of diameter in 5m height
(RiD5m <- bark(c(1,1), Dm = D5m, relH = c(5, 10)/H))
## find height of diameter without bark, which should be 5m
d5mub <- D5m - RiD5m
E_HDxoR_HmDm_Ht.f(DxoR = d5mub, Hm = Hm, Dm = Dm, mHt = H,
                  sHt = 0, par.lme = SKPar[[sppSK]])