deletions | additions       diff --git a/Method.tex b/Method.tex  index 5db429b..2ceb5af 100644  --- a/Method.tex  +++ b/Method.tex 
...
\end{equation}  where the $\sigma^2$ variances were obtained as the reciprocal of the precisions ($\tau$) mentioned in the model description.  All the models were fitted using rjags, an R \citep{RCORE} interface to JAGS (Just Another Gibbs Sampler, \citep{RJAGS}) which uses Gibbs Sampling to estimate the marginal posterior distributions for the parameters of interest. Approximately 35\% of the strain data was left-censored, and missing values due to censoring were imputed by JAGS as a random value below the limit of detection (using the function \verb|dinterval|; for details see \citet{lunn2012bugs}), but above -1.5 (equivalent to a negative strain of 1500 micro-strains). We used a family model instead of an animal model as, for the purposes of a simple one-generation pedigree the models are equivalent (\textit{sensu} identical expected values and variance, \cite{henderson1985equivalent}), \citep{henderson1985equivalent}),  while the family model is much less computationally intensive. The R/JAGS code is available as supplementary material.