HDmodels: Available 2- and 3- parameter H-D model functions to be used...

Description Usage Arguments Details Value Author(s) References Examples

Description

Nonlinear functions for modeling tree height on diameter. Usually called using fithd.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
HDnaslund(d, a, b, bh=1.3) 
HDcurtis(d, a, b, bh=1.3) 
HDmichailoff(d, a, b, bh=1.3) 
HDmeyer(d, a, b, bh=1.3) 
HDpower(d, a, b, bh=1.3)
HDnaslund2(d, a, b, bh=1.3)
HDnaslund3(d, a, b, bh=1.3)
HDnaslund4(d, a, b, bh=1.3)
HDmicment(d, a, b, bh=1.3) 
HDmicment2(d, a, b, bh=1.3) 
HDwykoff(d, a, b, bh=1.3) 


HDprodan(d, a, b, c, bh=1.3) 
HDlogistic(d, a, b, c, bh=1.3) 
HDrichards(d, a, b, c, bh=1.3) 
HDweibull(d, a, b, c, bh=1.3) 
HDgomperz(d, a, b, c, bh=1.3) 
HDsibbesen(d, a, b, c, bh=1.3) 
HDkorf(d, a, b, c, bh=1.3) 
HDratkowsky(d, a, b, c, bh=1.3) 
HDhossfeldIV(d, a, b, c, bh=1.3)

startHDnaslund(d, h, bh=1.3) 
startHDcurtis(d, h, bh=1.3) 
startHDmichailoff(d, h, bh=1.3) 
startHDmeyer(d, h, bh=1.3) 
startHDpower(d, h, bh=1.3)
startHDnaslund2(d, h, bh=1.3) 
startHDnaslund3(d, h, bh=1.3) 
startHDnaslund4(d, h, bh=1.3) 
startHDmicment(d, h, bh=1.3) 
startHDmicment2(d, h, bh=1.3) 
startHDwykoff(d, h, bh=1.3) 

startHDprodan(d, h, bh=1.3) 
startHDlogistic(d, h, bh=1.3) 
startHDrichards(d, h, bh=1.3, b=0.04) 
startHDweibull(d, h, bh=1.3) 
startHDgomperz(d, h, bh=1.3) 
startHDsibbesen(d, h, bh=1.3, a=0.5) 
startHDkorf(d, h, bh=1.3) 
startHDratkowsky(d, h, bh=1.3, c=5) 
startHDhossfeldIV(d, h, bh=1.3, c=5)

Arguments

d

A vector of tree diameters, usually in cm

h

A vector of tree heights, usually in m. The observed heights should be always above or equal to bh.

a, b, c

Parameters a, b (and c for 3- parameter functions) of the applied function. See details for expressions of different functions.

bh

The applied height for the measurement of tree diameter (so called breast height). Of the same unit as h.

Details

The available 2- parameter functions are

The available 3- parameter functions are

For each model, two functions are provided: one computing the value of the H-D model for given diameters using given values of parameters a, b (and c), and another returning the initial guesses of a, b (and c) for given h-d data.

The initial guesses are in most cases computed by fitting a linearized version of the model into the provided h-d data using lm. For some 3- parameter versions, no straightforward linearization is possible and one of the parameters is set to a fixed sensible constant. Those values can be seen as additional arguments in the corresponding startHD - functions. Details can be seen directly from the function definitions.

The user can define her own functions to be used with fithd. The case-sensitive naming of the functions should follow exactly the naming convention shown above. In addition, the names of the of arguments, as well as their order, should be the same as in the functions above.

The models are named according to references in

Suggestions on naming and references on the functions are welcome.

Value

For functions HDxxx, a vector of tree heights corresponding diameters d is returned. For functions startHDxxx, a named vector of initial estimates of a, b and (c).

Author(s)

Lauri Mehtatalo <lauri.mehtatalo@uef.fi>

References

Mehtatalo, L., Gregoire, T.G., and de Miguel, S. Modeling Height-diameter curves for height prediction. Canadian Journal of Forest Research, 45(7): 826-837, doi: 10.1139/cjfr-2015-0054

Examples

1
2
3
4
5
data(spati)
theta<-startHDnaslund(spati$d,spati$h)
plot(spati$d,spati$h)
d<-seq(0,50)
lines(d,HDnaslund(d,a=theta[1],b=theta[2]),col="red",lwd=5)

Example output

Loading required package: stats4
Loading required package: nlme
Loading required package: spatstat
Loading required package: spatstat.data
Loading required package: rpart

spatstat 1.59-0       (nickname: 'J'ai omis les oeufs de caille') 
For an introduction to spatstat, type 'beginner' 


Note: R version 3.4.4 (2018-03-15) is more than 9 months old; we strongly recommend upgrading to the latest version

lmfor documentation built on July 8, 2020, 5:49 p.m.