Nicholas Davies edited Method.tex  about 8 years ago

Commit id: 026ef608704016f74fa7585aac4724ded81fe191

deletions | additions      

       

\section{Materials and Methods}  An \textit{E. bosistoana} open-pollinated progeny trial was established at an irrigated nursery site in Harewood, Christchurch, New Zealand. The trial represented 40 families from two provenances, for a total of 423 seedlings planted into 100 L bags, which were coppiced after the first harvest, giving a total of 623 tested samples. Two separate plantings (or trial sections) occurred in 2010 and 2012. The 2010 families originated from South East Australia, were harvested and coppiced in 2012, and harvested again in December 2014. The 2012 families originated from higher altitude in New South Wales and were harvested in 2013 (this data was not included in the analysis, due to the magnitude of errors induced by small, malformed stems) and again in October 2015 (at age two). All seedlings were established following a completely randomised design.  When harvested, each sample was processed for growth strain, volumetric shrinkage (displacement method, before and after drying), stem diameter (measured under bark using digital calipers), basic density (mass and displacement method) and dry acoustic velocity (resonance). Growth-strain was measured using a modified version of \citet{Chauhan_2010} \cite{Chauhan_2010}  and \citet{Entwistle_2014} \cite{Entwistle_2014}  method. The newly developed "rapid splitting test" substantially reduces measurement time, enabling larger numbers of samples to be processed. The method involves stripping the bark and measuring the under-bark large-end diameter of a clear section of the stem, giving an over estimate of the average diameter used by \citet{Chauhan_2010}. \cite{Chauhan_2010}.  The sample is then cut lengthwise, with the slit length determined by the length of clear wood and diameter of the sample. Diameter and slit length are measured, recorded and the stem cut. The small-end of the sample is left intact with the large-end free to distort, which removes the need to clamp the two halves together. Finally, the opening is measured and recorded. The calculation of strain is unchanged (with the exception that average radius is now large-end radius) from \citet{Chauhan_2010} \cite{Chauhan_2010}  and calculated using equation \ref{gse}. It is important to note that the over estimate of radius slightly reduces the strain value, but does so linearly over all samples. \begin{equation}  \label{gse} 

  Where: \(Y_u\) is the deflection, \(\epsilon\) is the strain, \(L\) is the cut length and \(R\) is the big end cross-section radius.   Narrow-sense heritability for all properties was calculated using equation \ref{H2}. The constant of 2.5 used was suggested by \citep{Griffin1988Genetic} \cite{Griffin1988Genetic}  due to the unknown proportions of selfing, full-siblings and half-siblings within the open-pollinated families (loosely refereed to as half-sibling families in this article). \begin{equation}  \label{H2}  H^2 = \frac{2.5 \times Var(Family)}{Var(Family) + Var(Tree) + Var(Residual)}  \end{equation}  The analyses were conducted in R \citep{RCORE} \cite{RCORE}  and JAGS \citep{RJAGS}, \cite{RJAGS},  utilizing a Bayesian approach to estimate the posterior distributions for the heritability of growth-strain and other wood properties. The effect of coppicing and trial section were included as fixed effects, while half sibling family and the individual trees are included as random effects. Although all specimens were grown on the same site, they were grown during different time periods which are confounded with the effect of the two provenances and hence are included as the trial section fixed-effect. The tree effect accounts for the measurement unit, as the same trees were assessed as both seedlings and coppice. In addition, residuals were modeled separately for each section. Aproximatly 35\% of the strain data was left-censored, and the JAGS function \verb|dinterval| was used to simulate the truncated data, for details see \citet{lunn2012bugs}. \cite{lunn2012bugs}.