Demographic characteristics and opportunistic diseases associated with attrition during preparation for antiretroviral therapy in primary health centres in Kibera, Kenya.

Using routine data from HIV-positive adult patients eligible for antiretroviral therapy (ART), we report on routinely collected demographic characteristics and opportunistic diseases associated with pre-ART attrition (deaths and loss to follow-up). Among 2471 ART eligible patients, enrolled between January 2005 and November 2008, 446 (18%) were lost to attrition pre-ART. Adjusted risk factors significantly associated with pre-ART attrition included age <35 years (Odds Ratio, OR 1.4, 95% Confidence Interval, CI 1.1-1.8), severe malnutrition (OR 1.5, 95% CI 1.1-2.0), active pulmonary tuberculosis (OR 1.6, 95% CI 1.1-2.4), severe bacterial infections including severe bacterial pneumonia (OR 1.9, 95% CI 1.2-2.8) and prolonged unexplained fever (>1 month), (OR 2.6, 95% CI 1.3-5.2). This study highlights a number of clinical markers associated with pre-ART attrition that could serve as 'pointers' or screening tools to identify patients who merit fast-tracking onto ART and/or closer clinical attention and follow-up

Lie algebra solution of population models based on time-inhomogeneous Markov chains

Many natural populations are well modelled through time-inhomogeneous stochastic processes. Such processes have been analysed in the physical sciences using a method based on Lie algebras, but this methodology is not widely used for models with ecological, medical and social applications. This paper presents the Lie algebraic method, and applies it to three biologically well motivated examples. The result of this is a solution form that is often highly computationally advantageous.Comment: 10 pages; 1 figure; 2 tables. To appear in Applied Probabilit

Catching Super Massive Black Hole Binaries Without a Net

The gravitational wave signals from coalescing Supermassive Black Hole Binaries are prime targets for the Laser Interferometer Space Antenna (LISA). With optimal data processing techniques, the LISA observatory should be able to detect black hole mergers anywhere in the Universe. The challenge is to find ways to dig the signals out of a combination of instrument noise and the large foreground from stellar mass binaries in our own galaxy. The standard procedure of matched filtering against a grid of templates can be computationally prohibitive, especially when the black holes are spinning or the mass ratio is large. Here we develop an alternative approach based on Metropolis-Hastings sampling and simulated annealing that is orders of magnitude cheaper than a grid search. We demonstrate our approach on simulated LISA data streams that contain the signals from binary systems of Schwarzschild Black Holes, embedded in instrument noise and a foreground containing 26 million galactic binaries. The search algorithm is able to accurately recover the 9 parameters that describe the black hole binary without first having to remove any of the bright foreground sources, even when the black hole system has low signal-to-noise.Comment: 4 pages, 3 figures, Refined search algorithm, added low SNR exampl

Estimating population cardinal health state valuation models from individual ordinal (rank) health state preference data

Ranking exercises have routinely been used as warm-up exercises within health state valuation surveys. Very little use has been made of the information obtained in this process. Instead, research has focussed upon the analysis of health state valuation data obtained using the visual analogue scale, standard gamble and time trade off methods. Thurstone’s law of comparative judgement postulates a stable relationship between ordinal and cardinal preferences, based upon the information provided by pairwise choices. McFadden proposed that this relationship could be modelled by estimating conditional logistic regression models where alternatives had been ranked. In this paper we report the estimation of such models for the Health Utilities Index Mark 2 and the SF-6D. The results are compared to the conventional regression models estimated from standard gamble data, and to the observed mean standard gamble health state valuations. For both the HUI2 and the SF-6D, the models estimated using rank data are broadly comparable to the models estimated on standard gamble data and the predictive performance of these models is close to that of the standard gamble models. Our research indicates that rank data has the potential to provide useful insights into community health state preferences. However, important questions remain

Markov Chain Monte Carlo Method without Detailed Balance

We present a specific algorithm that generally satisfies the balance condition without imposing the detailed balance in the Markov chain Monte Carlo. In our algorithm, the average rejection rate is minimized, and even reduced to zero in many relevant cases. The absence of the detailed balance also introduces a net stochastic flow in a configuration space, which further boosts up the convergence. We demonstrate that the autocorrelation time of the Potts model becomes more than 6 times shorter than that by the conventional Metropolis algorithm. Based on the same concept, a bounce-free worm algorithm for generic quantum spin models is formulated as well.Comment: 5 pages, 5 figure

Bayesian Overdispersed Poisson Model and the Bornhuetter-Ferguson Claim Reserving Method

We consider the Bayesian overdispersed Poisson (ODP) model for claims reserving in general insurance. We choose two different types of prior distributions for the parameters and then study the different Bayesian predictors. This study leads, on the one hand, to the classical chain ladder predictor and, on the other hand, to Bornhuetter-Ferguson predictors. We highlight (either analytically or numerically) how these predictors are obtained and how their prediction uncertainty can be determined

Using Markov chain Monte Carlo methods for estimating parameters with gravitational radiation data

We present a Bayesian approach to the problem of determining parameters for coalescing binary systems observed with laser interferometric detectors. By applying a Markov Chain Monte Carlo (MCMC) algorithm, specifically the Gibbs sampler, we demonstrate the potential that MCMC techniques may hold for the computation of posterior distributions of parameters of the binary system that created the gravity radiation signal. We describe the use of the Gibbs sampler method, and present examples whereby signals are detected and analyzed from within noisy data.Comment: 21 pages, 10 figure

Estimating population cardinal health state valuation models from individual ordinal (rank) health state preference data

Ranking exercises have routinely been used as warm-up exercises within health state valuation surveys. Very little use has been made of the information obtained in this process. Instead, research has focussed upon the analysis of health state valuation data obtained using the visual analogue scale, standard gamble and time trade off methods. Thurstone’s law of comparative judgement postulates a stable relationship between ordinal and cardinal preferences, based upon the information provided by pairwise choices. McFadden proposed that this relationship could be modelled by estimating conditional logistic regression models where alternatives had been ranked. In this paper we report the estimation of such models for the Health Utilities Index Mark 2 and the SF-6D. The results are compared to the conventional regression models estimated from standard gamble data, and to the observed mean standard gamble health state valuations. For both the HUI2 and the SF-6D, the models estimated using rank data are broadly comparable to the models estimated on standard gamble data and the predictive performance of these models is close to that of the standard gamble models. Our research indicates that rank data has the potential to provide useful insights into community health state preferences. However, important questions remain

Metropolis Sampling

Monte Carlo (MC) sampling methods are widely applied in Bayesian inference, system simulation and optimization problems. The Markov Chain Monte Carlo (MCMC) algorithms are a well-known class of MC methods which generate a Markov chain with the desired invariant distribution. In this document, we focus on the Metropolis-Hastings (MH) sampler, which can be considered as the atom of the MCMC techniques, introducing the basic notions and different properties. We describe in details all the elements involved in the MH algorithm and the most relevant variants. Several improvements and recent extensions proposed in the literature are also briefly discussed, providing a quick but exhaustive overview of the current Metropolis-based sampling's world.Comment: Wiley StatsRef-Statistics Reference Online, 201