Search PubMed⌕ Search

Biomedical subjects

Ying C MacNab

Publications and source records attributed to Ying C MacNab.

8 recordsLinked to original sources

An innovative application of Bayesian disease mapping methods to patient safety research: a Canadian adverse medical event study.

Recently developed disease mapping and ecological regression methods have become important techniques in studies of disease epidemiology and in health services research. This increase in importance is partially a result of the development of Bayesian statistical methodologies that make it possible to study associations between health problems and risk factors at an aggregate (i.e. areal) level while taking into account such matters as unmeasured confounding and spatial relationships. In this paper we present a demonstration of the joint use of empirical Bayes (EB) and full Bayesian inferential techniques in a small area study of adverse medical events (also known as 'iatrogenic injury') in British Columbia, Canada. In particular, we illustrate a unified Bayesian hierarchical spatial modelling framework that enables simultaneous examinations of potential associations between adverse medical event occurrence and regional characteristics, age effects, residual variation and spatial autocorrelation. We propose an analytic strategy for complementary use of EB and FB inferential techniques for risk assessment and model selection, presenting an EB-FB combined approach that draws on the strengths of each method while minimizing inherent weaknesses. The work was motivated by the need to explore relatively efficient ways to analyse regional variations of health services outcomes and resource utilization when a considerable amount of statistical modelling and inference are required.

Adolescent↗

SNAP-II predicts mortality among infants with congenital diaphragmatic hernia.

OBJECTIVE: Outcomes analysis in congenital diaphragmatic hernia (CDH) requires a validated risk-adjustment tool. The purpose of this study was to use the Canadian Neonatal Network (CNN) database to validate the Score for Neonatal Acute Physiology, Version II (SNAP-II) for prediction of mortality among CDH infants admitted to a neonatal intensive care unit (NICU), and to compare this to the predictive equation recently developed by the Congenital Diaphragmatic Hernia Study Group (CDHSG). STUDY DESIGN: Infants with CDH in the CNN database were identified. Bivariate and multivariable logistic regression models were used to identify risk factors predictive of mortality. Model predictive performance and calibration were assessed using the area under the receiver operator characteristic curve and the technique of Hosmer-Lemeshow, respectively, and compared with the CDHSG predictive equation. RESULTS: There were 88 patients with CDH among 19,507 admissions to CNN hospitals. The mortality rate among CDH patients surviving to NICU admission was 17%, and 12.5% received extracorporeal membrane oxygenation therapy. Gestational age and admission SNAP-II score predicted mortality. Model predictive performance and calibration were optimized with these variables combined. The CDHSG equation was equally predictive of mortality, but was only marginally calibrated. CONCLUSIONS: SNAP-II is highly predictive of mortality among patients with CDH, and can be used to risk-adjust these patients.

Cause of Death↗

Bayesian spatial and ecological models for small-area accident and injury analysis.

In this article, recently developed Bayesian spatial and ecological regression models are applied to analyse small-area variation in accident and injury. This study serves to demonstrate how Bayesian modelling techniques can be implemented to assess potential risk factors measured at group (e.g. area) level. Presented here is a unified modelling framework that enables thorough investigations into associations between injury rates and regional characteristics, residual variation and spatial autocorrelation. Using hospital separation data for 83 local health areas in British Columbia (BC), Canada, in 1990-1999, we explore and examine ecological/contextual determinants of motor vehicle accident injury (MVAI) among male children and youth aged 0-24 and for those of six age groups (<1, 1-4, 5-9, 10-14, 15-19 and 20-24). Eighteen local health area characteristics are studied. They include a broad spectrum of socio-economic indicators, residential environment indicators (roads and parks), medical services availability and utilisation, population health, proportion of recent immigrants, crime rates, rates of speeding charge and rates of seatbelt violation. Our study indicates a large regional variation in MVAI in males aged 0-24 in British Columbia, Canada, in 1990-1999, and that adjusting for appropriate risk factors eliminates nearly all the variation observed. Socio-economic influence on MVAI was profoundly apparent in young males of all ages with the injury being more common in communities of lower socio-economic status. High adult male crime rates were significantly associated with high injury rates of boys aged 1-14. Seatbelt violations and excess speeding charges were found to be positively associated with the injury rates of young men aged 20-24. This and similar ecological studies shed light on reasons for regional variations in accident occurrence as well as in the resulting injuries and hospital utilisation. Thereby they are potentially useful in identifying priority areas for injury/accident prevention and in informing regional health planning and policy development.

Accidents↗

Estimation in Bayesian disease mapping.

Recent work on Bayesian inference of disease mapping models discusses the advantages of the fully Bayesian (FB) approach over its empirical Bayes (EB) counterpart, suggesting that FB posterior standard deviations of small-area relative risks are more reflective of the uncertainty associated with the relative risk estimation than counterparts based on EB inference, since the latter fail to account for the variability in the estimation of the hyperparameters. In this article, an EB bootstrap methodology for relative risk inference with accurate parametric EB confidence intervals is developed, illustrated, and contrasted with the hyperprior Bayes. We elucidate the close connection between the EB bootstrap methodology and hyperprior Bayes, present a comparison between FB inference via hybrid Markov chain Monte Carlo and EB inference via penalized quasi-likelihood, and illustrate the ability of parametric bootstrap procedures to adjust for the undercoverage in the "naive" EB interval estimates. We discuss the important roles that FB and EB methods play in risk inference, map interpretation, and real-life applications. The work is motivated by a recent analysis of small-area infant mortality rates in the province of British Columbia in Canada.

Algorithms↗

Hierarchical Bayesian spatial modelling of small-area rates of non-rare disease.

We present Bayesian hierarchical spatial models for the analysis of the geographical distribution of a non-rare disease or event. The work is motivated by the need for ascertaining regional variations in health services outcomes and resource use and for assessing the potential sources of these variations. The models discussed herein readily accommodate random spatial effects and covariate effects. We discuss Bayesian inferential framework and implementation of a hybrid Markov chain Monte Carlo method for full Bayesian model inference. The methods are illustrated through an analysis of regional variation in chronic lung disease (CLD) rates among neonatal intensive care unit (NICU) patients across Canada. Specifically, we first present a random effects binomial model for spatially correlated CLD rates, with random spatial effects accounting for latent or covariate effects. These random spatial effects depict regional or spatial variation in chronic lung disease occurrence. We then extend this model to include covariates. With this extension, we assess residual spatial effects and the extent to which risk factors such as illness severity at NICU admission, low birth weight, and very low birth weight influence the CLD rate variation.

Bayes Theorem↗

A Bayesian hierarchical model for accident and injury surveillance.

This article presents a recent study which applies Bayesian hierarchical methodology to model and analyse accident and injury surveillance data. A hierarchical Poisson random effects spatio-temporal model is introduced and an analysis of inter-regional variations and regional trends in hospitalisations due to motor vehicle accident injuries to boys aged 0-24 in the province of British Columbia, Canada, is presented. The objective of this article is to illustrate how the modelling technique can be implemented as part of an accident and injury surveillance and prevention system where transportation and/or health authorities may routinely examine accidents, injuries, and hospitalisations to target high-risk regions for prevention programs, to evaluate prevention strategies, and to assist in health planning and resource allocation. The innovation of the methodology is its ability to uncover and highlight important underlying structure of the data. Between 1987 and 1996, British Columbia hospital separation registry registered 10,599 motor vehicle traffic injury related hospitalisations among boys aged 0-24 who resided in British Columbia, of which majority (89%) of the injuries occurred to boys aged 15-24. The injuries were aggregated by three age groups (0-4, 5-14, and 15-24), 20 health regions (based of place-of-residence), and 10 calendar years (1987 to 1996) and the corresponding mid-year population estimates were used as 'at risk' population. An empirical Bayes inference technique using penalised quasi-likelihood estimation was implemented to model both rates and counts, with spline smoothing accommodating non-linear temporal effects. The results show that (a) crude rates and ratios at health region level are unstable, (b) the models with spline smoothing enable us to explore possible shapes of injury trends at both the provincial level and the regional level, and (c) the fitted models provide a wealth of information about the patterns (both over space and time) of the injury counts, rates and ratios. During the 10-year period, high injury risk ratios evolved from northwest to central-interior and the southeast [corrected].

Accidents, Traffic↗

Hierarchical Bayesian modeling of spatially correlated health service outcome and utilization rates.

We present Bayesian hierarchical spatial models for spatially correlated small-area health service outcome and utilization rates, with a particular emphasis on the estimation of both measured and unmeasured or unknown covariate effects. This Bayesian hierarchical model framework enables simultaneous modeling of fixed covariate effects and random residual effects. The random effects are modeled via Bayesian prior specifications reflecting spatial heterogeneity globally and relative homogeneity among neighboring areas. The model inference is implemented using Markov chain Monte Carlo methods. Specifically, a hybrid Markov chain Monte Carlo algorithm (Neal, 1995, Bayesian Learning for Neural Networks; Gustafson, MacNab, and Wen, 2003, Statistics and Computing, to appear) is used for posterior sampling of the random effects. To illustrate relevant problems, methods, and techniques, we present an analysis of regional variation in intraventricular hemorrhage incidence rates among neonatal intensive care unit patients across Canada.

Bayes Theorem↗

Spatio-temporal modelling of rates for the construction of disease maps.

There have been significant developments in disease mapping in the past few decades. The continual development of statistical methodology in this area is responsible for the growing popularity of disease mapping because of its potential usefulness in regional health planning, disease surveillance and intervention, and allocating health funding. Here we review the area of disease mapping where relative risks pertain to an event such as incidence or mortality over space and time. In particular we briefly discuss the use of generalized additive mixed models, an additive extension of generalized linear mixed models, for spatio-temporal analysis of disease rates. To illustrate the procedures, we present an in-depth analysis of infant mortality data in the province of British Columbia, Canada. The goals of the analysis are to produce more reliable small-area estimates of mortality rates, assess spatial patterns over time, and examine risk trends at both global (provincial) and local (local health area) levels.

British Columbia↗