The importance of modelling heterogeneity in complex disease: application to NIMH Schizophrenia Genetics Initiative data

As for other complex diseases, linkage analyses of schizophrenia (SZ) have produced evidence for numerous chromosomal regions, with inconsistent results reported across studies. The presence of locus heterogeneity appears likely and may reduce the power of linkage analyses if homogeneity is assumed. In addition, when multiple heterogeneous datasets are pooled, inter-sample variation in the proportion of linked families (alpha) may diminish the power of the pooled sample to detect susceptibility loci, in spite of the larger sample size obtained. We compare the significance of linkage findings obtained using allele-sharing LOD scores (LOD(exp))-which assume homogeneity-and heterogeneity LOD scores (HLOD) in European American and African American NIMH SZ families. We also pool these two samples and evaluate the relative power of the LOD(exp) and two different heterogeneity statistics. One of these (HLOD-P) estimates the heterogeneity parameter alpha only in aggregate data, while the second (HLOD-S) determines alpha separately for each sample. In separate and combined data, we show consistently improved performance of HLOD scores over LOD(exp). Notably, genome-wide significant evidence for linkage is obtained at chromosome 10p in the European American sample using a recessive HLOD score. When the two samples are combined, linkage at the 10p locus also achieves genome-wide significance under HLOD-S, but not HLOD-P. Using HLOD-S, improved evidence for linkage was also obtained for a previously reported region on chromosome 15q. In linkage analyses of complex disease, power may be maximised by routinely modelling locus heterogeneity within individual datasets, even when multiple datasets are combined to form larger samples.

Optimal dynamic inversion-based semi-active control of benchmark bridge using MR dampers

A new two-stage state feedback control design approach has been developed to monitor the voltage supplied to magnetorheological (MR) dampers for semi-active vibration control of the benchmark highway bridge. The first stage contains a primary controller, which provides the force required to obtain a desired closed-loop response of the system. In the second stage, an optimal dynamic inversion (ODI) approach has been developed to obtain the amount of voltage to be supplied to each of the MR dampers such that it provides the required force prescribed by the primary controller. ODI is formulated by optimization with dynamic inversion, such that an optimal voltage is supplied to each damper in a set. The proposed control design has been simulated for both phase-I and phase-II study of the recently developed benchmark highway bridge problem. The efficiency of the proposed controller is analyzed in terms of the performance indices defined in the benchmark problem definition. Simulation results demonstrate that the proposed approach generally reduces peak response quantities over those obtained from the sample semi-active controller, although some response quantities have been seen to be increasing. Overall, the proposed control approach is quite competitive as compared with the sample semi-active control approach.

Surface instability in a semi-infinite plasma in the presence of high frequency fields

The surface instability of a semi-infinite plasma immersed in a high frequency field is investigated. When the natural Langmuir frequency of the surface is nearly equal to the frequency of the high frequency field, the dispersion relation predicts build-up of oscillations with a growth rate comparable with the real part of the frequency. Threshold values above which the instability is possible are derived.

An explicit finite element modelling method for masonry walls under out-of-plane loading

An explicit finite element modelling method is formulated using a layered shell element to examine the behaviour of masonry walls subject to out-of-plane loading. Masonry is modelled as a homogenised material with distinct directional properties that are calibrated from datasets of a “C” shaped wall tested under pressure loading applied to its web. The predictions of the layered shell model have been validated using several out-of-plane experimental datasets reported in the literature. Profound influence of support conditions, aspect ratio, pre-compression and opening to the strength and ductility of masonry walls is exhibited from the sensitivity analyses performed using the model.

Modelling temperature, photoperiod and vernalization responses of Brunonia australis (Goodeniaceae) and Calandrinia sp (Portulacaceae) to predict flowering time

Crop models for herbaceous ornamental species typically include functions for temperature and photoperiod responses, but very few incorporate vernalization, which is a requirement of many traditional crops. This study investigated the development of floriculture crop models, which describe temperature responses, plus photoperiod or vernalization requirements, using Australian native ephemerals Brunonia australis and Calandrinia sp. A novel approach involved the use of a field crop modelling tool, DEVEL2. This optimization program estimates the parameters of selected functions within the development rate models using an iterative process that minimizes sum of squares residual between estimated and observed days for the phenological event. Parameter profiling and jack-knifing are included in DEVEL2 to remove bias from parameter estimates and introduce rigour into the parameter selection process. Development rate of B. australis from planting to first visible floral bud (VFB) was predicted using a multiplicative approach with a curvilinear function to describe temperature responses and a broken linear function to explain photoperiod responses. A similar model was used to describe the development rate of Calandrinia sp., except the photoperiod function was replaced with an exponential vernalization function, which explained a facultative cold requirement and included a coefficient for determining the vernalization ceiling temperature. Temperature was the main environmental factor influencing development rate for VFB to anthesis of both species and was predicted using a linear model. The phenology models for B. australis and Calandrinia sp. described development rate from planting to VFB and from VFB to anthesis in response to temperature and photoperiod or vernalization and may assist modelling efforts of other herbaceous ornamental plants. In addition to crop management, the vernalization function could be used to identify plant communities most at risk from predicted increases in temperature due to global warming.

A participatory whole farm modelling approach to understand impacts and increase preparedness to climate change in Australia

This study presents the use of a whole farm model in a participatory modelling research approach to examine the sensitivity of four contrasting case study farms to a likely climate change scenario. The newly generated information was used to support discussions with the participating farmers in the search for options to design more profitable and sustainable farming systems in Queensland Australia. The four case studies contrasted in key systems characteristics: opportunism in decision making, i.e. flexible versus rigid crop rotations; function, i.e. production of livestock or crops; and level of intensification, i.e. dryland versus irrigated agriculture. Tested tactical and strategic changes under a baseline and climate change scenario (CCS) involved changes in the allocation of land between cropping and grazing enterprises, alternative allocations of limited irrigation water across cropping enterprises, and different management rules for planting wheat and sorghum in rainfed cropping. The results show that expected impacts from a likely climate change scenario were evident in the following increasing order: the irrigated cropping farm case study, the cropping and grazing farm, the more opportunistic rainfed cropping farm and the least opportunistic rainfed cropping farm. We concluded that in most cases the participating farmers were operating close to the efficiency frontier (i.e. in the relationship between profits and risks). This indicated that options to adapt to climate change might need to evolve from investments in the development of more innovative cropping and grazing systems and/or transformational changes on existing farming systems. We expect that even though assimilating expected changes in climate seems to be rather intangible and premature for these farmers, as innovations are developed, adaptation is likely to follow quickly. The multiple interactions among farm management components in complex and dynamic farm businesses operating in a variable and changing climate, make the use of whole farm participatory modelling approaches valuable tools to quantify benefits and trade-offs from alternative farming systems designs in the search for improved profitability and resilience.

Individual-based modelling of the efficacy of fumigation tactics to control lesser grain borer (Rhyzopertha dominica) in stored grain

Increasing resistance to phosphine (PH 3) in insect pests, including lesser grain borer (Rhyzopertha dominica) has become a critical issue, and development of effective and sustainable strategies to manage resistance is crucial. In practice, the same grain store may be fumigated multiple times, but usually for the same exposure period and concentration. Simulating a single fumigation allows us to look more closely at the effects of this standard treatment.We used an individual-based, two-locus model to investigate three key questions about the use of phosphine fumigant in relation to the development of PH 3 resistance. First, which is more effective for insect control; long exposure time with a low concentration or short exposure period with a high concentration? Our results showed that extending exposure duration is a much more efficient control tactic than increasing the phosphine concentration. Second, how long should the fumigation period be extended to deal with higher frequencies of resistant insects in the grain? Our results indicated that if the original frequency of resistant insects is increased n times, then the fumigation needs to be extended, at most, n days to achieve the same level of insect control. The third question is how does the presence of varying numbers of insects inside grain storages impact the effectiveness of phosphine fumigation? We found that, for a given fumigation, as the initial population number was increased, the final survival of resistant insects increased proportionally. To control initial populations of insects that were n times larger, it was necessary to increase the fumigation time by about n days. Our results indicate that, in a 2-gene mediated resistance where dilution of resistance gene frequencies through immigration of susceptibles has greater effect, extending fumigation times to reduce survival of homozygous resistant insects will have a significant impact on delaying the development of resistance. © 2012 Elsevier Ltd.

Microclimate and gas emissions in dairy buildings : Instrumentation, theory and measurements

The aim of this thesis was to develop measurement techniques and systems for measuring air quality and to provide information about air quality conditions and the amount of gaseous emissions from semi-insulated and uninsulated dairy buildings in Finland and Estonia. Specialization and intensification in livestock farming, such as in dairy production, is usually accompanied by an increase in concentrated environmental emissions. In addition to high moisture, the presence of dust and corrosive gases, and widely varying gas concentrations in dairy buildings, Finland and Estonia experience winter temperatures reaching below -40 ºC and summer temperatures above +30 ºC. The adaptation of new technologies for long-term air quality monitoring and measurement remains relatively uncommon in dairy buildings because the construction and maintenance of accurate monitoring systems for long-term use are too expensive for the average dairy farmer to afford. Though the documentation of accurate air quality measurement systems intended mainly for research purposes have been made in the past, standardised methods and the documentation of affordable systems and simple methods for performing air quality and emissions measurements in dairy buildings are unavailable. In this study, we built three measurement systems: 1) a Stationary system with integrated affordable sensors for on-site measurements, 2) a Wireless system with affordable sensors for off-site measurements, and 3) a Mobile system consisting of expensive and accurate sensors for measuring air quality. In addition to assessing existing methods, we developed simplified methods for measuring ventilation and emission rates in dairy buildings. The three measurement systems were successfully used to measure air quality in uninsulated, semi-insulated, and fully-insulated dairy buildings between the years 2005 and 2007. When carefully calibrated, the affordable sensors in the systems gave reasonably accurate readings. The spatial air quality survey showed high variation in microclimate conditions in the dairy buildings measured. The average indoor air concentration for carbon dioxide was 950 ppm, for ammonia 5 ppm, for methane 48 ppm, for relative humidity 70%, and for inside air velocity 0.2 m/s. The average winter and summer indoor temperatures during the measurement period were -7º C and +24 ºC for the uninsulated, +3 ºC and +20 ºC for the semi-insulated and +10 ºC and +25 ºC for the fully-insulated dairy buildings. The measurement results showed that the uninsulated dairy buildings had lower indoor gas concentrations and emissions compared to fully insulated buildings. Although occasionally exceeded, the ventilation rates and average indoor air quality in the dairy buildings were largely within recommended limits. We assessed the traditional heat balance, moisture balance, carbon dioxide balance and direct airflow methods for estimating ventilation rates. The direct velocity measurement for the estimation of ventilation rate proved to be impractical for naturally ventilated buildings. Two methods were developed for estimating ventilation rates. The first method is applicable in buildings in which the ventilation can be stopped or completely closed. The second method is useful in naturally ventilated buildings with large openings and high ventilation rates where spatial gas concentrations are heterogeneously distributed. The two traditional methods (carbon dioxide and methane balances), and two newly developed methods (theoretical modelling using Fick s law and boundary layer theory, and the recirculation flux-chamber technique) were used to estimate ammonia emissions from the dairy buildings. Using the traditional carbon dioxide balance method, ammonia emissions per cow from the dairy buildings ranged from 7 g day-1 to 35 g day-1, and methane emissions per cow ranged from 96 g day-1 to 348 g day-1. The developed methods proved to be as equally accurate as the traditional methods. Variation between the mean emissions estimated with the traditional and the developed methods was less than 20%. The developed modelling procedure provided sound framework for examining the impact of production systems on ammonia emissions in dairy buildings.

A Paddock to reef monitoring and modelling framework for the Great Barrier Reef: Paddock and catchment component

Targets for improvements in water quality entering the Great Barrier Reef (GBR) have been set through the Reef Water Quality Protection Plan (Reef Plan). To measure and report on progress towards the targets set a program has been established that combines monitoring and modelling at paddock through to catchment and reef scales; the Paddock to Reef Integrated Monitoring, Modelling and Reporting Program (Paddock to Reef Program). This program aims to provide evidence of links between land management activities, water quality and reef health. Five lines of evidence are used: the effectiveness of management practices to improve water quality; the prevalence of management practice adoption and change in catchment indicators; long-term monitoring of catchment water quality; paddock & catchment modelling to provide a relative assessment of progress towards meeting targets; and finally marine monitoring of GBR water quality and reef ecosystem health. This paper outlines the first four lines of evidence. (C) 2011 Elsevier Ltd. All rights reserved.

The cost-effectiveness of exercise referral schemes

Background Exercise referral schemes (ERS) aim to identify inactive adults in the primary care setting. The primary care professional refers the patient to a third party service, with this service taking responsibility for prescribing and monitoring an exercise programme tailored to the needs of the patient. This paper examines the cost-effectiveness of ERS in promoting physical activity compared with usual care in primary care setting. Methods A decision analytic model was developed to estimate the cost-effectiveness of ERS from a UK NHS perspective. The costs and outcomes of ERS were modelled over the patient's lifetime. Data were derived from a systematic review of the literature on the clinical and cost-effectiveness of ERS, and on parameter inputs in the modelling framework. Outcomes were expressed as incremental cost per quality-adjusted life-year (QALY). Deterministic and probabilistic sensitivity analyses investigated the impact of varying ERS cost and effectiveness assumptions. Sub-group analyses explored the cost-effectiveness of ERS in sedentary people with an underlying condition. Results Compared with usual care, the mean incremental lifetime cost per patient for ERS was £169 and the mean incremental QALY was 0.008, generating a base-case incremental cost-effectiveness ratio (ICER) for ERS at £20,876 per QALY in sedentary individuals without a diagnosed medical condition. There was a 51% probability that ERS was cost-effective at £20,000 per QALY and 88% probability that ERS was cost-effective at £30,000 per QALY. In sub-group analyses, cost per QALY for ERS in sedentary obese individuals was £14,618, and in sedentary hypertensives and sedentary individuals with depression the estimated cost per QALY was £12,834 and £8,414 respectively. Incremental lifetime costs and benefits associated with ERS were small, reflecting the preventative public health context of the intervention, with this resulting in estimates of cost-effectiveness that are sensitive to variations in the relative risk of becoming physically active and cost of ERS. Conclusions ERS is associated with modest increase in lifetime costs and benefits. The cost-effectiveness of ERS is highly sensitive to small changes in the effectiveness and cost of ERS and is subject to some significant uncertainty mainly due to limitations in the clinical effectiveness evidence base.

Stochastic modelling of monthly rainfall

A transformation is suggested which can transform a non-Gaussian monthly hydrological time series into a Gaussian one. The suggested approach is verified with data of ten Indian rainfall time series. Incidentally, it is observed that once the deterministic trends are removed, the transformation leads to an uncorrelated process for monthly rainfall. The procedure for normalization is general enough in that it should be also applicable to river discharges. This is verified to a limited extent by considering data of two Indian river discharges.

A semi-empirical theory of molten salts

A semi-empirical model is presented for describing the interionic interactions in molten salts using the experimentally available structure data. An extension of Bertaut's method of non-overlapping charges is used to estimate the electrostatic interaction energy in ionic melts. It is shown, in agreement with earlier computer simulation studies, that this energy increases when an ionic salt melts. The repulsion between ions is described using a compressible ion theory which uses structure-independent parameters. The van der Waals interactions and the thermal free energy are also included in the total energy, which is minimised with respect to isostructural volume variations to calculate the equilibrium density. Detailed results are presented for three molten systems, NaCl, CaCl2 and ZnCl2, and are shown to be in satisfactory agreement with experiments. With reliable structural data now being reported for several other molten salts, the present study gains relevance.