Search

In this article, we present the novel approach of using a multi-state model to describe longitudinal changes in cognitive test scores. Scores are modelled according to a truncated Poisson distribution, conditional on survival to a fixed endpoint, with the Poisson mean dependent upon the baseline score and covariates. The model provides a unified treatment of the distribution of cognitive scores, taking into account baseline scores and survival. It offers a simple framework for the simultaneous estimation of the effect of covariates modulating these distributions, over different baseline scores. A distinguishing feature is that this approach permits estimation of the probabilities of transitions in different directions: improvements, declines and death. The basic model is characterised by four parameters, two of which represent cognitive transitions in survivors, both for individuals with no cognitive errors at baseline and for those with non-zero errors, within the range of test scores. The two other parameters represent corresponding likelihoods of death. The model is applied to an analysis of data from the Canadian Study of Health and Aging (1991-2001) to identify the risk of death, and of changes in cognitive function as assessed by errors in the Modified Mini-Mental State Examination. The model performance is compared with more conventional approaches, such as multivariate linear and polytomous regressions. This model can also be readily applied to a wide variety of other cognitive test scores and phenomena which change with age., (© The Author(s) 2011 Reprints and permissions: sagepub.co.uk/journalsPermissions.nav.)

In this article, we consider methods for Bayesian computation within the context of brain imaging studies. In such studies, the complexity of the resulting data often necessitates the use of sophisticated statistical models; however, the large size of these data can pose significant challenges for model fitting. We focus specifically on the neuroelectromagnetic inverse problem in electroencephalography, which involves estimating the neural activity within the brain from electrode-level data measured across the scalp. The relationship between the observed scalp-level data and the unobserved neural activity can be represented through an underdetermined dynamic linear model, and we discuss Bayesian computation for such models, where parameters represent the unknown neural sources of interest. We review the inverse problem and discuss variational approximations for fitting hierarchical models in this context. While variational methods have been widely adopted for model fitting in neuroimaging, they have received very little attention in the statistical literature, where Markov chain Monte Carlo is often used. We derive variational approximations for fitting two models: a simple distributed source model and a more complex spatiotemporal mixture model. We compare the approximations to Markov chain Monte Carlo using both synthetic data as well as through the analysis of a real electroencephalography dataset examining the evoked response related to face perception. The computational advantages of the variational method are demonstrated and the accuracy associated with the resulting approximations are clarified.

Obtaining accurate measures of exposure to forest fire smoke is important for the assessment of health risk. Estimating exposure from air quality monitors is challenging because of the sparseness of the monitoring networks in remote areas. However, satellite imagery offers a novel and data-rich tool to provide visual information on smoke plumes. We will discuss statistical techniques for obtaining estimates of forest fire smoke plumes using classification algorithms on data from satellite imagery in order to develop automated processes for identifying exposure. The aim is to identify whether such methods may offer a high-resolution approach that provides a reliable estimate of smoke and a more thorough caption of the spatial distribution of smoke from fires than is currently available.