Cost-effectiveness of pharmacotherapies for nicotine dependence in primary care settings: a multinational comparison.

OBJECTIVE: To estimate the incremental cost-effectiveness of the first-line pharmacotherapies (nicotine gum, patch, spray, inhaler, and bupropion) for smoking cessation across six Western countries-Canada, France, Spain, Switzerland, the United States, and the United Kingdom. DESIGN AND STUDY POPULATION: A Markov-chain cohort model to simulate two cohorts of smokers: (1) a reference cohort given brief cessation counselling by a general practitioner (GP); (2) a treatment cohort given counselling plus pharmacotherapy. Effectiveness expressed as odds ratios for quitting associated with pharmacotherapies. Costs based on the additional physician time required and retail prices of the medications. INTERVENTIONS: Addition of each first-line pharmacotherapy to GP cessation counselling. MAIN OUTCOME MEASURES: Cost per life-year saved associated with pharmacotherapies. RESULTS: The cost per life-year saved for counselling only ranged from US190 dollars in Spain to 773 dollars in the UK for men, and from 288 dollars in Spain to 1168 dollars in the UK for women. The incremental cost per life-year saved for gum ranged from 2230 dollars for men in Spain to 7643 dollars for women in the US; for patch from 1758 dollars for men in Spain to 5131 dollars for women in the UK; for spray from 1935 dollars for men in Spain to 7969 dollars for women in the US; for inhaler from 3480 dollars for men in Switzerland to 8700 dollars for women in France; and for bupropion from 792 dollars for men in Canada to 2922 dollars for women in the US. In sensitivity analysis, changes in discount rate, treatment effectiveness, and natural quit rate had the strongest influences on cost-effectiveness. CONCLUSIONS: The cost-effectiveness of the pharmacotherapies varied significantly across the six study countries, however, in each case, the results would be considered favourable as compared to other common preventive pharmacotherapies.

Evaluation of adding a second marker to overcome Staphylococcus aureus spa typing homoplasies.

The utility of sequencing a second highly variable locus in addition to the spa gene (e.g., double-locus sequence typing [DLST]) was investigated to overcome limitations of a Staphylococcus aureus single-locus typing method. Although adding a second locus seemed to increase discriminatory power, it was not sufficient to definitively infer evolutionary relationships within a single multilocus sequence type (ST-5).

Panel index VAR models: Specification, Estimation, Testing and Leading Indicators

This paper proposes a method to conduct inference in panel VAR models with cross unit interdependencies and time variations in the coefficients. The approach can be used to obtain multi-unit forecasts and leading indicators and to conduct policy analysis in a multiunit setups. The framework of analysis is Bayesian and MCMC methods are used to estimate the posterior distribution of the features of interest. The model is reparametrized to resemble an observable index model and specification searches are discussed. As an example, we construct leading indicators for inflation and GDP growth in the Euro area using G-7 information.

The Effect of Four-Lane to Three-Lane Conversion on the Number of Crashes and Crash Rates in Iowa Roads, June 2005

We analyze crash data collected by the Iowa Department of Transportation using Bayesian methods. The data set includes monthly crash numbers, estimated monthly traffic volumes, site length and other information collected at 30 paired sites in Iowa over more than 20 years during which an intervention experiment was set up. The intervention consisted in transforming 15 undivided road segments from four-lane to three lanes, while an additional 15 segments, thought to be comparable in terms of traffic safety-related characteristics were not converted. The main objective of this work is to find out whether the intervention reduces the number of crashes and the crash rates at the treated sites. We fitted a hierarchical Poisson regression model with a change-point to the number of monthly crashes per mile at each of the sites. Explanatory variables in the model included estimated monthly traffic volume, time, an indicator for intervention reflecting whether the site was a “treatment” or a “control” site, and various interactions. We accounted for seasonal effects in the number of crashes at a site by including smooth trigonometric functions with three different periods to reflect the four seasons of the year. A change-point at the month and year in which the intervention was completed for treated sites was also included. The number of crashes at a site can be thought to follow a Poisson distribution. To estimate the association between crashes and the explanatory variables, we used a log link function and added a random effect to account for overdispersion and for autocorrelation among observations obtained at the same site. We used proper but non-informative priors for all parameters in the model, and carried out all calculations using Markov chain Monte Carlo methods implemented in WinBUGS. We evaluated the effect of the four to three-lane conversion by comparing the expected number of crashes per year per mile during the years preceding the conversion and following the conversion for treatment and control sites. We estimated this difference using the observed traffic volumes at each site and also on a per 100,000,000 vehicles. We also conducted a prospective analysis to forecast the expected number of crashes per mile at each site in the study one year, three years and five years following the four to three-lane conversion. Posterior predictive distributions of the number of crashes, the crash rate and the percent reduction in crashes per mile were obtained for each site for the months of January and June one, three and five years after completion of the intervention. The model appears to fit the data well. We found that in most sites, the intervention was effective and reduced the number of crashes. Overall, and for the observed traffic volumes, the reduction in the expected number of crashes per year and mile at converted sites was 32.3% (31.4% to 33.5% with 95% probability) while at the control sites, the reduction was estimated to be 7.1% (5.7% to 8.2% with 95% probability). When the reduction in the expected number of crashes per year, mile and 100,000,000 AADT was computed, the estimates were 44.3% (43.9% to 44.6%) and 25.5% (24.6% to 26.0%) for converted and control sites, respectively. In both cases, the difference in the percent reduction in the expected number of crashes during the years following the conversion was significantly larger at converted sites than at control sites, even though the number of crashes appears to decline over time at all sites. Results indicate that the reduction in the expected number of sites per mile has a steeper negative slope at converted than at control sites. Consistent with this, the forecasted reduction in the number of crashes per year and mile during the years after completion of the conversion at converted sites is more pronounced than at control sites. Seasonal effects on the number of crashes have been well-documented. In this dataset, we found that, as expected, the expected number of monthly crashes per mile tends to be higher during winter months than during the rest of the year. Perhaps more interestingly, we found that there is an interaction between the four to three-lane conversion and season; the reduction in the number of crashes appears to be more pronounced during months, when the weather is nice than during other times of the year, even though a reduction was estimated for the entire year. Thus, it appears that the four to three-lane conversion, while effective year-round, is particularly effective in reducing the expected number of crashes in nice weather.

Estimating multi-country VAR models

This paper describes a methodology to estimate the coefficients, to test specification hypothesesand to conduct policy exercises in multi-country VAR models with cross unit interdependencies, unit specific dynamics and time variations in the coefficients. The framework of analysis is Bayesian: a prior flexibly reduces the dimensionality of the model and puts structure on the time variations; MCMC methods are used to obtain posterior distributions; and marginal likelihoods to check the fit of various specifications. Impulse responses and conditional forecasts are obtained with the output of MCMC routine. The transmission of certain shocks across countries is analyzed.

Characterization of a new clinical yeast species, Candida tunisiensis sp. nov., isolated from a strain collection from Tunisian hospitals.

From a collection of yeast isolates isolated from patients in Tunisian hospitals between September 2006 and July 2010, the yeast strain JEY63 (CBS 12513), isolated from a 50-year-old male that suffered from oral thrush, could not be identified to the species level using conventional methods used in clinical laboratories. These methods include matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS), germ tube formation, and the use of CHROMagar Candida and metabolic galleries. Sequence analysis of the nuclear rRNA (18S rRNA, 5.8S rRNA, and 26S rRNA) and internal transcribed spacer regions (ITS1 and ITS2) indicated that the ribosomal DNA sequences of this species were not yet reported. Multiple gene phylogenic analyses suggested that this isolate clustered at the base of the Dipodascaceae (Saccharomycetales, Saccharomycetes, and Ascomycota). JEY63 was named Candida tunisiensis sp. nov. according to several phenotypic criteria and its geographical origin. C. tunisiensis was able to grow at 42°C and does not form chlamydospores and hyphae but could grow as yeast and pseudohyphal forms. C. tunisiensis exhibited most probably a haploid genome with an estimated size of 10 Mb on at least three chromosomes. Using European Committee for Antimicrobial Susceptibility Testing (EUCAST) and Clinical and Laboratory Standards Institute (CLSI) Candida albicans susceptibility breakpoints as a reference, C. tunisiensis was resistant to fluconazole (MIC = 8 μg/ml), voriconazole (MIC = 0.5 μg/ml), itraconazole (MIC = 16 μg/ml), and amphotericin B (MIC = 4 μg/ml) but still susceptible to posaconazole (MIC = 0.008 μg/ml) and caspofungin (MIC = 0.5 μg/ml). In conclusion, MALDI-TOF MS permitted the early selection of an unusual isolate, which was still unreported in molecular databases but could not be unambiguously classified based on phylogenetic approaches.

Stochastic models and numerical algorithms for a class of regulatory gene networks.

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.

The cost effectiveness of pharmacological smoking cessation therapies in developing countries: a case study in the Seychelles

OBJECTIVE: To examine the incremental cost effectiveness of the five first line pharmacological smoking cessation therapies in the Seychelles and other developing countries. DESIGN: A Markov chain cohort simulation. SUBJECTS: Two simulated cohorts of smokers: (1) a reference cohort given physician counselling only; (2) a treatment cohort given counselling plus cessation therapy. INTERVENTION: Addition of each of the five pharmacological cessation therapies to physician provided smoking cessation counselling. MAIN OUTCOME MEASURES: Cost per life-year saved (LYS) associated with the five pharmacotherapies. Effectiveness expressed as odds ratios for quitting associated with pharmacotherapies. Costs based on the additional physician time required and retail prices of the medications. RESULTS: Based on prices for currently available generic medications on the global market, the incremental cost per LYS for a 45 year old in the Seychelles was 599 US dollars for gum and 227 dollars for bupropion. Assuming US treatment prices as a conservative estimate, the incremental cost per LYS was significantly higher, though still favourable in comparison to other common medical interventions: 3712 dollars for nicotine gum, 1982 dollars for nicotine patch, 4597 dollars for nicotine spray, 4291 dollars for nicotine inhaler, and 1324 dollars for bupropion. Cost per LYS increased significantly upon application of higher discount rates, which may be used to reflect relatively high opportunity costs for health expenditures in developing countries with highly constrained resources and high overall mortality. CONCLUSION: Pharmacological cessation therapy can be highly cost effective as compared to other common medical interventions in low mortality, middle income countries, particularly if medications can be procured at low prices.

Cell viability and noninvasive in vivo MRI tracking of 3D cell encapsulating self-assembled microcontainers.

Several molecular therapies require the implantation of cells that secrete biotherapeutic molecules and imaging the location and microenvironment of the cellular implant to ascertain its function. We demonstrate noninvasive in vivo magnetic resonance imaging (MRI) of self-assembled microcontainers that are capable of cell encapsulation. Negative contrast was obtained to discern the microcontainer with MRI; positive contrast was obtained in the complete absence of background signal. MRI on a clinical scanner highlights the translational nature of this research. The microcontainers were loaded with cells that were dispersed in an extracellular matrix, and implanted both subcutaneously and in human tumor xenografts in SCID mice. MRI was performed on the implants, and microcontainers retrieved postimplantation showed cell viability both within and proximal to the implant. The microcontainers are characterized by their small size, three dimensionality, controlled porosity, ease of parallel fabrication, chemical and mechanical stability, and noninvasive traceability in vivo.

RDM1, a novel RNA recognition motif (RRM)-containing protein involved in the cell response to cisplatin in vertebrates.

A variety of cellular proteins has the ability to recognize DNA lesions induced by the anti-cancer drug cisplatin, with diverse consequences on their repair and on the therapeutic effectiveness of this drug. We report a novel gene involved in the cell response to cisplatin in vertebrates. The RDM1 gene (for RAD52 Motif 1) was identified while searching databases for sequences showing similarities to RAD52, a protein involved in homologous recombination and DNA double-strand break repair. Ablation of RDM1 in the chicken B cell line DT40 led to a more than 3-fold increase in sensitivity to cisplatin. However, RDM1-/- cells were not hypersensitive to DNA damages caused by ionizing radiation, UV irradiation, or the alkylating agent methylmethane sulfonate. The RDM1 protein displays a nucleic acid binding domain of the RNA recognition motif (RRM) type. By using gel-shift assays and electron microscopy, we show that purified, recombinant chicken RDM1 protein interacts with single-stranded DNA as well as double-stranded DNA, on which it assembles filament-like structures. Notably, RDM1 recognizes DNA distortions induced by cisplatin-DNA adducts in vitro. Finally, human RDM1 transcripts are abundant in the testis, suggesting a possible role during spermatogenesis.

Effects of colonization asymmetries on metapopulation persistence.

Ocean currents, prevailing winds, and the hierarchical structures of river networks are known to create asymmetries in re-colonization between habitat patches. The impacts of such asymmetries on metapopulation persistence are seldom considered, especially rarely in theoretical studies. Considering three classical models (the island, the stepping stone and the distance-dependent model), we explore how metapopulation persistence is affected by (i) asymmetry in dispersal strength, in which the colonization rate between two patches differs in direction, and (ii) asymmetry in connectivity, in which the overall colonization pattern displays asymmetry (circulating or dendritic networks). Viability can be drastically reduced when directional bias in dispersal strength is higher than 25%. Re-colonization patterns that allow for strong local connectivity provide the highest persistence compared to systems that allow circulation. Finally, asymmetry has relatively weak effects when metapopulations maintain strong general connectivity.

El desarrollo humano: perfiles y perspectivas futuras

El PNUD calcula todos los años el Índice de Desarrollo Humano (IDH). El objetivo de este trabajo es analizar los antecedentes y las perspectivas futuras del desarrollo humano a partir de los datos de este índice durante el periodo 1970-2000. Esto es, comprobar, a partir de los datos pasados, si las diferencias entre el IDH de los países del mundo están aumentando y valorar la tendencia del IDH esperada para los próximos años. En definitiva, se trata de buscar respuesta a la siguiente pregunta ¿cómo se modificarán los niveles de desarrollo humano en el futuro? Para ello, se utiliza la metodología estadística del análisis dinámico de distribución mediante las cadenas de Markov.

Disentangling the formation of contrasting tree line physiognomies combining model selection and Bayesian parameterization for simulation models.

Alpine tree-line ecotones are characterized by marked changes at small spatial scales that may result in a variety of physiognomies. A set of alternative individual-based models was tested with data from four contrasting Pinus uncinata ecotones in the central Spanish Pyrenees to reveal the minimal subset of processes required for tree-line formation. A Bayesian approach combined with Markov chain Monte Carlo methods was employed to obtain the posterior distribution of model parameters, allowing the use of model selection procedures. The main features of real tree lines emerged only in models considering nonlinear responses in individual rates of growth or mortality with respect to the altitudinal gradient. Variation in tree-line physiognomy reflected mainly changes in the relative importance of these nonlinear responses, while other processes, such as dispersal limitation and facilitation, played a secondary role. Different nonlinear responses also determined the presence or absence of krummholz, in agreement with recent findings highlighting a different response of diffuse and abrupt or krummholz tree lines to climate change. The method presented here can be widely applied in individual-based simulation models and will turn model selection and evaluation in this type of models into a more transparent, effective, and efficient exercise.

Statistical Optimum Design of Heat Exchangers

The optimal design of a heat exchanger system is based on given model parameters together with given standard ranges for machine design variables. The goals set for minimizing the Life Cycle Cost (LCC) function which represents the price of the saved energy, for maximizing the momentary heat recovery output with given constraints satisfied and taking into account the uncertainty in the models were successfully done. Nondominated Sorting Genetic Algorithm II (NSGA-II) for the design optimization of a system is presented and implemented inMatlab environment. Markov ChainMonte Carlo (MCMC) methods are also used to take into account the uncertainty in themodels. Results show that the price of saved energy can be optimized. A wet heat exchanger is found to be more efficient and beneficial than a dry heat exchanger even though its construction is expensive (160 EUR/m2) compared to the construction of a dry heat exchanger (50 EUR/m2). It has been found that the longer lifetime weights higher CAPEX and lower OPEX and vice versa, and the effect of the uncertainty in the models has been identified in a simplified case of minimizing the area of a dry heat exchanger.

Stochastic variability of output levels in a pulp mill

Quite often, in the construction of a pulp mill involves establishing the size of tanks which will accommodate the material from the various processes in which case estimating the right tank size a priori would be vital. Hence, simulation of the whole production process would be worthwhile. Therefore, there is need to develop mathematical models that would mimic the behavior of the output from the various production units of the pulp mill to work as simulators. Markov chain models, Autoregressive moving average (ARMA) model, Mean reversion models with ensemble interaction together with Markov regime switching models are proposed for that purpose.