Modelling count data using the logratio-normal-multinomial distribution

The logratio-normal-multinomial distribution is a count data model resulting from compounding a multinomial distribution for the counts with a multivariate logratio-normal distribution for the multinomial event probabilities. However, the logratio-normal-multinomial probability mass function does no...

Descripción completa

Detalles Bibliográficos
Autores: Comas Cufí, Marc, Martín-Fernández, Josep Antoni, Mateu-Figueras, Glòria, Palarea-Albaladejo, Javier
Tipo de recurso: artículo
Fecha de publicación:2020
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Idioma:inglés
OAI Identifier:oai:upcommons.upc.edu:2117/362096
Acceso en línea:https://hdl.handle.net/2117/362096
https://dx.doi.org/10.2436/20.8080.02.96
Access Level:acceso abierto
Palabra clave:count data
compound probability distribution
Dirichlet multinomial
logratio coordinates
Monte Carlo method
simplex
Estadística matemàtica
Estadística matemàtica--Aplicacions
Classificació AMS::62 Statistics::62P Applications
Classificació AMS::62 Statistics::62F Parametric inference
Àrees temàtiques de la UPC::Matemàtiques i estadística::Estadística matemàtica
Descripción
Sumario:The logratio-normal-multinomial distribution is a count data model resulting from compounding a multinomial distribution for the counts with a multivariate logratio-normal distribution for the multinomial event probabilities. However, the logratio-normal-multinomial probability mass function does not admit a closed form expression and, consequently, numerical approximation is required for parameter estimation. In this work, different estimation approaches are introduced and evaluated. We concluded that estimation based on a quasi-Monte Carlo Expectation-Maximisation algorithm provides the best overall results. Building on this, the performances of the Dirichlet-multinomial and logratio-normal-multinomial models are compared through a number of examples using simulated and real count data.