The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

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).

Forecasting and turning point predictions in a Bayesian panel VAR model

We provide methods for forecasting variables and predicting turning points in panel Bayesian VARs. We specify a flexible model which accounts for both interdependencies in the cross section and time variations in the parameters. Posterior distributions for the parameters are obtained for a particular type of diffuse, for Minnesota-type and for hierarchical priors. Formulas for multistep, multiunit point and average forecasts are provided. An application to the problem of forecasting the growth rate of output and of predicting turning points in the G-7 illustrates the approach. A comparison with alternative forecasting methods is also provided.

Cost-effectiveness of low-molecular-weight heparin for treatment of pulmonary embolism.

BACKGROUND: Low-molecular-weight heparin (LMWH) appears to be safe and effective for treating pulmonary embolism (PE), but its cost-effectiveness has not been assessed. METHODS: We built a Markov state-transition model to evaluate the medical and economic outcomes of a 6-day course with fixed-dose LMWH or adjusted-dose unfractionated heparin (UFH) in a hypothetical cohort of 60-year-old patients with acute submassive PE. Probabilities for clinical outcomes were obtained from a meta-analysis of clinical trials. Cost estimates were derived from Medicare reimbursement data and other sources. The base-case analysis used an inpatient setting, whereas secondary analyses examined early discharge and outpatient treatment with LMWH. Using a societal perspective, strategies were compared based on lifetime costs, quality-adjusted life-years (QALYs), and the incremental cost-effectiveness ratio. RESULTS: Inpatient treatment costs were higher for LMWH treatment than for UFH (dollar 13,001 vs dollar 12,780), but LMWH yielded a greater number of QALYs than did UFH (7.677 QALYs vs 7.493 QALYs). The incremental costs of dollar 221 and the corresponding incremental effectiveness of 0.184 QALYs resulted in an incremental cost-effectiveness ratio of dollar 1,209/QALY. Our results were highly robust in sensitivity analyses. LMWH became cost-saving if the daily pharmacy costs for LMWH were < dollar 51, if > or = 8% of patients were eligible for early discharge, or if > or = 5% of patients could be treated entirely as outpatients. CONCLUSION: For inpatient treatment of PE, the use of LMWH is cost-effective compared to UFH. Early discharge or outpatient treatment in suitable patients with PE would lead to substantial cost savings.

Cost effectiveness and cost utility of risedronate for osteoporosis treatment and fracture prevention in women: a Swiss perspective.

OBJECTIVES: To assess the incremental cost-effectiveness ratio (ICER) and incremental cost-utility ratio (ICUR) of risedronate compared to no intervention in postmenopausal osteoporotic women in a Swiss perspective. METHODS: A previously validated Markov model was populated with epidemiological and cost data specific to Switzerland and published utility values, and run on a population of 1,000 women of 70 years with established osteoporosis and previous vertebral fracture, treated over 5 years with risedronate 35 mg weekly or no intervention (base case), and five cohorts (according to age at therapy start) with eight risk factor distributions and three lengths of residual effects. RESULTS: In the base case population, the ICER of averting a hip fracture and the ICUR per quality-adjusted life year gained were both dominant. In the presence of a previous vertebral fracture, the ICUR was below euro45,000 (pound30,000) in all the scenarios. For all osteoporotic women>or=70 years of age with at least one risk factor, the ICUR was below euro45,000 or the intervention may even be cost saving. Age at the start of therapy and the fracture risk profile had a significant impact on results. CONCLUSION: Assuming a 2-year residual effect, that ICUR of risedronate in women with postmenopausal osteoporosis is below accepted thresholds from the age of 65 and even cost saving above the age of 70 with at least one risk factor.

Cost-effectiveness analysis of a European primary-care physician training in smoking cessation counseling.

BACKGROUND: Physician training in smoking cessation counseling has been shown to be effective as a means to increase quit success. We assessed the cost-effectiveness ratio of a smoking cessation counseling training programme. Its effectiveness was previously demonstrated in a cluster randomized, control trial performed in two Swiss university outpatients clinics, in which residents were randomized to receive training in smoking interventions or a control educational intervention. DESIGN AND METHODS: We used a Markov simulation model for effectiveness analysis. This model incorporates the intervention efficacy, the natural quit rate, and the lifetime probability of relapse after 1-year abstinence. We used previously published results in addition to hospital service and outpatient clinic cost data. The time horizon was 1 year, and we opted for a third-party payer perspective. RESULTS: The incremental cost of the intervention amounted to US$2.58 per consultation by a smoker, translating into a cost per life-year saved of US$25.4 for men and 35.2 for women. One-way sensitivity analyses yielded a range of US$4.0-107.1 in men and US$9.7-148.6 in women. Variations in the quit rate of the control intervention, the length of training effectiveness, and the discount rate yielded moderately large effects on the outcome. Variations in the natural cessation rate, the lifetime probability of relapse, the cost of physician training, the counseling time, the cost per hour of physician time, and the cost of the booklets had little effect on the cost-effectiveness ratio. CONCLUSIONS: Training residents in smoking cessation counseling is a very cost-effective intervention and may be more efficient than currently accepted tobacco control interventions.

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.

Missed appointments in an outpatient clinic for adolescents, an approach to predict the risk of missing.

PURPOSE: To predict the risk of an adolescent patient to miss an appointment, based on the previous appointments and on the characteristics of the patient and the appointment. METHODS: Two thousand one hundred ninety-three (1873 females) patients aged 12 to 20 years having scheduled at least four appointments were included. We assessed the rate of missed nonexcused appointments of each patient. Second, a Markovian multilevel model was used to predict the risk of defaulting. RESULTS: Forty-five percent of the patients have not missed even once, and 14% of females and 17% of males have missed &gt;25% of their appointments. Females show two types of behaviors (an abstract concept that groups individuals based on a combination of their appointment-keeping and their recorded type of healthcare need) depending on the diagnosis. Somatic, gynecology, violence, and counseling diagnoses are mostly grouped together. In this group, having already missed and having an appointment with a paramedical provider increases the risk of missing. In the second group (eating disorders and psychiatric diagnoses) having already missed and a longer delay between appointments influence the risk of missing, although the risk is lower for this latter group. Males only show one type of behavior regarding missed appointments. Having missed a previous appointment, being older, having cancelled the next to last appointment and the type of diagnosis explain the risk of missing. CONCLUSIONS: Patients who have already defaulted have a higher risk of defaulting again. Means of control regarding missed appointments should consequently focus on defaulters, to decrease the associated workload. Reminders could be a solution for the follow-up appointments scheduled with a long delay.

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.

Monte Carlo Methods in Parameter Estimation of Nonlinear Models

Yksi keskeisimmistä tehtävistä matemaattisten mallien tilastollisessa analyysissä on mallien tuntemattomien parametrien estimointi. Tässä diplomityössä ollaan kiinnostuneita tuntemattomien parametrien jakaumista ja niiden muodostamiseen sopivista numeerisista menetelmistä, etenkin tapauksissa, joissa malli on epälineaarinen parametrien suhteen. Erilaisten numeeristen menetelmien osalta pääpaino on Markovin ketju Monte Carlo -menetelmissä (MCMC). Nämä laskentaintensiiviset menetelmät ovat viime aikoina kasvattaneet suosiotaan lähinnä kasvaneen laskentatehon vuoksi. Sekä Markovin ketjujen että Monte Carlo -simuloinnin teoriaa on esitelty työssä siinä määrin, että menetelmien toimivuus saadaan perusteltua. Viime aikoina kehitetyistä menetelmistä tarkastellaan etenkin adaptiivisia MCMC menetelmiä. Työn lähestymistapa on käytännönläheinen ja erilaisia MCMC -menetelmien toteutukseen liittyviä asioita korostetaan. Työn empiirisessä osuudessa tarkastellaan viiden esimerkkimallin tuntemattomien parametrien jakaumaa käyttäen hyväksi teoriaosassa esitettyjä menetelmiä. Mallit kuvaavat kemiallisia reaktioita ja kuvataan tavallisina differentiaaliyhtälöryhminä. Mallit on kerätty kemisteiltä Lappeenrannan teknillisestä yliopistosta ja Åbo Akademista, Turusta.

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.

Computational methods for tactical simulations

Tämä taktiikan tutkimus keskittyy tietokoneavusteisen simuloinnin laskennallisiin menetelmiin, joita voidaan käyttää taktisen tason sotapeleissä. Työn tärkeimmät tuotokset ovat laskennalliset mallit todennäköisyyspohjaisen analyysin mahdollistaviin taktisen tason taistelusimulaattoreihin, joita voidaan käyttää vertailevaan analyysiin joukkue-prikaatitason tarkastelutilanteissa. Laskentamallit keskittyvät vaikuttamiseen. Mallit liittyvät vahingoittavan osuman todennäköisyyteen, jonka perusteella vaikutus joukossa on mallinnettu tilakoneina ja Markovin ketjuina. Edelleen näiden tulokset siirretään tapahtumapuuanalyysiin operaation onnistumisen todennäköisyyden osalta. Pienimmän laskentayksikön mallinnustaso on joukkue- tai ryhmätasolla, jotta laskenta-aika prikaatitason sotapelitarkasteluissa pysyisi riittävän lyhyenä samalla, kun tulokset ovat riittävän tarkkoja suomalaiseen maastoon. Joukkueiden mies- ja asejärjestelmävahvuudet ovat jakaumamuodossa, eivätkä yksittäisiä lukuja. Simuloinnin integroinnissa voidaan käyttää asejärjestelmäkohtaisia predictor corrector –parametreja, mikä mahdollistaa aika-askelta lyhytaikaisempien taistelukentän ilmiöiden mallintamisen. Asemallien pohjana ovat aiemmat tutkimukset ja kenttäkokeet, joista osa kuuluu tähän väitöstutkimukseen. Laskentamallien ohjelmoitavuus ja käytettävyys osana simulointityökalua on osoitettu tekijän johtaman tutkijaryhmän ohjelmoiman ”Sandis”- taistelusimulointiohjelmiston avulla, jota on kehitetty ja käytetty Puolustusvoimien Teknillisessä Tutkimuslaitoksessa. Sandikseen on ohjelmoitu karttakäyttöliittymä ja taistelun kulkua simuloivia laskennallisia malleja. Käyttäjä tai käyttäjäryhmä tekee taktiset päätökset ja syöttää nämä karttakäyttöliittymän avulla simulointiin, jonka tuloksena saadaan kunkin joukkuetason peliyksikön tappioiden jakauma, keskimääräisten tappioiden osalta kunkin asejärjestelmän aiheuttamat tappiot kuhunkin maaliin, ammuskulutus ja radioyhteydet ja niiden tila sekä haavoittuneiden evakuointi-tilanne joukkuetasolta evakuointisairaalaan asti. Tutkimuksen keskeisiä tuloksia (kontribuutio) ovat 1) uusi prikaatitason sotapelitilanteiden laskentamalli, jonka pienin yksikkö on joukkue tai ryhmä; 2) joukon murtumispisteen määritys tappioiden ja haavoittuneiden evakuointiin sitoutuvien taistelijoiden avulla; 3) todennäköisyyspohjaisen riskianalyysin käyttömahdollisuus vertailevassa tutkimuksessa sekä 4) kokeellisesti testatut tulen vaikutusmallit ja 5) toimivat integrointiratkaisut. Työ rajataan maavoimien taistelun joukkuetason todennäköisyysjakaumat luovaan laskentamalliin, kenttälääkinnän malliin ja epäsuoran tulen malliin integrointimenetelmineen sekä niiden antamien tulosten sovellettavuuteen. Ilmasta ja mereltä maahan -asevaikutusta voidaan tarkastella, mutta ei ilma- ja meritaistelua. Menetelmiä soveltavan Sandis -ohjelmiston malleja, käyttötapaa ja ohjelmistotekniikkaa kehitetään edelleen. Merkittäviä jatkotutkimuskohteita mallinnukseen osalta ovat muun muassa kaupunkitaistelu, vaunujen kaksintaistelu ja maaston vaikutus tykistön tuleen sekä materiaalikulutuksen arviointi.