Estimation of wood fibre length distributions from censored data through an EM algorithm

An expectation maximization (EM) algorithm is proposed to find fibre length distributions in standing trees. The available data come from cylindric wood samples (increment cores). The sample contains uncut fibres as well as fibres cut once or twice. The sample contains not only fibres, but also other cells, the so-called 'fines'. The lengths are measured by an automatic fibre-analyser, which is not able to distinguish fines from fibres and cannot tell if a cell has been cut. The data thus come from a censored version of a mixture of the fine and fibre length distributions in the tree. The parameters of the length distributions are estimated by a stochastic version of the EM algorithm, and an estimate of the corresponding covariance matrix is derived. The method is applied to data from northern Sweden. A simulation study is also presented. The method works well for sample sizes commonly obtained from increment cores.