Parametric accelerated failure time models with random effects and an application to kidney transplant survival

Accelerated failure time models with a shared random component are described, and are used to evaluate the effect of explanatory factors and different transplant centres on survival times following kidney transplantation. Different combinations of the distribution of the random effects and baseline hazard function are considered and the fit of such models to the transplant data is critically assessed. A mixture model that combines short- and long-term components of a hazard function is then developed, which provides a more flexible model for the hazard function. The model can incorporate different explanatory variables and random effects in each component. The model is straightforward to fit using standard statistical software, and is shown to be a good fit to the transplant data. Copyright (C) 2004 John Wiley Sons, Ltd.

Demographic models and the management of endangered species: a case study of the critically endangered Seychelles magpie robin

1. Demographic models are assuming an important role in management decisions for endangered species. Elasticity analysis and scope for management analysis are two such applications. Elasticity analysis determines the vital rates that have the greatest impact on population growth. Scope for management analysis examines the effects that feasible management might have on vital rates and population growth. Both methods target management in an attempt to maximize population growth. 2. The Seychelles magpie robin Copsychus sechellarum is a critically endangered island endemic, the population of which underwent significant growth in the early 1990s following the implementation of a recovery programme. We examined how the formal use of elasticity and scope for management analyses might have shaped management in the recovery programme, and assessed their effectiveness by comparison with the actual population growth achieved. 3. The magpie robin population doubled from about 25 birds in 1990 to more than 50 by 1995. A simple two-stage demographic model showed that this growth was driven primarily by a significant increase in the annual survival probability of first-year birds and an increase in the birth rate. Neither the annual survival probability of adults nor the probability of a female breeding at age 1 changed significantly over time. 4. Elasticity analysis showed that the annual survival probability of adults had the greatest impact on population growth. There was some scope to use management to increase survival, but because survival rates were already high (> 0.9) this had a negligible effect on population growth. Scope for management analysis showed that significant population growth could have been achieved by targeting management measures at the birth rate and survival probability of first-year birds, although predicted growth rates were lower than those achieved by the recovery programme when all management measures were in place (i.e. 1992-95). 5. Synthesis and applications. We argue that scope for management analysis can provide a useful basis for management but will inevitably be limited to some extent by a lack of data, as our study shows. This means that identifying perceived ecological problems and designing management to alleviate them must be an important component of endangered species management. The corollary of this is that it will not be possible or wise to consider only management options for which there is a demonstrable ecological benefit. Given these constraints, we see little role for elasticity analysis because, when data are available, a scope for management analysis will always be of greater practical value and, when data are lacking, precautionary management demands that as many perceived ecological problems as possible are tackled.

Evaluation of data for atmospheric models: master equation/RRKM calculations on the combination reaction, BrO+NO2 -> BrONO2, a conundrum

Experimental data for the title reaction were modeled using master equation (ME)/RRKM methods based on the Multiwell suite of programs. The starting point for the exercise was the empirical fitting provided by the NASA (Sander, S. P.; Finlayson-Pitts, B. J.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Molina, M. J.; Moortgat, G. K.; Orkin, V. L.; Ravishankara, A. R. Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15; Jet Propulsion Laboratory: Pasadena, California, 2006)(1) and IUPAC (Atkinson, R.; Baulch, D. L.; Cox, R. A.: R. F. Hampson, J.; Kerr, J. A.; Rossi, M. J.; Troe, J. J. Phys. Chem. Ref. Data. 2000, 29, 167) 2 data evaluation panels, which represents the data in the experimental pressure ranges rather well. Despite the availability of quite reliable parameters for these calculations (molecular vibrational frequencies (Parthiban, S.; Lee, T. J. J. Chem. Phys. 2000, 113, 145)3 and a. value (Orlando, J. J.; Tyndall, G. S. J. Phys. Chem. 1996, 100,. 19398)4 of the bond dissociation energy, D-298(BrO-NO2) = 118 kJ mol(-1), corresponding to Delta H-0(circle) = 114.3 kJ mol(-1) at 0 K) and the use of RRKM/ME methods, fitting calculations to the reported data or the empirical equations was anything but straightforward. Using these molecular parameters resulted in a discrepancy between the calculations and the database of rate constants of a factor of ca. 4 at, or close to, the low-pressure limit. Agreement between calculation and experiment could be achieved in two ways, either by increasing Delta H-0(circle) to an unrealistically high value (149.3 kJ mol(-1)) or by increasing , the average energy transferred in a downward collision, to an unusually large value (> 5000 cm(-1)). The discrepancy could also be reduced by making all overall rotations fully active. The system was relatively insensitive to changing the moments of inertia in the transition state to increase the centrifugal effect. The possibility of involvement of BrOONO was tested and cannot account for the difficulties of fitting the data.

In vitro fermentation of mixed linkage glucooligosaccharides produced by Gluconobacter oxydans NCIMB 4943 by the human colonic microflora

The aim of this study was to develop selectively fermented (prebiotic) carbohydrate molecules which would also result in the generation of butyric acid. Glucooligosaccharides produced by Gluconobacter oxydans NCIMB 4943 from various types of maltodextrins were evaluated for their fermentation by mixed cultures of human colonic microflora. The selectivity of growth of desirable bacteria (bifidobacteria, lactobacilli) was studied in stirred pH-controlled (6.8) batch cultures. Bacterial populations were enumerated using fluorescent in situ hybridization (FISH). Gluco-oligosaccharides resulted in significantly (P<0.05) increased numbers of bifidobacteria and lactobacilli within 24 hours. Bacteroides, clostridial and eubacterial populations were slightly decreased at 48 h. There was very little difference in selectivity between the maltodextrin substrates and the products, although maltodextrin displayed a slightly less selective fermentation than the gluco-oligosaccharide products, also stimulating the growth of bacteroides, clostridia and eubacteria. Gluco-oligosaccharides, produced from G19 maltodextrin, resulted in the best prebiotic effect with the highest prebiotic index (PI) of 5.90 at 48 hours. Acetate, propionate and butyrate were all produced from glucooligosaccharides, derived from G19 maltodextrin, at 48 hours but no lactate or formate were detected.

Effect of forage type and proportion of concentrate in the diet on milk fatty acid composition in cows given sunflower oil and fish oil

Based on the potential benefits of cis-9, trans- 11 conjugated linoleic acid (CLA) for human health there is a need to develop effective strategies for enhancing milk fat CLA concentrations. In this experiment, the effect of forage type and level of concentrate in the diet on milk fatty acid composition was examined in cows given a mixture of fish oil and sunflower oil. Four late lactation Holstein-British Friesian cows were used in a 4 x 4 Latin-square experiment with a 2 x 2 factorial arrangement of treatments and 21-day experimental periods. Treatments consisted of grass (G) or maize (M) silage supplemented with low (L) or high (H) levels of concentrates (65: 35 and 35: 65; forage: concentrate ratio, on a dry matter (DM) basis, respectively) offered as a total mixed ration at a restricted level of intake (20 kg DM per day). Lipid supplements (30 g/kg DM) containing fish oil and sunflower oil (2: 3 w/w) were offered during the last 14 days of each experimental period. Treatments had no effect on total DM intake, milk yield, milk constituent output or milk fat content, but milk protein concentrations were lower (P<0.05) for G than M diets (mean 43.0 and 47.3 g/kg, respectively). Compared with grass silage, milk fat contained higher (P<0.05) amounts Of C-12: 0, C-14: 0, trans C-18:1 and long chain >= C20 (n-3) polyunsaturated fatty acids (PUFA) and lower (P<0.05) levels Of C-18:0 and trans C-18:2 when maize silage was offered. Increases in the proportion of concentrate in the diet elevated (P<0.05) C-18:2 (n-6) and long chain >= C20 (n-3) PUFA content, but reduced (P<0.05) the amount Of C-18:3 (n-3). Concentrations of trans-11 C-18:1 in milk were independent of forage type, but tended (P<0.10) to be lower for high concentrate diets (mean 7.2 and 4.0 g/100 g fatty acids, for L and H respectively). Concentrations of trans-10 C-18:1 were higher (P<0.05) in milk from maize compared with grass silage (mean 10.3 and 4.1 g/100 g fatty acids, respectively) and increased in response to high levels of concentrates in the diet (mean 4.1 and 10.3 g/100 g fatty acids, for L and H, respectively). Forage type had no effect (P>0.05) on total milk conjugated linoleic acid (CLA) (2.7 and 2.8 g/100 g fatty acids, for M and G, respectively) or cis-9, trans-11 CLA content (2.2 and 2.4 g/100 g fatty acids). Feeding high concentrate diets tended (P<0.10) to decrease total CLA (3.3 and 2.2 g/100 g fatty acids, for L and H, respectively) and cis-9, trans-11 CLA (2.9 and 1/7 g/100 g fatty acids) concentrations and increase milk trans-9, cis-11 CLA and trans-10, cis-12 CLA content. In conclusion, the basal diet is an important determinant of milk fatty acid composition when a supplement of fish oil and sunflower oil is given.

Effect of storage temperature and media composition on the survivability of Bifidobacterium breve NCIMB 702257 in a malt hydrolisate

The aim of this study was to evaluate the survivability of Bifidobacterium breve NCIMB 702257 in a three malt-based media supplemented with cysteine and yeast extract, and to determine the protective effect of these growth factors. A number of parameterised mathematical models were used to predict of kinetics of viability and total acidity during storage at different temperatures. Results demonstrated a good fit to the experimental mathematical model. The Arrhenius equations showed only reasonable fits and the polynomial plots contained a large area without data between 4 and 25 degrees C. In addition, it was shown that cysteine promotes growth and acid production by bifidobacteria, but does not extend survivability. On the other hand, increasing the yeast extract content of the fermentation media enhances the survivability of B. breve. To our knowledge, this is the first study to address the modelling of the survivability of probiotic bacteria in a cereal based fermentation media at different temperatures, introducing a more quantitative approach to the study of the shelf-life of a probiotic product. (C) 2009 Elsevier B.V. All rights reserved.

Fermentation of resistant starches: influence of in vitro models on colon carcinogenesis

Resistant starch type 2 (RS2) and type 3 (RS3) containing preparations were digested using a batch (a) and a dynamic in vitro model (b). Furthermore, in vivo obtained indigestible fractions from ileostomy patients were used (c). Subsequently these samples were fermented with human feces with a batch and a dynamic in vitro method. The fermentation supernatants were used to treat CAC02 cells. Cytotoxicity, anti-genotoxicity against hydrogen peroxide (comet assay) and the effect on barrier function measured by trans-epithelial electrical resistance were determine. Dynamically fermented samples led to high cytotoxic activity, probably due to additional compounds added during in vitro fermentation. As a consequence only batch fermented samples were investigated further. Batch fermentation of RS resulted in an anti-genotoxic activity ranging from 9-30% decrease in DNA damage for all the samples, except for RS2-b. It is assumed that the changes in RS2 structures due to dynamic digestion resulted in a different fermentation profile not leading to any anti-genotoxic effect. Additionally, in vitro batch fermentation of RS caused an improvement in integrity across the intestinal barrier by approximately 22% for all the samples. We have demonstrated that batch in vitro fermentation of RS2 and RS3 preparations differently pre-digested are capable of inhibiting the initiation and promotion stage in colon carcinogenesis in vitro.

Reversing the emotional Stroop effect reveals that it is not what it seems: the role of fast and slow components

The relative contributions of slow and fast (online) components in a modified emotional Stroop task were evaluated. The slow component, neglected in previous research, was shown to lead to the prediction of a reversed emotional intrusion effect using pseudorandomly mixed negative and neutral stimuli. This prediction was supported in Experiments 1 and 2. In Experiments 3 and 4, a new paradigm was developed that allowed a more direct observation of the nature of disruptive effects from negative stimuli. The results provided a clear demonstration of the presence of the slow component. The fast component, which has generally been assumed to be the source of the interference, was shown, in fact, to have little or no role in the disruption.

Comparison of subspace and ARX models of a waveguide's terahertz transient response after optimal wavelet filtering

A quasi-optical deembedding technique for characterizing waveguides is demonstrated using wide-band time-resolved terahertz spectroscopy. A transfer function representation is adopted for the description of the signal in the input and output port of the waveguides. The time-domain responses were discretized and the waveguide transfer function was obtained through a parametric approach in the z-domain after describing the system with an AutoRegressive with eXogenous input (ARX), as well as with a state-space model. Prior to the identification procedure, filtering was performed in the wavelet domain to minimize both signal distortion, as well as the noise propagating in the ARX and subspace models. The optimal filtering procedure used in the wavelet domain for the recorded time-domain signatures is described in detail. The effect of filtering prior to the identification procedures is elucidated with the aid of pole-zero diagrams. Models derived from measurements of terahertz transients in a precision WR-8 waveguide adjustable short are presented.

Effect of incremental levels of fish oil in the diet on ruminal lipid metabolism in growing steers

Based on the potential benefits to human health, there is interest in developing sustainable nutritional strategies to enhance the concentration of long-chain n-3 fatty acids in ruminant-derived foods. Four Aberdeen Angus steers fitted with rumen and duodenal cannulae were used in a 4 × 4 Latin square experiment with 21 d experimental periods to examine the potential of fish oil (FO) in the diet to enhance the supply of 20 : 5n-3 and 22 : 6n-3 available for absorption in growing cattle. Treatments consisted of total mixed rations based on maize silage fed at a rate of 85 g DM/kg live weight0·75/d containing 0, 8, 16 and 24 g FO/kg diet DM. Supplements of FO reduced linearly (P < 0·01) DM intake and shifted (P < 0·01) rumen fermentation towards propionate at the expense of acetate and butyrate. FO in the diet enhanced linearly (P < 0·05) the flow of trans-16 : 1, trans-18 : 1, trans-18 : 2, 20 : 5n-3 and 22 : 6n-3, and decreased linearly (P < 0·05) 18 : 0 and 18 : 3n-3 at the duodenum. Increases in the flow of trans-18 : 1 were isomer dependent and were determined primarily by higher amounts of trans-11 reaching the duodenum. In conclusion, FO alters ruminal lipid metabolism of growing cattle in a dose-dependent manner consistent with an inhibition of ruminal biohydrogenation, and enhances the amount of long-chain n-3 fatty acids at the duodenum, but the increases are marginal due to extensive biohydrogenation in the rumen.

Emulating AOGCM results using simple climate models

Three simple climate models (SCMs) are calibrated using simulations from atmosphere ocean general circulation models (AOGCMs). In addition to using two conventional SCMs, results from a third simpler model developed specifically for this study are obtained. An easy to implement and comprehensive iterative procedure is applied that optimises the SCM emulation of global-mean surface temperature and total ocean heat content, and, if available in the SCM, of surface temperature over land, over the ocean and in both hemispheres, and of the global-mean ocean temperature profile. The method gives best-fit estimates as well as uncertainty intervals for the different SCM parameters. For the calibration, AOGCM simulations with two different types of forcing scenarios are used: pulse forcing simulations performed with 2 AOGCMs and gradually changing forcing simulations from 15 AOGCMs obtained within the framework of the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. The method is found to work well. For all possible combinations of SCMs and AOGCMs the emulation of AOGCM results could be improved. The obtained SCM parameters depend both on the AOGCM data and the type of forcing scenario. SCMs with a poor representation of the atmosphere thermal inertia are better able to emulate AOGCM results from gradually changing forcing than from pulse forcing simulations. Correct simultaneous emulation of both atmospheric temperatures and the ocean temperature profile by the SCMs strongly depends on the representation of the temperature gradient between the atmosphere and the mixed layer. Introducing climate sensitivities that are dependent on the forcing mechanism in the SCMs allows the emulation of AOGCM responses to carbon dioxide and solar insolation forcings equally well. Also, some SCM parameters are found to be very insensitive to the fitting, and the reduction of their uncertainty through the fitting procedure is only marginal, while other parameters change considerably. The very simple SCM is found to reproduce the AOGCM results as well as the other two comparably more sophisticated SCMs.

Temperature dependence of protein dynamics simulated with three different water models

The effect of variation of the water model on the temperature dependence of protein and hydration water dynamics is examined by performing molecular dynamics simulations of myoglobin with the TIP3P, TIP4P, and TIP5P water models and the CHARMM protein force field at temperatures between 20 and 300 K. The atomic mean-square displacements, solvent reorientational relaxation times, pair angular correlations between surface water molecules, and time-averaged structures of the protein are all found to be similar, and the protein dynamical transition is described almost indistinguishably for the three water potentials. The results provide evidence that for some purposes changing the water model in protein simulations without a loss of accuracy may be possible.

Second meal effect: modified sham feeding does not provoke the release of stored triacylglycerol from a previous high-fat meal

The present study was carried out to determine whether cephalic stimulation, associated with eating a meal, was sufficient stimulus to provoke the release of stored triacylglycerol (TAG) from a previous high-fat meal. Ten subjects were studied on three separate occasions. Following a 12 h overnight fast, subjects were given a standard mixed test meal which contained 56 g fat. Blood samples were taken before the meal and for 5 h after the meal when the subjects were randomly allocated to receive either water (control) or were modified sham fed a low-fat (6 g fat) or moderate-fat (38 g fat) meal. Blood samples were collected for a further 3 h. Compared with the control, modified sham feeding a low- or moderate-fat meal did not provoke an early entry of TAG, analysed in either plasma or TAG-rich lipoprotein (TRL) fraction (density ,1´006 kg/l). The TRL-retinyl ester data showed similar findings. A cephalic phase secretion of pancreatic polypeptide, without a significant increase in cholecystokinin levels, was observed on modified sham feeding. Although these data indicate that modified sham feeding was carried out successfully, analysis of the fat content of the expectorant showed that our subjects may have accidentally ingested a small amount of fat (0´7 g for the low-fat meal and 2´4 g for the moderate-fat meal). Nevertheless, an early TAG peak following modified sham feeding was not demonstrated in the present study, suggesting that significant ingestion of food, and not just orosensory stimulation, is necessary to provoke the release of any TAG stored from a previous meal.

Impacts of pollution and climate change on ombrotrophic Sphagnum species in the UK: analysis of uncertainties in two empirical niche models

A significant challenge in the prediction of climate change impacts on ecosystems and biodiversity is quantifying the sources of uncertainty that emerge within and between different models. Statistical species niche models have grown in popularity, yet no single best technique has been identified reflecting differing performance in different situations. Our aim was to quantify uncertainties associated with the application of 2 complimentary modelling techniques. Generalised linear mixed models (GLMM) and generalised additive mixed models (GAMM) were used to model the realised niche of ombrotrophic Sphagnum species in British peatlands. These models were then used to predict changes in Sphagnum cover between 2020 and 2050 based on projections of climate change and atmospheric deposition of nitrogen and sulphur. Over 90% of the variation in the GLMM predictions was due to niche model parameter uncertainty, dropping to 14% for the GAMM. After having covaried out other factors, average variation in predicted values of Sphagnum cover across UK peatlands was the next largest source of variation (8% for the GLMM and 86% for the GAMM). The better performance of the GAMM needs to be weighed against its tendency to overfit the training data. While our niche models are only a first approximation, we used them to undertake a preliminary evaluation of the relative importance of climate change and nitrogen and sulphur deposition and the geographic locations of the largest expected changes in Sphagnum cover. Predicted changes in cover were all small (generally <1% in an average 4 m2 unit area) but also highly uncertain. Peatlands expected to be most affected by climate change in combination with atmospheric pollution were Dartmoor, Brecon Beacons and the western Lake District.

Global climate change and soil carbon stocks: predictions from two contrasting models for the turnover of organic carbon in soil

Enhanced release of CO2 to the atmosphere from soil organic carbon as a result of increased temperatures may lead to a positive feedback between climate change and the carbon cycle, resulting in much higher CO2 levels and accelerated lobal warming. However, the magnitude of this effect is uncertain and critically dependent on how the decomposition of soil organic C (heterotrophic respiration) responds to changes in climate. Previous studies with the Hadley Centre’s coupled climate–carbon cycle general circulation model (GCM) (HadCM3LC) used a simple, single-pool soil carbon model to simulate the response. Here we present results from numerical simulations that use the more sophisticated ‘RothC’ multipool soil carbon model, driven with the same climate data. The results show strong similarities in the behaviour of the two models, although RothC tends to simulate slightly smaller changes in global soil carbon stocks for the same forcing. RothC simulates global soil carbon stocks decreasing by 54 GtC by 2100 in a climate change simulation compared with an 80 GtC decrease in HadCM3LC. The multipool carbon dynamics of RothC cause it to exhibit a slower magnitude of transient response to both increased organic carbon inputs and changes in climate. We conclude that the projection of a positive feedback between climate and carbon cycle is robust, but the magnitude of the feedback is dependent on the structure of the soil carbon model.