Weibull distribution for characterising stem-diameter structure

The Weibull density distribution is known as

Latex formula,

where Latex formula is the location, Latex formula the scale and Latex formula is the shape parameter, i.e. the parameters of the Weibull distribution are interpretable, which is always a good property of models.
The cumulative distribution function, i.e. the integral of density function, is much simpler:

Latex formula

What does it mean?

The Weibull density distribution allows characterising tree stem-diameter distributions by providing trend curves but more importantly by summarising stem diameters by means of three parameters that can be interpreted.

The shape parameter of the Weibull distribution, Latex formula is of particular interest in this context and the following interpretation aid can be used (Burkhart and Tomé, 2012, p. 198):

Latex formula Skewed to the right
Latex formula Skewed to the left
Latex formula Negative exponential, reverse j-shaped

When Latex formula is less than 1, the distribution is reverse j-shaped found in uneven-aged forest stands and when Latex formula equals 1, a negative exponential distribution results. If Latex formula, the Weibull distribution approximates a normal distribution and this value divides left- and right-skewed curves. In general, Latex formula gives bell-shapes typical of even-aged forest stands. The location parameter is directly related to the minimum diameter in a stand (Burkhart and Tomé, 2012, p. 265 ). In the context of diameter distributions, all model parameters must be positive.

Histogram of example dbh data overlaid by Weibull density functions with different shape parameters but the same location and scale parameters.

Where does it come from?

The Weibull distribution is named after Swedish mathematician Waloddi Weibull, who described it in detail in 1951, although it was first identified by Fréchet (1927) and first applied by Rosin and Rammler (1933) to describe a particle size distribution.

Why is it important?

The Weibull distribution is one of the most flexible and most commonly applied model for tree stem-diameters. The model parameters are comparatively easy to estimate and the distribution is easy to apply.

How can it be estimated?

Robinson and Hamann (2010, p. 164ff.) described in detail how the three parameters of the Weibull distribution can be estimated using R and the maximum-likelihood method. The grey box below gives an adaptation of that method:

dweibull3 <- function(x, gamma, beta, alpha) {
   (gamma/beta) * ((x -alpha)/beta)^(gamma - 1) * (exp(-((x - alpha)/beta)^gamma))
 }

loss.w3 <- function(p, data)
   sum(log(dweibull3(data, p[1], p[2], p[3])))
 
mle.w3.nm <- optim(c(gamma = 1, beta = 5, alpha = 10),
                    loss.w3, data = myData$dbh, hessian = TRUE,
                    control = list(fnscale = -1)) 

mle.w3.nm$par # model parameters

# Check whether the curve looks OK
xx <- seq(10, 60, 1)
hist(myData$dbh, freq = F, breaks = 50, xlim = c(10, 60))
lines(xx, dweibull3(xx, mle.w3.nm$par[1], mle.w3.nm$par[2],
mle.w3.nm$par[3]), lty = 1, col = "red")

Since R only provides implementations for the two-parameter version, the code starts with a new function dweibull3() implementing the three-parameter version. This is followed by a maximum-likelihood loss function using the previously defined function dweibull3(). This loss function in turn now forms one of the arguments of the optim() function used for carrying out the regression. myData is a data frame that includes a vector of stem diameters that can be addressed by myData$dbh.

An alternative to nonlinear maximum-likelihood based regression are percentile estimations. The idea of this approach is to estimate the three parameters of the Weibull distribution from selected points of the distribution, e.g. the 63rd or 95th percentile. The theory of percentile estimators is well explained in Clutter et al. (1983, p. 127ff.). However, percentile estimations can be very valuable where maximum-likelihood regression does not produce any or unsuitable results. Some methods of percentile estimation have also been linked to straightforward sampling methods, so that the parameters of the Weibull function can almost be directly sampled in the field without much effort. Also, percentile methods can be used to identify starting values for nonlinear regression. Below one example method of percentile estimations is provided.

According to Wenk et al. (1990, p. 198f.) and Burkhart and Tomé (2012, p. 265)

Latex formula.

Latex formula and can be interpreted as the diameter where approximately 63% of all trees are smaller in diameter. Latex formula , i.e. the minimum diameter in a tree population or forest stand.

Finally, Gerold(1988) suggested that Latex formula can be estimated from Latex formula and the diameter where approximately 95\% of all trees are smaller in diameter.

Latex formula

Latex formula, Latex formula and Latex formula can be estimated from any empirical diameter distribution but also by employing a simple sampling procedure. Based on a systematic sampling grid approximately ten sample points need to be identified in every forest stand along with the first twelve tree neighbours nearest to each sample point. Römisch (1983) found that Latex formula can be estimated from the diameter of the fifth largest tree out of twelve sample trees and that Latex formula can be estimated from the largest diameter tree out of ten sample trees nearest to the sample point. Latex formula and Latex formula are then calculated as the arithmetic means of all ten samples. Latex formula is the smallest diameter of all sample trees (Gerold, 1988).

Literature

Burkhart, H. and Tomé, M., 2012. Modeling forest trees and stands. Springer, Dordrecht.

Clutter, J. L, Fortson, J. C., Pienaar, L. V., Brister, G. H. and Bailey, R. L., 1983. Timber management. A quantitative approach. John Wiley & Sons, New York.

Fréchet, M., 1927. Sur la loi de probabilité de l’écart maximum. Annales de la Société Polonaise de Mathematique 6: 93-116.

Gerold, D., 1988. Describing stem diameter structure and its development by using the Weibull distribution. Wissenschaftliche Zeitschrift der Technischen Universität Dresden 37: 221-224.

Rosin, P. and Rammler, E., 1933. The laws governing the fineness of powdered Coal. Journal of the Institute of Fuel 7: 29–36.

Robinson , A. P. and Hamann , J. D., 2010. Forest analytics with R. An introduction. Use R! Springer. New York, 339p.

Römisch, K., 1983. A mathematical model for simulating growth and thinnings of even-aged pure stands. PhD thesis Technical University Dresden. Dresden, 197p.

Wenk, G., Antanaitis, V. and Šmelko, Š ., 1990. Forest growth and yield science. Deutscher Landwirtschaftsverlag. Berlin, 448p.

Weibull, W., 1951. A statistical distribution function of wide application. Journal of Applied Mechanics 18: 293–297.

2 thoughts on “Weibull distribution for characterising stem-diameter structure

  1. Hej, I tried this code on my dataset, but had this error, which I have no idea how to fix:

    > dweibull3
    > loss.w3
    > mle.w3.nm <- optim(c(gamma = 1, beta = 5, alpha = 10),
    + loss.w3, data = p4$dia, hessian = TRUE,
    + control = list(fnscale = -1))
    Error in fn(par, …) : unused argument (data = p4$dia)

    • dweibull3 <- function(x, gamma, beta, alpha) {
      (gamma/beta) * ((x -alpha)/beta)^(gamma – 1) * (exp(-((x – alpha)/beta)^gamma))
      }

      loss.w3 <- function(dia, p3)
      sum(log(dweibull3(data, p[1], p[2], p[3])))

      mle.w3.nm <- optim(c(gamma = 1, beta = 5, alpha = 10),
      loss.w3, data = p4$dia, hessian = TRUE,
      control = list(fnscale = -1))

Leave a Reply

Your email address will not be published. Required fields are marked *