Bayesian parametric bootstrap for models with intractable likelihoods

In this paper it is demonstrated how the Bayesian parametric bootstrap can be adapted to models with intractable likelihoods. The approach is most appealing when the semi-automatic approximate Bayesian computation (ABC) summary statistics are selected. After a pilot run of ABC, the likelihood-free parametric bootstrap approach requires very few model simulations to produce an approximate posterior, which can be a useful approximation in its own right. An alternative is to use this approximation as a proposal distribution in ABC algorithms to make them more efficient. In this paper, the parametric bootstrap approximation is used to form the initial importance distribution for the sequential Monte Carlo and the ABC importance and rejection sampling algorithms. The new approach is illustrated through a simulation study of the univariate g-and- k quantile distribution, and is used to infer parameter values of a stochastic model describing expanding melanoma cell colonies.

Evaluation of TRMM PR Sampling Error Over a Subtropical Basin Using Bootstrap Technique

Quantitative use of satellite-derived rainfall products for various scientific applications often requires them to be accompanied with an error estimate. Rainfall estimates inferred from low earth orbiting satellites like the Tropical Rainfall Measuring Mission (TRMM) will be subjected to sampling errors of nonnegligible proportions owing to the narrow swath of satellite sensors coupled with a lack of continuous coverage due to infrequent satellite visits. The authors investigate sampling uncertainty of seasonal rainfall estimates from the active sensor of TRMM, namely, Precipitation Radar (PR), based on 11 years of PR 2A25 data product over the Indian subcontinent. In this paper, a statistical bootstrap technique is investigated to estimate the relative sampling errors using the PR data themselves. Results verify power law scaling characteristics of relative sampling errors with respect to space-time scale of measurement. Sampling uncertainty estimates for mean seasonal rainfall were found to exhibit seasonal variations. To give a practical example of the implications of the bootstrap technique, PR relative sampling errors over a subtropical river basin of Mahanadi, India, are examined. Results reveal that the bootstrap technique incurs relative sampling errors < 33% (for the 2 degrees grid), < 36% (for the 1 degrees grid), < 45% (for the 0.5 degrees grid), and < 57% (for the 0.25 degrees grid). With respect to rainfall type, overall sampling uncertainty was found to be dominated by sampling uncertainty due to stratiform rainfall over the basin. The study compares resulting error estimates to those obtained from latin hypercube sampling. Based on this study, the authors conclude that the bootstrap approach can be successfully used for ascertaining relative sampling errors offered by TRMM-like satellites over gauged or ungauged basins lacking in situ validation data. This technique has wider implications for decision making before incorporating microwave orbital data products in basin-scale hydrologic modeling.

An important pitfall of pseudo-static finite element analysis

Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified methods of dynamic analysis and time domain analysis. The pseudo-static FE approach can accurately model the in situ, stresses prior to seismic loading (when it follows a static analysis simulating the construction sequence) is relatively simple and not as computationally expensive as the time domain approach. However this method should be used with caution as the results can be sensitive to the choice of the mesh dimensions. In this paper two simple examples of pseudo-static finite element analysis are examined parametrically, a homogeneous slope and a cantilever retaining wall, exploring the sensitivity of the pseudo-static analysis results on the adopted mesh size. The mesh dependence was found to be more pronounced for problems with high critical seismic coefficients values (e.g. gentle slopes or small walls), as in these cases a generalised layer failure mechanism is developed simultaneously with the slope or wall mechanism. In general the mesh width was found not to affect notably the predicted value of critical seismic coefficient but to have a major impact on the predicted movements. © 2012 Elsevier Ltd.

Heller's myotomy and pneumatic dilatation in the treatment of achalasia: a population-based case–control study assessing long-term quality of life

Long-term health-related quality-of-life (HRQL) outcomes have not been widely reported in the
treatment of achalasia. The aims of this study were to examine long-term disease-specific and general HRQL in
achalasia patients using a population-based case–control method, and to assess HRQL between treatment interventions.
Manometrically diagnosed achalasia cases (n = 120) were identified and matched with controls (n = 115)
using a population-based approach. Participants completed general (SF-12) and disease-specific (Achalasia Severity
Questionnaire [ASQ]) HRQL questionnaires, as appropriate, in a structured interview. Mean composite scores
for SF-12 (Mental Component Summary score [MCS-12] and Physical Component Summary score [PCS-12]) and
ASQ were compared between cases and controls, or between intervention groups, using an independent t-test.
Adjusted mean differences in HRQL scores were evaluated using a linear regression model. Achalasia cases were
treated with a Heller’s myotomy (n = 43), pneumatic dilatation (n = 44), or both modalities (n = 33). The median
time from last treatment to HRQL assessment was 5.7 years (interquartile range 2.4–11.5). Comparing achalasia
patients with controls, PCS-12 was significantly worse (40.9 vs. 44.2, P = 0.01), but MCS-12 was similar. However,
both PCS-12 (39.9 vs. 44.2, P = 0.03) and MCS-12 (46.7 vs. 53.5, P = 0.004) were significantly impaired in those
requiring dual treatment compared with controls. Overall however, there was no difference in adjusted HRQL
between patients treated with Heller’s myotomy, pneumatic dilatation or both treatment modalities. In summary,
despite treatment achalasia patients have significantly worse long-term physical HRQL compared with population
controls. No HRQL differences were observed between the treatment modalities to suggest a benefit of one
treatment over another.

Optimisation de l’administration des médicaments chez les enfants transplantés grâce à la pharmacocinétique de population

Ce travail de thèse porte sur l’application de la pharmacocinétique de population dans le but d’optimiser l’utilisation de certains médicaments chez les enfants immunosupprimés et subissant une greffe. Parmi les différents médicaments utilisés chez les enfants immunosupprimés, l’utilisation du busulfan, du tacrolimus et du voriconazole reste problématique, notamment à cause d’une très grande variabilité interindividuelle de leur pharmacocinétique rendant nécessaire l’individualisation des doses par le suivi thérapeutique pharmacologique. De plus, ces médicaments n’ont pas fait l’objet d’études chez les enfants et les doses sont adaptées à partir des adultes. Cette dernière pratique ne prend pas en compte les particularités pharmacologiques qui caractérisent l’enfant tout au long de son développement et rend illusoire l’extrapolation aux enfants des données acquises chez les adultes. Les travaux effectués dans le cadre de cette thèse ont étudié successivement la pharmacocinétique du busulfan, du voriconazole et du tacrolimus par une approche de population en une étape (modèles non-linéaires à effets mixtes). Ces modèles ont permis d’identifier les principales sources de variabilités interindividuelles sur les paramètres pharmacocinétiques. Les covariables identifiées sont la surface corporelle et le poids. Ces résultats confirment l’importance de tenir en compte l’effet de la croissance en pédiatrie. Ces paramètres ont été inclus de façon allométrique dans les modèles. Cette approche permet de séparer l’effet de la mesure anthropométrique d’autres covariables et permet la comparaison des paramètres pharmacocinétiques en pédiatrie avec ceux des adultes. La prise en compte de ces covariables explicatives devrait permettre d’améliorer la prise en charge a priori des patients. Ces modèles développés ont été évalués pour confirmer leur stabilité, leur performance de simulation et leur capacité à répondre aux objectifs initiaux de la modélisation. Dans le cas du busulfan, le modèle validé a été utilisé pour proposer par simulation une posologie qui améliorerait l’atteinte de l’exposition cible, diminuerait l’échec thérapeutique et les risques de toxicité. Le modèle développé pour le voriconazole, a permis de confirmer la grande variabilité interindividuelle dans sa pharmacocinétique chez les enfants immunosupprimés. Le nombre limité de patients n’a pas permis d’identifier des covariables expliquant cette variabilité. Sur la base du modèle de pharmacocinétique de population du tacrolimus, un estimateur Bayesien a été mis au point, qui est le premier dans cette population de transplantés hépatiques pédiatriques. Cet estimateur permet de prédire les paramètres pharmacocinétiques et l’exposition individuelle au tacrolimus sur la base d’un nombre limité de prélèvements. En conclusion, les travaux de cette thèse ont permis d’appliquer la pharmacocinétique de population en pédiatrie pour explorer les caractéristiques propres à cette population, de décrire la variabilité pharmacocinétique des médicaments utilisés chez les enfants immunosupprimés, en vue de l’individualisation du traitement. Les outils pharmacocinétiques développés s’inscrivent dans une démarche visant à diminuer le taux d'échec thérapeutique et l’incidence des effets indésirables ou toxiques chez les enfants immunosupprimés suite à une transplantation.

Population genetic models can be used to study the evolution of the interacting behaviours of parents and their progeny

An example of the evolution of the interacting behaviours of parents and progeny is studied using iterative equations linking the frequencies of the gametes produced by the progeny to the frequencies of the gametes in the parental generation. This population genetics approach shows that a model in which both behaviours are determined by a single locus can lead to a stable equilibrium in which the two behaviours continue to segregate. A model in which the behaviours are determined by genes at two separate loci leads eventually to fixation of the alleles at both loci but this can take many generations of selection. Models of the type described in this paper will be needed to understand the evolution of complex behaviour when genomic or experimental information is available about the genetic determinants of behaviour and the selective values of different genomes. (c) 2007 Elsevier Inc. All rights reserved.

A novel MAP-MRF approach for multispectral image contextual classification using combination of suboptimal iterative algorithms

In this paper we present a novel approach for multispectral image contextual classification by combining iterative combinatorial optimization algorithms. The pixel-wise decision rule is defined using a Bayesian approach to combine two MRF models: a Gaussian Markov Random Field (GMRF) for the observations (likelihood) and a Potts model for the a priori knowledge, to regularize the solution in the presence of noisy data. Hence, the classification problem is stated according to a Maximum a Posteriori (MAP) framework. In order to approximate the MAP solution we apply several combinatorial optimization methods using multiple simultaneous initializations, making the solution less sensitive to the initial conditions and reducing both computational cost and time in comparison to Simulated Annealing, often unfeasible in many real image processing applications. Markov Random Field model parameters are estimated by Maximum Pseudo-Likelihood (MPL) approach, avoiding manual adjustments in the choice of the regularization parameters. Asymptotic evaluations assess the accuracy of the proposed parameter estimation procedure. To test and evaluate the proposed classification method, we adopt metrics for quantitative performance assessment (Cohen`s Kappa coefficient), allowing a robust and accurate statistical analysis. The obtained results clearly show that combining sub-optimal contextual algorithms significantly improves the classification performance, indicating the effectiveness of the proposed methodology. (C) 2010 Elsevier B.V. All rights reserved.

Suicide in breast cancer patients: an individual-centered approach provides insight beyond epidemiology

Epidemiologic studies have identified increased suicide rates among breast cancer (BC) patients. The population-based approach, however, has considerable methodic shortcomings. None of the studies have been carried out in a prospective manner and none reported suicide rates from a country in which physician-assisted suicide (PAS) is legal.

Effect of the Legume Plant Host on Rhizobial Soil Population Dynamics

Rhizobium leguminosarum bv viciae (Rlv) is a soil bacterium able to establish specific root-nodule symbioses with legumes of four different genera: Pisum, Vicia, Lens and Lathyrus. Rlv isolates from nodules of any of these legumes can nodulate any of them; however, it has been shown that plants select specific rhizobial genotypes from those present in the soil (1,2). We have previously shown this at the genomic level by following a population genomics approach. Pool genomic sequences from 100 isolates from each of four plant species: P. sativum, L. culinaris, V. faba and V. sativa, show different, specific profiles at the single nucleotide polymorphism (SNP) level for relevant genes. In this work, the extent of Rlv selection from a well-characterized soil population by different legume plant hosts: P. sativum, L. culinaris, V. faba and V. sativa, after a medium-term mesocosm study is described. Direct soil isolates from each of these mesocosm studies have been tested for specific rhizobial genes (glnII and fnrN) and symbiotic genes (nodC and nifH). Different populations were characterized further by Sanger sequencing of both the rpoB phylogenetic marker gene and the symbiotic genes nodC and nifH. The distribution and size of the rhizobial population for each legume host showed changes during the medium-term mesocosm study. Particularly, a non-symbiotic group of rhizobia was enriched by all four hosts, in contrast to the symbiotic rhizobia profile, which was specific for each legume plant host.

Physical and environmental determinants

This chapter describes physical and environmental determinants of the health of Australians, providing a background to the development of successful public health activity. Health determinants are the biomedical, genetic, behavioural, socio-economic and environmental factors that impact on health and wellbeing. These determinants can be influenced by interventions and by resources and systems (AIHW 2006). Many factors combine to affect the health of individuals and communities. People’s circumstances and the environment determine whether the population is healthy or not. Factors such as where people live, the state of their environment, genetics, their education level and income, and their relationships with friends and family, all are likely to impact on their health. The determinants of population health reflect the context of people’s lives; however, people are very unlikely to be able to control many of these determinants (WHO 2007). This chapter and Chapter 6 illustrate how various determinants can relate to, and influence other determinants, as well as health and wellbeing. We believe it is particularly important to provide an understanding of determinants and their relationship to health and illness in order to provide a structure in which a broader conceptualisation of health can be placed. Determinants of health do not exist in isolation from one another. More frequently they work together in a complex system. What is clear to anyone who works in public health is that many factors impact on the health and wellbeing of people. For example, in the next chapter we discuss factors such as living and working conditions, social support, ethnicity and class, income, housing, work stress and the impact of education on the length and quality of people’s lives. In 1974, the influential ‘Lalonde Report’ (Lalonde 1974) described key factors that impact on health status. These factors included lifestyle, environment, human biology and health services. Taking a population health approach builds on the Lalonde Report, and recognises that a range of factors, such as living and working conditions and the distribution of wealth in society, interact to determine the health status of a population. Tackling health determinants has great potential to reduce the burden of disease and promote the health of the general population. In summary, we understand very clearly now that health is determined by the complex interactions between individual characteristics, social and economic factors and physical environments; the entire range of factors that impact on health must be addressed if we are to make significant gains in population health, and focussing interventions on the health of the population or significant sub-populations can achieve important health gains. In 2007, the Australian Government included in the list of National Health Priority Areas the following health issues: cancer control, injury prevention and control, cardiovascular health, diabetes mellitus, mental health, asthma, and arthritis and musculoskeletal conditions. The National Health Priority Areas set the agenda for the Commonwealth, States and Territories, Local Governments and not-for-profit organisations to place attention on those areas considered to be the major foci for action. Many of these health issues are discussed in this chapter and the following chapter.

Prevalence and potential fitness cost of weak phosphine resistance in Tribolium castaneum (Herbst) in eastern Australia

The prevalence of resistance to phosphine in the rust-red flour beetle, Tribolium castaneum, from eastern Australia was investigated, as well as the potential fitness cost of this type of resistance. Discriminating dose tests on 115 population samples collected from farms from 2006 to 2010 showed that populations containing insects with the weakly resistant phenotype are common in eastern Australia (65.2 of samples), although the frequency of resistant phenotypes within samples was typically low (median of 2.3). The population cage approach was used to investigate the possibility that carrying the alleles for weak resistance incurs a fitness cost. Hybridized populations were initiated using a resistant strain and either of two different susceptible strains. There was no evidence of a fitness cost based on the frequency of susceptible phenotypes in hybridized populations that were reared for seven generations without exposure to phosphine. This suggests that resistant alleles will tend to persist in field populations that have undergone selection even if selection pressure is removed. The prevalence of resistance is a warning that this species has been subject to considerable selection pressure and that effective resistance management practices are needed to address this problem. The resistance prevalence data also provide a basis against which to measure management success.

Modeling of agglomeration inside a droplet with nanosuspensions in an acoustic field

This short communication reports results of particle agglomeration details of an acoustically levitated nanosilica droplet. The droplet undergoes thermo-physical and morphological changes under external heating load (convective or radiative) forming different solid structures due to particle agglomeration. We report an agglomeration model based on population balance approach coupled with species and energy conservation equations in the liquid phase and compare it with the experimentally observed structure formations using high speed photography. The analysis is able to predict similar spherical bowl shaped morphologies as observed experimentally using scanning electron microscopy and laser induced fluorescence. (C) 2012 Elsevier Ltd. All rights reserved.

Algoritmos multi-objetivos para la optimización de una secuencia de matrices generadoras en computación cuántica topológica

Quantum Computing is a relatively modern field which simulates quantum computation conditions. Moreover, it can be used to estimate which quasiparticles would endure better in a quantum environment. Topological Quantum Computing (TQC) is an approximation for reducing the quantum decoherence problem1, which is responsible for error appearance in the representation of information. This project tackles specific instances of TQC problems using MOEAs (Multi-objective Optimization Evolutionary Algorithms). A MOEA is a type of algorithm which will optimize two or more objectives of a problem simultaneously, using a population based approach. We have implemented MOEAs that use probabilistic procedures found in EDAs (Estimation of Distribution Algorithms), since in general, EDAs have found better solutions than ordinary EAs (Evolutionary Algorithms), even though they are more costly. Both, EDAs and MOEAs are population-based algorithms. The objective of this project was to use a multi-objective approach in order to find good solutions for several instances of a TQC problem. In particular, the objectives considered in the project were the error approximation and the length of a solution. The tool we used to solve the instances of the problem was the multi-objective framework PISA. Because PISA has not too much documentation available, we had to go through a process of reverse-engineering of the framework to understand its modules and the way they communicate with each other. Once its functioning was understood, we began working on a module dedicated to the braid problem. Finally, we submitted this module to an exhaustive experimentation phase and collected results.

Bone mineral density in relation to medical and lifestyle risk factors for osteoporosis in premenopausal, menopausal and postmenopausal women in general practice

Background: Interest in the prevention of osteoporosis is increasing and thus there is a need for an acceptable osteoporosis prevention programme in general practice. AIM. A study was undertaken to identify a cohort of middle-aged women attending a general practice who would be eligible for a longitudinal study looking at bone mineral density, osteoporosis and the effectiveness of hormone replacement therapy. This study aimed to describe the relationship between medical and lifestyle risk factors for osteoporosis and the initial bone density measurements in this group of women. METHOD. A health visitor administered a questionnaire to women aged between 48 and 52 years registered with a Belfast general practice. The main outcome measures were menopausal status, presence of medical and lifestyle risk factors and bone mineral density measurements. RESULTS. A total of 358 women our of 472 (76%) took part in the study which was conducted in 1991 and 1992. A highly significant difference was found between the mean bone mineral density of premenopausal, menopausal and postmenopausal women within the narrow study age range, postmenopausal women having the lowest bone mineral density. A significant relationship was found between body mass index and bone mineral density, a greater bone mineral density being found among women with a higher body mass index. Risk factors such as smoking and sedentary lifestyle were common (reported by approximately one third of respondents) but a poor relationship was found between these two and all the other risk factors and bone mineral density in this age group. CONCLUSION. Risk of osteoporosis cannot be identified by the presence of risk factors in women aged between 48 and 52 years. In terms of a current prevention strategy for general practice it would be better to take a population-based approach except for those women known to be at high risk of osteoporosis: women with early menopause or those who have had an oophorectomy.

Parenting Programmes: The Best Available Evidence

Parenting programmes have been provided to a wide range of child and parent groups across a number of countries, but are they effective? This aim of this paper is to examine the findings from a number of systematic reviews that summarise the best available research evidence on the impact of these programmes on a range of parental and child outcomes. In addition to examining the findings from systematic reviews, the paper also takes a selective look at the uptake of parenting programmes in the United Kingdom, the evidence for effectiveness and the efficacy of adopting a population-based approach to parent education.

The findings from systematic reviews indicate that parenting programmes can have a positive impact on a range of outcomes, including improved child behaviour, increased maternal self-esteem and relationship adjustment, improved mother–child interaction and knowledge and decreased maternal depression and stress. While there is a need for greater evaluation of the long-term impact of these programmes, preliminary evidence indicates that these positive results are maintained over time, with group-based, behaviourally orientated programmes tending to be more effective.

While several recent trials indicate that that these programmes can be effective within the United Kingdom, high drop-out rates may mean that they only reach a minority of parents. However, multi-level parent education strategies such as the Australian Triple P Positive Parenting Strategy that incorporate an array of mediums aimed at different levels of need may provide an opportunity to reach a wider range of parents. This approach is currently being evaluated in order to ascertain whether it is effective in improving child outcomes in the general population.

While there is no coherent strategy for parent training across the United Kingdom, within the Northern Ireland context there is a move towards the development of a family support strategy. While uptake of parent education and training is currently unknown the best available evidence highlights the positive impact that parent training can have, suggesting the importance of including parent education as one aspect of this strategy