Hierarchical modeling of entangled polymers

We demonstrate that it is possible to link multi-chain molecular dynamics simulations with the tube model using a single chain slip-links model as a bridge. This hierarchical approach allows significant speed up of simulations, permitting us to span the time scales relevant for a comparison with the tube theory. Fitting the mean-square displacement of individual monomers in molecular dynamics simulations with the slip-spring model, we show that it is possible to predict the stress relaxation. Then, we analyze the stress relaxation from slip-spring simulations in the framework of the tube theory. In the absence of constraint release, we establish that the relaxation modulus can be decomposed as the sum of contributions from fast and longitudinal Rouse modes, and tube survival. Finally, we discuss some open questions regarding possible future directions that could be profitable in rendering the tube model quantitative, even for mildly entangled polymers

Estimating the spatial scales of regionalized variables by nested sampling, hierarchical analysis of variance and residual maximum likelihood

The variogram is essential for local estimation and mapping of any variable by kriging. The variogram itself must usually be estimated from sample data. The sampling density is a compromise between precision and cost, but it must be sufficiently dense to encompass the principal spatial sources of variance. A nested, multi-stage, sampling with separating distances increasing in geometric progression from stage to stage will do that. The data may then be analyzed by a hierarchical analysis of variance to estimate the components of variance for every stage, and hence lag. By accumulating the components starting from the shortest lag one obtains a rough variogram for modest effort. For balanced designs the analysis of variance is optimal; for unbalanced ones, however, these estimators are not necessarily the best, and the analysis by residual maximum likelihood (REML) will usually be preferable. The paper summarizes the underlying theory and illustrates its application with data from three surveys, one in which the design had four stages and was balanced and two implemented with unbalanced designs to economize when there were more stages. A Fortran program is available for the analysis of variance, and code for the REML analysis is listed in the paper. (c) 2005 Elsevier Ltd. All rights reserved.

Estimating the spatial scale of herbicide and soil interactions by nested sampling, hierarchical analysis of variance and residual maximum likelihood

An unbalanced nested sampling design was used to investigate the spatial scale of soil and herbicide interactions at the field scale. A hierarchical analysis of variance based on residual maximum likelihood (REML) was used to analyse the data and provide a first estimate of the variogram. Soil samples were taken at 108 locations at a range of separating distances in a 9 ha field to explore small and medium scale spatial variation. Soil organic matter content, pH, particle size distribution, microbial biomass and the degradation and sorption of the herbicide, isoproturon, were determined for each soil sample. A large proportion of the spatial variation in isoproturon degradation and sorption occurred at sampling intervals less than 60 m, however, the sampling design did not resolve the variation present at scales greater than this. A sampling interval of 20-25 m should ensure that the main spatial structures are identified for isoproturon degradation rate and sorption without too great a loss of information in this field.

Estimating the spatial scales of regionalized variables by nested sampling, hierarchical analysis of variance and residual maximum likelihood

The variogram is essential for local estimation and mapping of any variable by kriging. The variogram itself must usually be estimated from sample data. The sampling density is a compromise between precision and cost, but it must be sufficiently dense to encompass the principal spatial sources of variance. A nested, multi-stage, sampling with separating distances increasing in geometric progression from stage to stage will do that. The data may then be analyzed by a hierarchical analysis of variance to estimate the components of variance for every stage, and hence lag. By accumulating the components starting from the shortest lag one obtains a rough variogram for modest effort. For balanced designs the analysis of variance is optimal; for unbalanced ones, however, these estimators are not necessarily the best, and the analysis by residual maximum likelihood (REML) will usually be preferable. The paper summarizes the underlying theory and illustrates its application with data from three surveys, one in which the design had four stages and was balanced and two implemented with unbalanced designs to economize when there were more stages. A Fortran program is available for the analysis of variance, and code for the REML analysis is listed in the paper. (c) 2005 Elsevier Ltd. All rights reserved.

The seasonal forecast of electricity demand: A hierarchical Bayesian model with climatological weather generator

In this paper we focus on the one year ahead prediction of the electricity peak-demand daily trajectory during the winter season in Central England and Wales. We define a Bayesian hierarchical model for predicting the winter trajectories and present results based on the past observed weather. Thanks to the flexibility of the Bayesian approach, we are able to produce the marginal posterior distributions of all the predictands of interest. This is a fundamental progress with respect to the classical methods. The results are encouraging in both skill and representation of uncertainty. Further extensions are straightforward at least in principle. The main two of those consist in conditioning the weather generator model with respect to additional information like the knowledge of the first part of the winter and/or the seasonal weather forecast. Copyright (C) 2006 John Wiley & Sons, Ltd.

A hierarchical distributed approach for mining molecular fragments

Recently, two approaches have been introduced that distribute the molecular fragment mining problem. The first approach applies a master/worker topology, the second approach, a completely distributed peer-to-peer system, solves the scalability problem due to the bottleneck at the master node. However, in many real world scenarios the participating computing nodes cannot communicate directly due to administrative policies such as security restrictions. Thus, potential computing power is not accessible to accelerate the mining run. To solve this shortcoming, this work introduces a hierarchical topology of computing resources, which distributes the management over several levels and adapts to the natural structure of those multi-domain architectures. The most important aspect is the load balancing scheme, which has been designed and optimized for the hierarchical structure. The approach allows dynamic aggregation of heterogenous computing resources and is applied to wide area network scenarios.

Clinical and subclinical mastitis in smallholder dairy farms in Tanzania: Risk, intervention and knowledge transfer

In a cross-sectional study of 400 randomly selected smallholder dairy farms in the Tanga and Iringa regions of Tanzania, 14.2% (95% confidence interval (CI) = 11.6-17.3) of cows had developed clinical mastitis during the previous year. The point prevalence of subclinical mastitis, defined as a quarter positive by the California Mastitis Test (CMT) or by bacteriological culture, was 46.2% (95% Cl = 43.6-48.8) and 24.3% (95% Cl = 22.2-26.6), respectively. In a longitudinal disease study in Iringa, the incidence of clinical mastitis was 31.7 cases per 100 cow-years. A randomised intervention trial indicated that intramammary antibiotics significantly reduced the proportion of bacteriologically positive quarters in the short-term (14 days post-infusion) but teat dipping had no detectable effect on bacteriological infection and CMT positive quarters. Other risk and protective factors were identified from both the cross-sectional and longitudinal included animals with Boran breeding (odds ratio (OR) = 3,40, 95% CI = 1.00-11.57, P < 0.05 for clinical mastitis, and OR = 3.51, 95% CI = 1.299.55, P < 0.01 for a CMT positive quarter), while the practice of residual calf suckling was protective for a bacteriologically positive quarter (OR = 0.63, 95% Cl = 0.48-0.81, P <= 0.001) and for a CMT positive quarter (OR = 0.69, 95% Cl = 0.63-0.75, P < 0.001). A mastitis training course for farmers and extension officers was held, and the knowledge gained and use of different methods of dissemination were assessed over time. In a subsequent randomised controlled trial, there were strong associations between knowledge gained and both the individual question asked and the combination of dissemination methods (village meeting, video and handout) used. This study demonstrated that both clinical and subclinical mastitis is common in smallholder dairying in Tanzania, and that some of the risk and protective factors for mastitis can be addressed by practical management of dairy cows following effective knowledge transfer. (c) 2006 Elsevier B.V. All rights reserved.

Clinical and subclinical mastitis in smallholder dairy farms in Tanzania: Risk, intervention and knowledge transfer

In a cross-sectional study of 400 randomly selected smallholder dairy farms in the Tanga and Iringa regions of Tanzania, 14.2% (95% confidence interval (CI) = 11.6-17.3) of cows had developed clinical mastitis during the previous year. The point prevalence of subclinical mastitis, defined as a quarter positive by the California Mastitis Test (CMT) or by bacteriological culture, was 46.2% (95% Cl = 43.6-48.8) and 24.3% (95% Cl = 22.2-26.6), respectively. In a longitudinal disease study in Iringa, the incidence of clinical mastitis was 31.7 cases per 100 cow-years. A randomised intervention trial indicated that intramammary antibiotics significantly reduced the proportion of bacteriologically positive quarters in the short-term (14 days post-infusion) but teat dipping had no detectable effect on bacteriological infection and CMT positive quarters. Other risk and protective factors were identified from both the cross-sectional and longitudinal included animals with Boran breeding (odds ratio (OR) = 3,40, 95% CI = 1.00-11.57, P < 0.05 for clinical mastitis, and OR = 3.51, 95% CI = 1.299.55, P < 0.01 for a CMT positive quarter), while the practice of residual calf suckling was protective for a bacteriologically positive quarter (OR = 0.63, 95% Cl = 0.48-0.81, P <= 0.001) and for a CMT positive quarter (OR = 0.69, 95% Cl = 0.63-0.75, P < 0.001). A mastitis training course for farmers and extension officers was held, and the knowledge gained and use of different methods of dissemination were assessed over time. In a subsequent randomised controlled trial, there were strong associations between knowledge gained and both the individual question asked and the combination of dissemination methods (village meeting, video and handout) used. This study demonstrated that both clinical and subclinical mastitis is common in smallholder dairying in Tanzania, and that some of the risk and protective factors for mastitis can be addressed by practical management of dairy cows following effective knowledge transfer. (c) 2006 Elsevier B.V. All rights reserved.

Phylogeny and the hierarchical organization of plant diversity

R. H. Whittaker's idea that plant diversity can be divided into a hierarchy of spatial components from alpha at the within-habitat scale through beta for the turnover of species between habitats to gamma along regional gradients implies the underlying existence of alpha, beta, and gamma niches. We explore the hypothesis that the evolution of a, (3, and gamma niches is also hierarchical, with traits that define the a niche being labile, while those defining a and 7 niches are conservative. At the a level we find support for the hypothesis in the lack of close significant phylogenetic relationship between meadow species that have similar a niches. In a second test, a niche overlap based on a variety of traits is compared between congeners and noncongeners in several communities; here, too, there is no evidence of a correlation between a niche and phylogeny. To test whether beta and gamma niches evolve conservatively, we reconstructed the evolution of relevant traits on evolutionary trees for 14 different clades. Tests against null models revealed a number of instances, including some in island radiations, in which habitat (beta niche) and elevational maximum (an aspect of the gamma niche) showed evolutionary conservatism.

A hierarchical Bayesian model for predicting the functional consequences of amino-acid polymorphisms

Genetic polymorphisms in deoxyribonucleic acid coding regions may have a phenotypic effect on the carrier, e.g. by influencing susceptibility to disease. Detection of deleterious mutations via association studies is hampered by the large number of candidate sites; therefore methods are needed to narrow down the search to the most promising sites. For this, a possible approach is to use structural and sequence-based information of the encoded protein to predict whether a mutation at a particular site is likely to disrupt the functionality of the protein itself. We propose a hierarchical Bayesian multivariate adaptive regression spline (BMARS) model for supervised learning in this context and assess its predictive performance by using data from mutagenesis experiments on lac repressor and lysozyme proteins. In these experiments, about 12 amino-acid substitutions were performed at each native amino-acid position and the effect on protein functionality was assessed. The training data thus consist of repeated observations at each position, which the hierarchical framework is needed to account for. The model is trained on the lac repressor data and tested on the lysozyme mutations and vice versa. In particular, we show that the hierarchical BMARS model, by allowing for the clustered nature of the data, yields lower out-of-sample misclassification rates compared with both a BMARS and a frequen-tist MARS model, a support vector machine classifier and an optimally pruned classification tree.

The seasonal forecast of electricity demand: A hierarchical Bayesian model with climatological weather generator

In this paper we focus on the one year ahead prediction of the electricity peak-demand daily trajectory during the winter season in Central England and Wales. We define a Bayesian hierarchical model for predicting the winter trajectories and present results based on the past observed weather. Thanks to the flexibility of the Bayesian approach, we are able to produce the marginal posterior distributions of all the predictands of interest. This is a fundamental progress with respect to the classical methods. The results are encouraging in both skill and representation of uncertainty. Further extensions are straightforward at least in principle. The main two of those consist in conditioning the weather generator model with respect to additional information like the knowledge of the first part of the winter and/or the seasonal weather forecast. Copyright (C) 2006 John Wiley & Sons, Ltd.

Achievement of dietary fatty acid intakes in long-term controlled intervention studies: approach and methodology

Objective: To describe the calculations and approaches used to design experimental diets of differing saturated fatty acid (SFA) and monounsaturated fatty acid (MUFA) compositions for use in a long-term dietary intervention study, and to evaluate the degree to which the dietary targets were met. Design, setting and subjects: Fifty-one students living in a university hall of residence consumed a reference (SFA) diet for 8 weeks followed by either a moderate MUFA (MM) diet or a high MUFA (HM) diet for 16 weeks. The three diets were designed to differ only in their proportions of SFA and MUFA, while keeping total fat, polyunsaturated fatty acids (PUFA), trans-fatty acids, and the ratio of palmitic to stearic acid, and n-6 to n-3 PUFA, unchanged. Results: Using habitual diet records and a standardised database for food fatty acid compositions, a sequential process of theoretical fat substitutions enabled suitable fat sources for use in the three diets to be identified, and experimental margarines for baking, spreading and the manufacture of snack foods to be designed. The dietary intervention was largely successful in achieving the fatty acid targets of the three diets, although unintended differences between the original target and the analysed fatty acid composition of the experimental margarines resulted in a lower than anticipated MUFA intake on the HM diet, and a lower ratio of palmitic to stearic acid compared with the reference or MM diet. Conclusions: This study has revealed important theoretical considerations that should be taken into account when designing diets of specific fatty acid composition, as well as practical issues of implementation.

Fecal water as a non-invasive biomarker in nutritional intervention: comparison of preparation methods and refinement of different endpoints

The assessment of cellular effects by the aqueous phase of human feces (fecal water, FW) is a useful biomarker approach to study cancer risks and protective activities of food. In order to refine and develop the biomarker, different protocols of preparing FW were compared. Fecal waters were prepared by 3 methods: (A) direct centrifugation; (B) extraction of feces in PBS before centrifugation; and (C) centrifugation of lyophilized and reconstituted feces. Genotoxicity was determined in colon cells using the Comet assay. Selected samples were investigated for additional parameters related to carcinogenesis. Two of 7 FWs obtained by methods A and B were similarly genotoxic. Method B, however, yielded higher volumes of FW, allowing sterile filtration for long-term culture experiments. Four of 7 samples were non-genotoxic when prepared according to all 3 methods. FW from lyophilized feces and from fresh samples were equally genotoxic. FWs modulated cytotoxicity, paracellular permeability, and invasion, independent of their genotoxicity. All 3 methods of FW preparation can be used to assess genotoxicity. The higher volumes of FWobtained by preparation method B greatly enhance the perspectives of measuring different types of biological parameters and using these to disclose activities related to cancer development.