SK_EBLUP_LME.f: Evaluate fitted taper curve
In jonibio/TapeR: Flexible Tree Taper Curves Based on Semiparametric Mixed Models

Description Usage Arguments Details Value Author(s) See Also Examples

This is the actual function to estimate diameters according to the fitted mixed B-splines model.

1	SK_EBLUP_LME.f(xm, ym, xp, par.lme, R0 = FALSE, ...)

`xm`	relative heights for which measurements are available
`ym`	corresponding diameter measurements in height `xm`
`xp`	relative heights for which predictions are required
`par.lme`	Fitted model object, return of `TapeR_FIT_LME.f`
`R0`	boolean, should residual variance-covariance matrix be set to zero on off-diagonal elements? Defaults to FALSE. See details.

This function is the actual working horse for prediction using the fitted taper model. Based on the model par.lme and the measured diameters ym and corresponding (relative) heights xm of a specific tree (there might be just one measurement), the random effect parameters and subsequently diameters are estimated. Depending on the boolean parameter R0, the calibrated taper curve is forced through the given diameter ym (R0 = TRUE), or calibrated using the complete residual variance-covariance information (R0 = FALSE). See Kublin et al. (2013) p. 987 for more details.

a list holding nine elements:

b_fix fixed effects parameter of taper model
b_rnd random effects parameter given tree (posterior mean b_k)
yp estimated diameter in height xp
KOV_Mean variance-covariance matrix of expected value
KOV_Pred variance-covariance matrix of prediction
CI_Mean mean and limits of confidence interval
MSE_Mean mean squared error of expected value
MSE_Pred mean squared error of prediction
CI_Pred mean and limits of prediction interval

Edgar Kublin

E_DHx_HmDm_HT.f, E_VOL_AB_HmDm_HT.f

data("SK.par.lme")
TapeR:::SK_EBLUP_LME.f(1.3/27, 30, 1.3/27, SK.par.lme)
## using empirical best linear unbiased estimator, estimate != 30
TapeR:::SK_EBLUP_LME.f(1.3/27, 30, 1.3/27, SK.par.lme, R0=FALSE)$yp
## interpolate / force through given diameter, estimate  == 30
TapeR:::SK_EBLUP_LME.f(1.3/27, 30, 1.3/27, SK.par.lme, R0=TRUE)$yp
TapeR:::SK_EBLUP_LME.f(1.3/27, 30, c(1.3, 5)/27, SK.par.lme)