Cohesive framework for modelling plant cover class data

The study of plant distribution and abundance is a fundamental pursuit in ecology and conservation biology. Measuring plant abundance by visually assessing percent cover and recording a cover class is a common field method that yields ordinal data. Statistical models for ordinal data exist but entail cumbersome interpretations and sometimes restrictive assumptions. We propose a Bayesian hierarchical framework for analysing cover class data that allows for linking ordinal observations to a latent beta distribution and accounts for zero inflation. Harnessing a latent beta distribution supports interpreting changes in abundance in terms of mean percent cover rather than odds ratios of cumulative cover classes as for cumulative link models. The zero augmentation allows for simultaneous inferences on both occurrence (distribution) and abundance. We show how our model can account for true and false zeros, misclassification of cover classes, multiple species and hierarchical sampling designs, using empirical examples and simulations. Simulated observation errors, when ignored, led to models overestimating abundance and underestimating occurrence. Based on simulations, we found no substantial difference between mean percent cover estimates when analyzing ordinal cover classes versus continuous percent cover as the response. Our empirical datasets displayed high probability of detection (>0.85 on average for all species), likely due to the sampling design used and training of observers. Probability of occurrence was slightly underestimated for bare ground, Artemisia tridentata, Elycap medusae, and Poa secunda using a model that ignored imperfect detection. Estimated mean percent cover was not substantially impacted by ignoring measurement error for five plant species and bare ground. Our modelling framework for cover class data allows for an explicit separation of distribution from abundance and, importantly, allows for interpreting species-environment relationships in terms of variation in mean percent cover as compared to cumulative odds ratios. The beta distribution inherently accommodates heteroscedasticity and skewness, statistical properties that are a consequence of spatially aggregated patterns common to plant survey data. Recording cover classes provides a reliable, efficient way to measure plants and our simulations suggest little loss of information compared to assuming continuous percent cover. We provide JAGS and Stan model code for implementation.