Estimating the cost-effectiveness of modern screening mammography programmes

Context: The debate about the balance of risks and benefits of mammography screening has prompted a comprehensive review by an independent panel in the UK. However, the panel's remit did not cover the important economic dimension of breast cancer screening. Methods: The life histories of two cohort studies of 50-year-old women, who would be eligible within the National Health Service (NHS) breast screening programme (NHSBSP), were simulated over 35 years, using a life table approach. One cohort participant was offered screening at age 50 and triennially thereafter until age 70, assuming 75% attendance, while the other received no screening. Based on the findings from the panel's report, the cost-effectiveness of the NHSBSP was assessed for various scenarios of screening effect on breast cancer incidence (base case scenario: screening advances diagnosis by 5 years; 10% incidence reduction after screening stops).

One kind of lawlessness: estimating the welfare cost of Somali piracy

This paper estimates the effect of piracy attacks on shipping costs using a unique data set on shipping contracts in the dry bulk market. We look at shipping routes whose shortest path exposes them to piracy attacks and find that the increase in attacks in 2008 lead to around a ten percent increase in shipping costs. We use this estimate to get a sense of the welfare loss imposed by piracy. Our intermediate estimate suggests that the creation of $120 million of revenue for pirates in the Somalia area led to a welfare loss of over $1.5 billion.

Estimating national-level syringe availability to injecting drug users and injection coverage: Switzerland, 1996-2006.

BACKGROUND: Measuring syringe availability and coverage is essential in the assessment of HIV/AIDS risk reduction policies. Estimates of syringe availability and coverage were produced for the years 1996 and 2006, based on all relevant available national-level aggregated data from published sources. METHODS: We defined availability as the total monthly number of syringes provided by harm reduction system divided by the estimated number of injecting drug users (IDU), and defined coverage as the proportion of injections performed with a new syringe, at national level (total supply over total demand). Estimates of supply of syringes were derived from the national monitoring system, including needle and syringe programmes (NSP), pharmacies, and medically prescribed heroin programmes. Estimates of syringe demand were based on the number of injections performed by IDU derived from surveys of low threshold facilities for drug users (LTF) with NSP combined with the number of IDU. This number was estimated by two methods combining estimates of heroin users (multiple estimation method) and (a) the number of IDU in methadone treatment (MT) (non-injectors) or (b) the proportion of injectors amongst LTF attendees. Central estimates and ranges were obtained for availability and coverage. RESULTS: The estimated number of IDU decreased markedly according to both methods. The MT-based method (from 14,818 to 4809) showed a much greater decrease and smaller size of the IDU population compared to the LTF-based method (from 24,510 to 12,320). Availability and coverage estimates are higher with the MT-based method. For 1996, central estimates of syringe availability were 30.5 and 18.4 per IDU per month; for 2006, they were 76.5 and 29.9. There were 4 central estimates of coverage. For 1996 they ranged from 24.3% to 43.3%, and for 2006, from 50.5% to 134.3%. CONCLUSION: Although 2006 estimates overlap 1996 estimates, the results suggest a shift to improved syringe availability and coverage over time.

Prediction of fat-free mass using bioelectrical impedance analysis in young adults from five populations of African origin.

Background/objectives:Bioelectrical impedance analysis (BIA) is used in population and clinical studies as a technique for estimating body composition. Because of significant under-representation in existing literature, we sought to develop and validate predictive equation(s) for BIA for studies in populations of African origin.Subjects/methods:Among five cohorts of the Modeling the Epidemiologic Transition Study, height, weight, waist circumference and body composition, using isotope dilution, were measured in 362 adults, ages 25-45 with mean body mass indexes ranging from 24 to 32. BIA measures of resistance and reactance were measured using tetrapolar placement of electrodes and the same model of analyzer across sites (BIA 101Q, RJL Systems). Multiple linear regression analysis was used to develop equations for predicting fat-free mass (FFM), as measured by isotope dilution; covariates included sex, age, waist, reactance and height(2)/resistance, along with dummy variables for each site. Developed equations were then tested in a validation sample; FFM predicted by previously published equations were tested in the total sample.Results:A site-combined equation and site-specific equations were developed. The mean differences between FFM (reference) and FFM predicted by the study-derived equations were between 0.4 and 0.6âeuro0/00kg (that is, 1% difference between the actual and predicted FFM), and the measured and predicted values were highly correlated. The site-combined equation performed slightly better than the site-specific equations and the previously published equations.Conclusions:Relatively small differences exist between BIA equations to estimate FFM, whether study-derived or published equations, although the site-combined equation performed slightly better than others. The study-derived equations provide an important tool for research in these understudied populations.

Estimating arterial blood gases with the EABC® System (“Earlobe Arterialized Blood Collector”) in critically ill patients. A validation study.

The Earlobe Arterialized Blood Collector® is a minimally invasive system able to perform arterialized capillary blood gas analysis from the earlobe (EL). A prospective validation study was performed in 55 critical ill patients. Sampling failure rate was high (53.6%). Risk factors were age > 65 years, diabetes, vasoactive drug therapy and noradrenaline (NA) doses above 0.22 μg / kg / min. Multivariate analysis showed age > 65 years was the only factor independently associated with failure. Concordance analysis with arterial blood gases and Bland-Altman agreement evaluation were insufficient for validating the new system for all gasometrical variables.

Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Hyperbolic reaction-diffusion equations for a forest fire model

Forest fire models have been widely studied from the context of self-organized criticality and from the ecological properties of the forest and combustion. On the other hand, reaction-diffusion equations have interesting applications in biology and physics. We propose here a model for fire propagation in a forest by using hyperbolic reaction-diffusion equations. The dynamical and thermodynamical aspects of the model are analyzed in detail

Speed of wave-front solutions to hyperbolic reaction-diffusion equations

The asymptotic speed problem of front solutions to hyperbolic reaction-diffusion (HRD) equations is studied in detail. We perform linear and variational analyses to obtain bounds for the speed. In contrast to what has been done in previous work, here we derive upper bounds in addition to lower ones in such a way that we can obtain improved bounds. For some functions it is possible to determine the speed without any uncertainty. This is also achieved for some systems of HRD (i.e., time-delayed Lotka-Volterra) equations that take into account the interaction among different species. An analytical analysis is performed for several systems of biological interest, and we find good agreement with the results of numerical simulations as well as with available observations for a system discussed recently

S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Estimation of inelastic seismic material properties of a surgical sea-bed from multi-component marine seismic data

To date, state-of-the-art seismic material parameter estimates from multi-component sea-bed seismic data are based on the assumption that the sea-bed consists of a fully elastic half-space. In reality, however, the shallow sea-bed generally consists of soft, unconsolidated sediments that are characterized by strong to very strong seismic attenuation. To explore the potential implications, we apply a state-of-the-art elastic decomposition algorithm to synthetic data for a range of canonical sea-bed models consisting of a viscoelastic half-space of varying attenuation. We find that in the presence of strong seismic attenuation, as quantified by Q-values of 10 or less, significant errors arise in the conventional elastic estimation of seismic properties. Tests on synthetic data indicate that these errors can be largely avoided by accounting for the inherent attenuation of the seafloor when estimating the seismic parameters. This can be achieved by replacing the real-valued expressions for the elastic moduli in the governing equations in the parameter estimation by their complex-valued viscoelastic equivalents. The practical application of our parameter procedure yields realistic estimates of the elastic seismic material properties of the shallow sea-bed, while the corresponding Q-estimates seem to be biased towards too low values, particularly for S-waves. Given that the estimation of inelastic material parameters is notoriously difficult, particularly in the immediate vicinity of the sea-bed, this is expected to be of interest and importance for civil and ocean engineering purposes.

A comparison of random sequence reads versus 16S rDNA sequences for estimating the biodiversity of a metagenomic library

The construction of metagenomic libraries has permitted the study of microorganisms resistant to isolation and the analysis of 16S rDNA sequences has been used for over two decades to examine bacterial biodiversity. Here, we show that the analysis of random sequence reads (RSRs) instead of 16S is a suitable shortcut to estimate the biodiversity of a bacterial community from metagenomic libraries. We generated 10,010 RSRs from a metagenomic library of microorganisms found in human faecal samples. Then searched them using the program BLASTN against a prokaryotic sequence database to assign a taxon to each RSR. The results were compared with those obtained by screening and analysing the clones containing 16S rDNA sequences in the whole library. We found that the biodiversity observed by RSR analysis is consistent with that obtained by 16S rDNA. We also show that RSRs are suitable to compare the biodiversity between different metagenomic libraries. RSRs can thus provide a good estimate of the biodiversity of a metagenomic library and, as an alternative to 16S, this approach is both faster and cheaper.

Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.

Spirometry reference equations for central European populations from school age to old age.

BACKGROUND: Spirometry reference values are important for the interpretation of spirometry results. Reference values should be updated regularly, derived from a population as similar to the population for which they are to be used and span across all ages. Such spirometry reference equations are currently lacking for central European populations. OBJECTIVE: To develop spirometry reference equations for central European populations between 8 and 90 years of age. MATERIALS: We used data collected between January 1993 and December 2010 from a central European population. The data was modelled using "Generalized Additive Models for Location, Scale and Shape" (GAMLSS). RESULTS: The spirometry reference equations were derived from 118'891 individuals consisting of 60'624 (51%) females and 58'267 (49%) males. Altogether, there were 18'211 (15.3%) children under the age of 18 years. CONCLUSION: We developed spirometry reference equations for a central European population between 8 and 90 years of age that can be implemented in a wide range of clinical settings.

A subsampling approach to estimating the distribution of diversing statistics with application to assessing financial market risks

In this paper we propose a subsampling estimator for the distribution ofstatistics diverging at either known rates when the underlying timeseries in strictly stationary abd strong mixing. Based on our results weprovide a detailed discussion how to estimate extreme order statisticswith dependent data and present two applications to assessing financialmarket risk. Our method performs well in estimating Value at Risk andprovides a superior alternative to Hill's estimator in operationalizingSafety First portofolio selection.

Estimating the effect of transitory economic shocks on civil conflict

This note tries to clarify some remaining issues in the debate on the effect of income shocks oncivil conflict. Section 1 discusses the discrepant findings on the effect of rainfall shocks oncivil conflict in Miguel and Satyanath (2010, 2011) and Ciccone (2011). Section 2 develops aninstrumental variables approach to estimate the effect of transitory (rainfall-driven) incomeshocks on civil conflict and contrasts the conclusions with those of Miguel, Satyanath, andSergenti (2004) and Miguel and Satyanath (2010, 2011). Throughout, the note uses the data ofMiguel, Satyanath, and Sergenti to focus on the methodological issues at the core of the debate(for results using the latest data see Ciccone, 2011).