Hierarchical linear mixed models in multi-stage sampling soil studies

Abstract

The issue of variances of different soil variables prevailing at different sampling scales is addressed. This topic is relevant for soil science, agronomy and landscape ecology. In multi-stage sampling there are randomness components in each stage of sampling which can be taken into account by introducing random effects in analysis through the use of hierarchical linear mixed models (HLMM). Due to the nested sampling scheme, there are several hierarchical sub-models. The selection of the best model can be carried out through likelihood ratio tests (LRTs) or Wald tests, which are asymptotically equivalent under standard conditions. However, when the comparison leads to a restricted hypothesis of variance components, standard conditions are not maintained, which leads to more elaborated versions of LRTs. These versions are not disseminated among environmental scientists. The present study shows the modeling of soil data from a sampling where sites, fields within sites, transects within fields, and sampling points within transects were selected in order to take samples from different vegetation types (open and shade). For soil data, several sub-models were compared using Wald tests, classic LRTs and adjusted LRTs where the distribution of the test statistic under the null hypothesis is the Chi-square mixture of Chi-square distributions. The inclusion of random effects via HLMM and suggested by the latest version of LRT allowed us to detect effects of vegetation type on soil properties that were not detected under a classical ANOVA.