Statistical Modelling

Assessing associations between a response of interest and a set of covariates in spatial areal models is the leitmotiv of ecological regression. However, the presence of spatially correlated random effects can mask or even bias estimates of such associations due to confounding effects if they are not carefully handled. Though potentially harmful, confounding issues have often been ignored in practice leading to wrong conclusions about the underlying associations between the response and the covariates. In spatio-temporal areal models, the temporal dimension may emerge as a new source of confounding, and the problem may be even worse. In this work, we propose two approaches to deal with confounding of fixed effects by spatial and temporal random effects, while obtaining good model predictions. In particular, restricted regression and an apparently—though in fact not—equivalent procedure using constraints are proposed within both fully Bayes and empirical Bayes approaches. The methods are compared in terms of fixed-effect estimates and model selection criteria. The techniques are used to assess the association between dowry deaths and certain socio-demographic covariates in the districts of Uttar Pradesh, India.

This article develops an approximate N-mixture model for infectious disease counts that accounts for under-reporting as well as spatial dependence induced by person-to-person spread of disease. We employ the model to estimate actual case counts in Oregon of chlamydia, an easily-treated but usually asymptomatic sexually transmitted disease. We describe a combined parametric bootstrap to account for uncertainty in parameter estimates as well as sampling variability in actual case counts. A simulation study illustrates that our method performs well in many scenarios when the model is correctly specified, and also gives reasonable results when the model is misspecified, and no spatial dependence exists.

It is known that the estimating equations for quantile regression (QR) can be solved using an EM algorithm in which the M-step is computed via weighted least squares, with weights computed at the E-step as the expectation of independent generalized inverse-Gaussian variables. This fact is exploited here to extend QR to allow for random effects in the linear predictor. Convergence of the algorithm in this setting is established by showing that it is a generalized alternating minimization (GAM) procedure. Another modification of the EM algorithm also allows us to adapt a recently proposed method for variable selection in mean regression models to the QR setting. Simulations show that the resulting method significantly outperforms variable selection in QR models using the lasso penalty. Applications to real data include a frailty QR analysis of hospital stays, and variable selection for age at onset of lung cancer and for riboflavin production rate using high-dimensional gene expression arrays for prediction.

There is an increasing interest in models for discrete valued time series. Among them, the integer autoregressive conditional heteroscedastic (INGARCH) is a model that has found several applications. In the present article, we study the problem of model selection for this family of models. Namely we consider that an observation conditional on the past follows a Poisson distribution where its mean depends on its past mean values and on past observations. We consider both linear and log-linear models. Our purpose is to select the most appropriate order of such models, using a trans-dimensional Bayesian approach that allows jumps between competing models. A small simulation experiment supports the usage of the method. We apply the methodology to real datasets to illustrate the potential of the approach.