Tropospheric planetary wave dynamics and mixture modeling: Two preferred regimes and a regime shift

Investigation of preferred structures of planetary wave dynamics is addressed using multivariate Gaussian mixture models. The number of components in the mixture is obtained using order statistics of the mixing proportions, hence avoiding previous difficulties related to sample sizes and independence issues. The method is first applied to a few low-order stochastic dynamical systems and data from a general circulation model. The method is next applied to winter daily 500-hPa heights from 1949 to 2003 over the Northern Hemisphere. A spatial clustering algorithm is first applied to the leading two principal components (PCs) and shows significant clustering. The clustering is particularly robust for the first half of the record and less for the second half. The mixture model is then used to identify the clusters. Two highly significant extratropical planetary-scale preferred structures are obtained within the first two to four EOF state space. The first pattern shows a Pacific-North American (PNA) pattern and a negative North Atlantic Oscillation (NAO), and the second pattern is nearly opposite to the first one. It is also observed that some subspaces show multivariate Gaussianity, compatible with linearity, whereas others show multivariate non-Gaussianity. The same analysis is also applied to two subperiods, before and after 1978, and shows a similar regime behavior, with a slight stronger support for the first subperiod. In addition a significant regime shift is also observed between the two periods as well as a change in the shape of the distribution. The patterns associated with the regime shifts reflect essentially a PNA pattern and an NAO pattern consistent with the observed global warming effect on climate and the observed shift in sea surface temperature around the mid-1970s.

Biosensor-based diagnostics of contaminated groundwater: assessment and remediation strategy

Shallow groundwater beneath a former airfield site in southern England has been heavily contaminated with a wide range of chlorinated solvents. The feasibility of using bacterial biosensors to complement chemical analysis and enable cost-effective, and focussed sampling has been assessed as part of a site evaluation programme. Five different biosensors, three metabolic (Vibrio fischeri, Pseudomonas fluorescens 10568 and Escherichia coli HB101) and two catabolic (Pseudomonas putida TVA8 and E. coli DH5alpha), were employed to identify areas where the availability and toxicity of pollutants is of most immediate environmental concern. The biosensors used showed different sensitivities to each other and to the groundwater samples tested. There was generally a good agreement with chemical analyses. The potential efficacy of remediation strategies was explored by coupling sample manipulation to biosensor tests. Manipulation involved sparging and charcoal treatment procedures to simulate remediative engineering solutions. Sparging was sufficient at most locations. (C) 2004 Elsevier Ltd. All rights reserved.

The development and contribution of surface charge by crop residues in two Malawian acid soils

The effects of maize and soya bean residues on the pH and charge of a loamy sand (Kawalazi) and a sandy clay loam (Naming'omba) from Malawi were measured to determine both the indirect effect of the residues on soil charge through the changes in pH, and the direct contribution of charge carried on the residue surfaces. The soils had pH values (10 mM CaCl2) of 4.3 and 5.0 and organic matter contents were 1.4% and 2.7%, respectively. The clay fractions were dominated by kaolinite and goethite, and mica was present in both samples. The soils were incubated for 28 days with maize (Zea mays) and soya bean (Glycine max) residues. The maximum addition of residue (12.0%) in the Kawalazi and Naming'omba soils increased the pH from 4.3 and 5.0 to 4.8 and 5.3 (maize) and to 9.0 and 8.8 (soya bean), respectively. Negative charge increased from 2.1 and 4.7 cmol(c) kg(-1) to 3.8 and 7.5 (maize) and to 5.3 and 9.3 cmol(c) kg(-1) (soya bean). Positive charge increased from 0.72 and 0.62 to 0.87 and 0.85 cmol(c) kg(-1) (maize) and to 0.75 and 0.68 (soya bean). The charge contribution by the residues was calculated by difference between the charge on a sample incubated with residue and the charge on a soil without residue limed to the same pH value. Up to 100 cmolc negative charge and 10 cmol(c) of positive charge per kg of residue were directly contributed to the soil-residue mixture, the amounts depending on the type of residue, the extent to which the residue was decomposed in the soil and the pH of the mixture. The Anderson and Sposito method [Soil Sci. Soc. Am. J. 55 (1991) 1569] was used to partition the permanent negative charge (holding Cs+) from variable negative charge (holding Li+). In the pH range 3.7-6.5 the maize residue contributed between 3 and 26 cmol(c) of variable charge per kg of residue in the Kawalazi soil and between 6 and 25 cmol(c) per kg of residue in the Naming'omba soil. For soya bean the values were between I and 28 and between 4 and 68 cmolc per kg of residue, respectively. At a given pH value, the charge tended to increase with time of incubation and for a given addition of residue, pH decreased during incubation. Addition of residues contributed no permanent negative charge and the charge on the soil measured by Cs adsorption was independent of pH change caused by the residue showing that the method is valid for soil-residue mixtures. With time there was a decrease in the amount of permanent charge probably due to masking as humic material become adsorbed on mineral surfaces. (C) 2003 Elsevier Science B.V. All rights reserved.

A comparison of two methods for developing the linearization of a shallow-water model

We develop the linearization of a semi-implicit semi-Lagrangian model of the one-dimensional shallow-water equations using two different methods. The usual tangent linear model, formed by linearizing the discrete nonlinear model, is compared with a model formed by first linearizing the continuous nonlinear equations and then discretizing. Both models are shown to perform equally well for finite perturbations. However, the asymptotic behaviour of the two models differs as the perturbation size is reduced. This leads to difficulties in showing that the models are correctly coded using the standard tests. To overcome this difficulty we propose a new method for testing linear models, which we demonstrate both theoretically and numerically. © Crown copyright, 2003. Royal Meteorological Society

The effect of consumption volume on profile and liking of oral nutritional supplements of varied sweetness: sequential profiling and boredom tests

Oral nutrition supplements (ONS) are routinely prescribed to those with, or at risk of, malnutrition. Previous research identified poor compliance due to taste and sweetness. This paper investigates taste and hedonic liking of ONS, of varying sweetness and metallic levels, over consumption volume; an important consideration as patients are prescribed large volumes of ONS daily. A sequential descriptive profile was developed to determine the perception of sensory attributes over repeat consumption of ONS. Changes in liking of ONS following repeat consumption were characterised by a boredom test. Certain flavour (metallic taste, soya milk flavour) and mouthfeel (mouthdrying, mouthcoating) attributes built up over increased consumption volume (p 0.002). Hedonic liking data from two cohorts, healthy older volunteers (n = 32, median age 73) and patients (n = 28, median age 85), suggested such build-up was disliked. Efforts made to improve the palatability of ONS must take account of the build up of taste and mouthfeel characteristics over increased consumption volume.

Phyto-activity of biocides used to manipulate herbivory: tests of three pesticides on fourteen plant species

1. We tested three pesticides used for field manipulations of herbivory for direct phytoactive effects on the germination and growth of 14 herbaceous plant species selected to provide a range of life-history strategies and functional groups. 2. We report three companion experiments: (A) Two insecticides, chlorpyrifos (granular soil insecticide) and dimethoate (foliar spray), were applied in fully-factorial combination to pot-germinated individuals of 12 species. (B) The same fully-factorial design was used to test for direct effects on the germination of four herbaceous legumes. (C) The molluscicide, metaldehyde, was tested for direct effects on the germination and growth of six plant species. 3. The insecticides had few significant effects on growth and germination. Dimethoate acted only on growth stimulating Anisantha sterilis, Sonchus asper and Stellaria graminea. In contrast, chlorpyrifos acted on germination increasing the germination of Trifolium dubium and Trifolium pratense. There was also a significant interactive effect of chlorpyrifos and dimethoate on the germination of T pratense. However, all. effects were relatively small in magnitude and explanatory power. The molluscicide had no significant effect on plant germination or growth. 4. The small number and size of direct effects of the pesticides on plant performance is encouraging for the use of these pesticides in manipulative experiments on herbivory, especially for the molluscicide. However, a smatt number of direct (positive) effects of the insecticides on some plant species need to be taken into account when interpreting field manipulations of herbivory with these compounds, and emphasises the importance of conducting tests for direct phyto-active effects. (C) 2004 Elsevier GmbH. All rights reserved.

Comparison of sample size formulae for 2x2 cross-over designs applied to bioequivalence studies

We consider the comparison of two formulations in terms of average bioequivalence using the 2 × 2 cross-over design. In a bioequivalence study, the primary outcome is a pharmacokinetic measure, such as the area under the plasma concentration by time curve, which is usually assumed to have a lognormal distribution. The criterion typically used for claiming bioequivalence is that the 90% confidence interval for the ratio of the means should lie within the interval (0.80, 1.25), or equivalently the 90% confidence interval for the differences in the means on the natural log scale should be within the interval (-0.2231, 0.2231). We compare the gold standard method for calculation of the sample size based on the non-central t distribution with those based on the central t and normal distributions. In practice, the differences between the various approaches are likely to be small. Further approximations to the power function are sometimes used to simplify the calculations. These approximations should be used with caution, because the sample size required for a desirable level of power might be under- or overestimated compared to the gold standard method. However, in some situations the approximate methods produce very similar sample sizes to the gold standard method. Copyright © 2005 John Wiley & Sons, Ltd.

Comparing correlations of continuous observations from two independent populations using a sequential approach

A sequential study design generally makes more efficient use of available information than a fixed sample counterpart of equal power. This feature is gradually being exploited by researchers in genetic and epidemiological investigations that utilize banked biological resources and in studies where time, cost and ethics are prominent considerations. Recent work in this area has focussed on the sequential analysis of matched case-control studies with a dichotomous trait. In this paper, we extend the sequential approach to a comparison of the associations within two independent groups of paired continuous observations. Such a comparison is particularly relevant in familial studies of phenotypic correlation using twins. We develop a sequential twin method based on the intraclass correlation and show that use of sequential methodology can lead to a substantial reduction in the number of observations without compromising the study error rates. Additionally, our approach permits straightforward allowance for other explanatory factors in the analysis. We illustrate our method in a sequential heritability study of dysplasia that allows for the effect of body mass index and compares monozygotes with pairs of singleton sisters. Copyright (c) 2006 John Wiley & Sons, Ltd.

A simple two-stage design for quantitative responses with application to a study in diabetic neuropathic pain

Two-stage designs offer substantial advantages for early phase II studies. The interim analysis following the first stage allows the study to he stopped for futility, or more positively, it might lead to early progression to the trials needed for late phase H and phase III. If the study is to continue to its second stage, then there is an opportunity for a revision of the total sample size. Two-stage designs have been implemented widely in oncology studies in which there is a single treatment arm and patient responses are binary. In this paper the case of two-arm comparative studies in which responses are quantitative is considered. This setting is common in therapeutic areas other than oncology. It will be assumed that observations are normally distributed, but that there is some doubt concerning their standard deviation, motivating the need for sample size review. The work reported has been motivated by a study in diabetic neuropathic pain, and the development of the design for that trial is described in detail. Copyright (C) 2008 John Wiley & Sons, Ltd.

Bayesian sample size for exploratory clinical trials incorporating historical data

This paper presents a simple Bayesian approach to sample size determination in clinical trials. It is required that the trial should be large enough to ensure that the data collected will provide convincing evidence either that an experimental treatment is better than a control or that it fails to improve upon control by some clinically relevant difference. The method resembles standard frequentist formulations of the problem, and indeed in certain circumstances involving 'non-informative' prior information it leads to identical answers. In particular, unlike many Bayesian approaches to sample size determination, use is made of an alternative hypothesis that an experimental treatment is better than a control treatment by some specified magnitude. The approach is introduced in the context of testing whether a single stream of binary observations are consistent with a given success rate p(0). Next the case of comparing two independent streams of normally distributed responses is considered, first under the assumption that their common variance is known and then for unknown variance. Finally, the more general situation in which a large sample is to be collected and analysed according to the asymptotic properties of the score statistic is explored. Copyright (C) 2007 John Wiley & Sons, Ltd.

Detection of frailty in Weibull lifetime data using Outlier tests

Heterogeneity in lifetime data may be modelled by multiplying an individual's hazard by an unobserved frailty. We test for the presence of frailty of this kind in univariate and bivariate data with Weibull distributed lifetimes, using statistics based on the ordered Cox-Snell residuals from the null model of no frailty. The form of the statistics is suggested by outlier testing in the gamma distribution. We find through simulation that the sum of the k largest or k smallest order statistics, for suitably chosen k , provides a powerful test when the frailty distribution is assumed to be gamma or positive stable, respectively. We provide recommended values of k for sample sizes up to 100 and simple formulae for estimated critical values for tests at the 5% level.

Melt structure and its transformation by sequential crystallization of the two blocks within poly(L-lactide)-block-poly(epsilon-caprolactone) double crystalline diblock copolymers

Sequential crystallization of poly(L-lactide) (PLLA) followed by poly(epsilon-caprolactone) (PCL) in double crystalline PLLA-b-PCL diblock copolymers is studied by differential scanning calorimetry (DSC), polarized optical microscopy (POM), wide-angle X-ray scattering (WAXS) and small-angle X-ray scattering (SAXS). Three samples with different compositions are studied. The sample with the shortest PLLA block (32 wt.-% PLLA) crystallizes from a homogeneous melt, the other two (with 44 and 60% PLLA) from microphase separated structures. The microphase structure of the melt is changed as PLLA crystallizes at 122 degrees C (a temperature at which the PCL block is molten) forming spherulites regardless of composition, even with 32% PLLA. SAXS indicates that a lamellar structure with a different periodicity than that obtained in the melt forms (for melt segregated samples). Where PCL is the majority block, PCL crystallization at 42 degrees C following PLLA crystallization leads to rearrangement of the lamellar structure, as observed by SAXS, possibly due to local melting at the interphases between domains. POM results showed that PCL crystallizes within previously formed PLLA spherulites. WAXS data indicate that the PLLA unit cell is modified by crystallization of PCL, at least for the two majority PCL samples. The PCL minority sample did not crystallize at 42 degrees C (well below the PCL homopolymer crystallization temperature), pointing to the influence of pre-crystallization of PLLA on PCL crystallization, although it did crystallize at lower temperature. Crystallization kinetics were examined by DSC and WAXS, with good agreement in general. The crystallization rate of PLLA decreased with increase in PCL content in the copolymers. The crystallization rate of PCL decreased with increasing PLLA content. The Avrami exponents were in general depressed for both components in the block copolymers compared to the parent homopolymers. Polarized optical micrographs during isothermal crystalli zation of (a) homo-PLLA, (b) homo-PCL, (c) and (d) block copolymer after 30 min at 122 degrees C and after 15 min at 42 degrees C.

A comparison of two time-domain analysis procedures in the determination of VO2 kinetics by pseudorandom binary sequence exercise testing

The purpose of this study was to apply and compare two time-domain analysis procedures in the determination of oxygen uptake (VO2) kinetics in response to a pseudorandom binary sequence (PRBS) exercise test. PRBS exercise tests have typically been analysed in the frequency domain. However, the complex interpretation of frequency responses may have limited the application of this procedure in both sporting and clinical contexts, where a single time measurement would facilitate subject comparison. The relative potential of both a mean response time (MRT) and a peak cross-correlation time (PCCT) was investigated. This study was divided into two parts: a test-retest reliability study (part A), in which 10 healthy male subjects completed two identical PRBS exercise tests, and a comparison of the VO2 kinetics of 12 elite endurance runners (ER) and 12 elite sprinters (SR; part B). In part A, 95% limits of agreement were calculated for comparison between MRT and PCCT. The results of part A showed no significant difference between test and retest as assessed by MRT [mean (SD) 42.2 (4.2) s and 43.8 (6.9) s] or by PCCT [21.8 (3.7) s and 22.7 (4.5) s]. Measurement error (%) was lower for MRT in comparison with PCCT (16% and 25%, respectively). In part B of the study, the VO2 kinetics of ER were significantly faster than those of SR, as assessed by MRT [33.4 (3.4) s and 39.9 (7.1) s, respectively; P<0.01] and PCCT [20.9 (3.8) s and 24.8 (4.5) s; P < 0.05]. It is possible that either analysis procedure could provide a single test measurement Of VO2 kinetics; however, the greater reliability of the MRT data suggests that this method has more potential for development in the assessment Of VO2 kinetics by PRBS exercise testing.

PASSCLAIM - gut health and immunity

Background The gut and immune system form a complex integrated structure that has evolved to provide effective digestion and defence against ingested toxins and pathogenic bacteria. However, great variation exists in what is considered normal healthy gut and immune function. Thus, whilst it is possible to measure many aspects of digestion and immunity, it is more difficult to interpret the benefits to individuals of variation within what is considered to be a normal range. Nevertheless, it is important to set standards for optimal function for use both by the consumer, industry and those concerned with the public health. The digestive tract is most frequently the object of functional and health claims and a large market already exists for gut-functional foods worldwide. Aim To define normal function of the gut and immune system and describe available methods of measuring it. Results We have defined normal bowel habit and transit time, identified their role as risk factors for disease and how they may be measured. Similarly, we have tried to define what is a healthy gut flora in terms of the dominant genera and their metabolism and listed the many, varied and novel methods for determining these parameters. It has proved less easy to provide boundaries for what constitutes optimal or improved gastric emptying, gut motility, nutrient and water absorption and the function of organs such as the liver, gallbladder and pancreas. The many tests of these functions are described. We have discussed gastrointestinal well being. Sensations arising from the gut can be both pleasant and unpleasant. However, the characteristics of well being are ill defined and merge imperceptibly from acceptable to unacceptable, a state that is subjective. Nevertheless, we feel this is an important area for future work and method development. The immune system is even more difficult to make quantitative judgements about. When it is defective, then clinical problems ensure, but this is an uncommon state. The innate and adaptive immune systems work synergistically together and comprise many cellular and humoral factors. The adaptive system is extremely sophisticated and between the two arms of immunity there is great redundancy, which provides robust defences. New aspects of immune function are discovered regularly. It is not clear whether immune function can be "improved". Measuring aspects of immune function is possible but there is no one test that will define either the status or functional capacity of the immune system. Human studies are often limited by the ability to sample only blood or secretions such as saliva but it should be remembered that only 2% of lymphocytes circulate at any given time, which limits interpretation of data. We recommend assessing the functional capacity of the immune system by: measuring specific cell functions ex vivo, measuring in vivo responses to challenge, e. g. change in antibody in blood or response to antigens, determining the incidence and severity of infection in target populations during naturally occurring episodes or in response to attenuated pathogens.

Expressive versus receptive language skills in specific reading disorder

Despite ample research into the language skills of children with specific reading disorder no studies so far have investigated whether there may be a difference between expressive and receptive language skills in this population. Yet, neuro-anatomical models would predict that children who have specific reading disorder which is not associated with movement or attention difficulties, would have lower receptive language skills than expressive. This study investigates the difference between expressive and receptive language skills in a sample of 17 children with specific reading difficulty aged between 7 and 12 years. They were administered a battery of two receptive and two expressive language measures. The results showed that as the neuro-anatomical model would predict, the children scored significantly lower on tests of receptive than on tests of expressive language skills.