Tsagris 2018 - Dirichlet Regression with Zeros - Overview

Summary

Tsagris and Stewart modify the Dirichlet regression log-likelihood so compositions containing exact zeros can be modelled without imputing them: each zero pattern is a “population” whose non-zero components follow a Dirichlet with the same coefficients and precision. The resulting ZADR model competes with Aitchison’s imputed log-ratio regression and the Kent/square-root model, and wins when zeros are numerous.

Overview

Question. Aitchison’s log-ratio regression and Dirichlet regression both fail when a component is exactly zero. Imputation (Palarea-Albaladejo and Martin-Fernandez 2008) assumes zeros are rounded values, alters every other component of an affected vector, and is sensitive to the assumed detection limit.

Contribution. ZADR (simple and full) built on Stewart and Field (2011) and Bear and Billheimer (2015), plus a four-step maximisation using the Aitchison fit as starting values.

Key result. Simple ZADR beats Aitchison+imputation on Arctic-lake data with injected zeros and on foraminifera (KL 3.126 vs 3.556), and in simulation reaches proportion 1 of runs with lower KL as grows.

Motivating PyMC examples

E70 (Aitchison’s Arctic lake sediments: sand/silt/clay vs depth; Dirichlet regression in mean-precision form) at https://redam94.github.io/pymc-challenges/examples.html; see Greenacre 2023 - Compositional Outcome Models in PyMC.

Main Content — Source Map

Relevance to marketing measurement

Share outcomes with structural zeros (a channel not run, a competitor absent in a market) are common. ZADR-style handling keeps observed shares unchanged rather than imputing an arbitrary epsilon.

Connections

Foundation in Greenacre 2023 - Wide Data, Zeros and Quasi-Coherence; Generalized Linear Models (multinomial-logit link).

See Also

Cluster notes under Greenacre 2023 - Aitchison CoDA 40 Years On - Overview.