TapeR (version 0.5.3)

SK_VOLab_EBLUP_LME.f: taper volume estimation

Description

Internal function not usually called by users

Usage

SK_VOLab_EBLUP_LME.f(
  xm,
  ym,
  a = 0,
  b = 1,
  Ht,
  par.lme,
  Rfn = list(fn = "sig2"),
  IntPolOpt = TRUE,
  ...
)

Value

List with two elements, the estimated volume and its variance

Arguments

xm

relative heights for which measurements are available

ym

corresponding diameter measurements in height xm

a

relative height of lower threshold of stem section

b

relative height of upper threshold of stem section

Ht

tree height

par.lme

List of taper model parameters obtained by TapeR_FIT_LME.f

Rfn

list with function name to provide estimated or assumed residual variances for the given measurements, optionally parameters for such functions

IntPolOpt

option for method of interpolation, if TRUE using a natural interpolating spline (splinefun), if FALSE using a smoothing spline (smooth.spline); defaults to TRUE

...

not currently used

Author

Edgar Kublin

Details

with Rfn=list(fn="zero") one can decide whether the measured diameters are forced to lie exactly on the taper curve; this interferes somewhat with the IntPolOpt, which determines the method of taper curve point interpolation for integration. The default TRUE (used throughout all function calls) applies natural interpolating splines, hence this does not contradict the optional use of Rfn=list(fn="zero").