Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

High nuclear genetic diversity, high levels of outcrossing and low differentiation among remnant populations of Quercus petraea at the margin of its range in Ireland

Background and Aims Quercus petraea colonized Ireland after the last glaciation from refugia on mainland Europe. Deforestation. however. beginning in Neolithic times, has resulted in small, scattered forest fragments, now covering less than 12 000 ha. Methods Plastid (three fragments) and microsatellite variation (13 loci) were characterized in seven Irish populations sampled along a north-south gradient. Using Bayesian approaches and Wright's F-statistics, the effects of colonization and fragmentation on the genetic structure and mating patterns of extant oak populations were investigated. Key-Results All Populations possessed cytotypes common to the Iberian Peninsula. Despite the distance from the refugial core and the extensive deforestation in Ireland, nuclear genetic variation was high and comparable to mainland Europe. Low population differentiation was observed within Ireland and populations showed no evidence for isolation by distance. As expected of a marker with an effective Population size of one-quarter relative to the nuclear genome, plastid variation indicated higher differentiation. Individual inbreeding coefficients indicated high levels of outcrossing. Conclusions Consistent with a large effective Population size in the historical migrant gene pool and/or with high gene flow among populations, high within-population diversity and low population differentiation was observred within Ireland. It is proposed that native Q. petraea populations in Ireland share a common phylogeographic history and that the present genetic structure does not reflect founder effects. (C) 2004 Annals of Botany Company.

Analysis of population pharmacokinetic data using NONMEM and WinBUGS

The aim of this report is to describe the use of WinBUGS for two datasets that arise from typical population pharmacokinetic studies. The first dataset relates to gentamicin concentration-time data that arose as part of routine clinical care of 55 neonates. The second dataset incorporated data from 96 patients receiving enoxaparin. Both datasets were originally analyzed by using NONMEM. In the first instance, although NONMEM provided reasonable estimates of the fixed effects parameters it was unable to provide satisfactory estimates of the between-subject variance. In the second instance, the use of NONMEM resulted in the development of a successful model, albeit with limited available information on the between-subject variability of the pharmacokinetic parameters. WinBUGS was used to develop a model for both of these datasets. Model comparison for the enoxaparin dataset was performed by using the posterior distribution of the log-likelihood and a posterior predictive check. The use of WinBUGS supported the same structural models tried in NONMEM. For the gentamicin dataset a one-compartment model with intravenous infusion was developed, and the population parameters including the full between-subject variance-covariance matrix were available. Analysis of the enoxaparin dataset supported a two compartment model as superior to the one-compartment model, based on the posterior predictive check. Again, the full between-subject variance-covariance matrix parameters were available. Fully Bayesian approaches using MCMC methods, via WinBUGS, can offer added value for analysis of population pharmacokinetic data.

Bayesian feedback versus Markovian feedback in a two-level atom

We compare two different approaches to the control of the dynamics of a continuously monitored open quantum system. The first is Markovian feedback, as introduced in quantum optics by Wiseman and Milburn [Phys. Rev. Lett. 70, 548 (1993)]. The second is feedback based on an estimate of the system state, developed recently by Doherty and Jacobs [Phys. Rev. A 60, 2700 (1999)]. Here we choose to call it, for brevity, Bayesian feedback. For systems with nonlinear dynamics, we expect these two methods of feedback control to give markedly different results. The simplest possible nonlinear system is a driven and damped two-level atom, so we choose this as our model system. The monitoring is taken to be homodyne detection of the atomic fluorescence, and the control is by modulating the driving. The aim of the feedback in both cases is to stabilize the internal state of the atom as close as possible to an arbitrarily chosen pure state, in the presence of inefficient detection and other forms of decoherence. Our results (obtained without recourse to stochastic simulations) prove that Bayesian feedback is never inferior, and is usually superior, to Markovian feedback. However, it would be far more difficult to implement than Markovian feedback and it loses its superiority when obvious simplifying approximations are made. It is thus not clear which form of feedback would be better in the face of inevitable experimental imperfections.

Analysis of Melanoma Onset: Assessing Familial Aggregation by Using Estimating Equations and Fitting Variance Components via Bayesian Random Effects Models

We investigate whether relative contributions of genetic and shared environmental factors are associated with an increased risk in melanoma. Data from the Queensland Familial Melanoma Project comprising 15,907 subjects arising from 1912 families were analyzed to estimate the additive genetic, common and unique environmental contributions to variation in the age at onset of melanoma. Two complementary approaches for analyzing correlated time-to-onset family data were considered: the generalized estimating equations (GEE) method in which one can estimate relationship-specific dependence simultaneously with regression coefficients that describe the average population response to changing covariates; and a subject-specific Bayesian mixed model in which heterogeneity in regression parameters is explicitly modeled and the different components of variation may be estimated directly. The proportional hazards and Weibull models were utilized, as both produce natural frameworks for estimating relative risks while adjusting for simultaneous effects of other covariates. A simple Markov Chain Monte Carlo method for covariate imputation of missing data was used and the actual implementation of the Bayesian model was based on Gibbs sampling using the free ware package BUGS. In addition, we also used a Bayesian model to investigate the relative contribution of genetic and environmental effects on the expression of naevi and freckles, which are known risk factors for melanoma.

A two-phase approach for detecting recombination in nucleotide sequences

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. Delineating recombination events is important in the study of molecular evolution, as inference of such events provides a clearer picture of the phylogenetic relationships among different gene sequences or genomes. Nevertheless, detecting recombination events can be a daunting task, as the performance of different recombination-detecting approaches can vary, depending on evolutionary events that take place after recombination. We recently evaluated the effects of post-recombination events on the prediction accuracy of recombination-detecting approaches using simulated nucleotide sequence data. The main conclusion, supported by other studies, is that one should not depend on a single method when searching for recombination events. In this paper, we introduce a two-phase strategy, applying three statistical measures to detect the occurrence of recombination events, and a Bayesian phylogenetic approach in delineating breakpoints of such events in nucleotide sequences. We evaluate the performance of these approaches using simulated data, and demonstrate the applicability of this strategy to empirical data. The two-phase strategy proves to be time-efficient when applied to large datasets, and yields high-confidence results.

Management Approaches to Merger Evoked Cultural Change and Acculturation Outcomes

In a study of merger-evoked cultural change in three organizations, quantitative and qualitative data were collected from individuals at all employment levels in both merger partners within each organization. Results were that most individuals perceived that the merger had impacted significantly on them personally. There was, however, a perceived lack of congruence between the organizational cultures of merging partners, resulting in cultural clashes and significant changes to the organizations' organizational cultures. More specifically, outcomes for both individuals and the subsequent acculturation following the mergers were related to the approach adopted to manage the merger process: incremental, immediate, or indifferent.

A two-phase strategy for detecting recombination in nucleotide sequences

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. Delineating recombination events is important in the study of molecular evolution, as inference of such events provides a clearer picture of the phylogenetic relationships among different gene sequences or genomes. Nevertheless, detecting recombination events can be a daunting task, as the performance of different recombination-detecting approaches can vary, depending on evolutionary events that take place after recombination. We previously evaluated the effects of post-recombination events on the prediction accuracy of recombination-detecting approaches using simulated nucleotide sequence data. The main conclusion, supported by other studies, is that one should not depend on a single method when searching for recombination events. In this paper, we introduce a two-phase strategy, applying three statistical measures to detect the occurrence of recombination events, and a Bayesian phylogenetic approach to delineate breakpoints of such events in nucleotide sequences. We evaluate the performance of these approaches using simulated data, and demonstrate the applicability of this strategy to empirical data. The two-phase strategy proves to be time-efficient when applied to large datasets, and yields high-confidence results.

The Swiss heroin trials: testing alternative approaches

Over half a million heroin misusers receive oral methadone maintenance treatment world-wide1 but the maintenance prescription of injectable opioid drugs, like heroin, remains controversial. In 1992 Switzerland began a large scale evaluation of heroin and other injectable opiate prescribing that eventually involved 1035 misusers. 2 3 The results of the evaluation have recently been reported.4 These show that it was feasible to provide heroin by intravenous injection at a clinic, up to three times a day, for seven days a week. This was done while maintaining good drug control, good order, client safety, and staff morale. Patients were stabilised on 500 to 600 mg heroin daily without evidence of increasing tolerance. Retention in treatment was 89% at six months and 69% at 18 months.4 The self reported use of non-prescribed heroin fell signifianctly, but other drug use was minimally affected. The death rate was 1% per year, and there were no deaths from overdose among participants . . . [Full text of this article]

Challenges and approaches to estimating mortality attributable to the use of selected illicit drugs

A number of unique challenges are faced when attempting to estimate mortality attributable to illicit drugs. The hidden nature of illicit drug use creates difficulties in quantifying the prevalence of such use; identifying adverse health effects associated with exposure, and calculating the risk of these effects. The use of cohort studies of drug users allows the identification of causes of mortality associated with drug use and the determination of the risk of these causes. This risk estimate can then be used in conjunction with estimates of the prevalence of drug use to, extrapolate the burden of mortality. We identify a number of such studies and present some solutions to the major challenges faced when attempting to estimate the global burden of mortality attributable to illicit drug use. Copyright (C) 2001 S. Karger AG, Basel.

A context noise model of episodic word recognition

Item noise models of recognition assert that interference at retrieval is generated by the words from the study list. Context noise models of recognition assert that interference at retrieval is generated by the contexts in which the test word has appeared. The authors introduce the bind cue decide model of episodic memory, a Bayesian context noise model, and demonstrate how it can account for data from the item noise and dual-processing approaches to recognition memory. From the item noise perspective, list strength and list length effects, the mirror effect for word frequency and concreteness, and the effects of the similarity of other words in a list are considered. From the dual-processing perspective, process dissociation data on the effects of length, temporal separation of lists, strength, and diagnosticity of context are examined. The authors conclude that the context noise approach to recognition is a viable alternative to existing approaches.