Enhancing the selenium content of bovine milk through alteration of the form and concentration of selenium in the diet of the dairy cow

Since estimated dietary selenium intake in the UK has declined steadily from around 60 mug day(-1) in 1975 to 34 mug day(-1) in 1997, there is a need to increase selenium intake from staple foods such as milk and milk products. An experiment was therefore done to investigate the relationship between dietary source and concentration of selenium and the selenium content of bovine milk. In a 3 x 3 factorial design, 90 mid-lactation Holstein dairy cows were supplemented over 8 weeks with either sodium selenite (S), a chelated selenium product (Selenium Metasolate(TM)) (C) or a selenium yeast (Sel-plex(TM)) (Y) at three different dietary inclusion levels of 0.38 (L), 0.76 (M) and 1.14 (H) mg kg(-1) dry matter (DM). Significant increases in milk selenium concentration were observed for all three sources with increasing inclusion level in the diet, but Y gave a much greater response (up to +65 mug l(-1)) than the other two sources of selenium (S and C up to +4 and +6 mug l(-1) respectively). The Y source also resulted in a substantially higher apparent efficiency of transfer of selenium from diet to milk than S or C. Feeding Y at the lowest dietary concentration, and thus within the maximum level permitted under EU regulations, resulted in milk with a selenium concentration of 28 mug l(-1). If the selenium concentration of milk in the UK was increased to this value, it would, at current consumption rates, provide an extra 8.7 mug selenium day(-1), or 11 and 14% of daily recommended national intake for men and women respectively. (C) 2004 Society of Chemical Industry.

Pea-barley intercropping and short-term subsequent crop effects across European organic cropping conditions

Grain legumes are known to increase the soil mineral nitrogen (N) content, reduce the infection pressure of soil borne pathogens, and hence enhance subsequent cereals yields. Replicated field experiments were performed throughout W. Europe (Denmark, United Kingdom, France, Germany and Italy) to asses the effect of intercropping pea and barley on the N supply to subsequent wheat in organic cropping systems. Pea and barley were grown either as sole crops at the recommended plant density (P100 and B100, respectively) or in replacement (P50B50) or additive (P100B50) intercropping designs. In the replacement design the total relative plant density is kept constant, while the additive design uses the optimal sole crop density for pea supplementing with 'extra' barley plants. The pea and barley crops were followed by winter wheat with and without N application. Additional experiments in Denmark and the United Kingdom included subsequent spring wheat with grass-clover as catch crops. The experiment was repeated over the three cropping seasons of 2003, 2004 and 2005. Irrespective of sites and intercrop design pea-barley intercropping improved the plant resource utilization (water, light, nutrients) to grain N yield with 25-30% using the Land Equivalent ratio. In terms of absolute quantities, sole cropped pea accumulated more N in the grains as compared to the additive design followed by the replacement design and then sole cropped barley. The post harvest soil mineral N content was unaffected by the preceding crops. Under the following winter wheat, the lowest mineral N content was generally found in early spring. Variation in soil mineral N content under the winter wheat between sites and seasons indicated a greater influence of regional climatic conditions and long-term cropping history than annual preceding crop and residue quality. Just as with the soil mineral N, the subsequent crop response to preceding crop was negligible. Soil N balances showed general negative values in the 2-year period, indicating depletion of N independent of preceding crop and cropping strategy. It is recommended to develop more rotational approaches to determine subsequent crop effects in organic cropping systems, since preceding crop effects, especially when including legumes, can occur over several years of cropping.

A field experiment to recreate species rich hay meadows using regional seed mixtures

Species rich semi-natural grasslands are an important but threatened habitat throughout Europe and much of the former area has been lost since the 1950s. However, in some countries large areas have been preserved and the demand for meadow recreation by sowing seed mixtures is increasing. In the White Carpathians Protected Landscape Area (Czech Republic) the use of commercial seed mixtures is undesirable and the use of regional mixtures has been investigated. The costs for seeding large areas are high and lower cost techniques are needed. In 1999 a field experiment was set up to investigate the establishment of hay meadow vegetation comparing sowing a regional mixture all over a plot with sowing narrow 2.5 In strips of regional seed mixtures into a matrix of a commercial grass mixture or into natural regeneration. The results after five seasons showed good establishment of the sown species in the meadow treatment. Spread of sown species from the sown strips into the surrounding matrix occurred but the cover of species was lower in the commercial grass matrix compared with the natural regeneration matrix. Colonisation of some plots by unsown desirable grassland species from adjacent grassland habitats also occurred, but more species colonised the natural regeneration matrix than the commercial grasses or the sown meadow matrix itself. Overall, the results indicate that, in appropriate situations, sown strips can provide a lower cost but slower and longer-term alternative to field scale sowing of regional seed mixtures for recreation of hay meadow vegetation. (C) 2007 Elsevier Ltd. All rights reserved.

Heterosis for yield and its physiological determinants in wheat

Heterosis in hybrid wheat varieties produced using a chemical hybridising agent was assessed in field experiments. Hyno Esta and its parents were compared in factorial combinations of four-seed rates (25-300 seeds m(-2)) and two nitrogen fertilizer rates (0 and 200 kg N ha(-1)) in 2001/02 and again in 2002/03. Hyno Rista and Hyno Renta and their parents were compared at two-seed rates in 2001/02. Hyno Rista and its parents were added factorially to the Hyno Esta experiment in 2002/03, while Hyno Renta and Hybred and their parents were compared at two seed rates in 2002/03. Mid parent heterosis for grain yield was found in three hybrids and two of these showed high parent heterosis. High parent heterosis in Hyno Esta over a range of sowing densities was mostly exhibited in total biomass but also, in one of two years, in harvest index. High parent heterosis in Hyno Renta was associated more with harvest index than with biomass. The heterosis for biomass in Hyno Esta resulted from greater interception of photosynthetically active radiation (PAR) than the male parent, with better radiation use efficiency than the female parent. In both seasons Hyno Esta achieved grain numbers per ear at least as high as the high parent for this trait (Audace), and combined this with mean grain weights at least as heavy as the high parent for mean grain weight (Estica). Much of the increased biomass and grain yield in the hybrid came late in the season as high parent heterosis was expressed for both maximum grain filling rate and grain filling duration. Heterosis was higher when nitrogen was applied than when withheld; only greater at lower seed rates when expressed in proportionate terms (e.g. as a percentage of the parents), rather than in absolute terms (e.g. t ha(-1)); and greater in the year with the cooler and wetter summer.

Choice of designs and doses for early phase trials

This paper reviews state-of-art statistical designs for dose-escalation procedures in first-into-man studies. The main focus will be on studies in oncology, as most statistical procedures for phase I trials have been proposed in this context. Extensions to situations such as the observation of bivariate outcomes and healthy volunteer studies are also discussed. The number of dose levels and cohort sizes used in early phase trials are considered. Finally, this paper raises some practical issues for dose-escalation procedures.

An evaluation of Bayesian designs for dose-escalation studies in healthy volunteers

In this paper, Bayesian decision procedures previously proposed for dose-escalation studies in healthy volunteers are reviewed and evaluated. Modifications are made to the expression of the prior distribution in order to make the procedure simpler to implement and a more relevant criterion for optimality is introduced. The results of an extensive simulation exercise to establish the proper-ties of the procedure and to aid choice between designs are summarized, and the way in which readers can use simulation to choose a design for their own trials is described. The influence of the value of the within-subject correlation on the procedure is investigated and the use of a simple prior to reflect uncertainty about the correlation is explored. Copyright (c) 2005 John Wiley & Sons, Ltd.

An adaptive approach to implementing bivariate group sequential clinical trial designs

In clinical trials, situations often arise where more than one response from each patient is of interest; and it is required that any decision to stop the study be based upon some or all of these measures simultaneously. Theory for the design of sequential experiments with simultaneous bivariate responses is described by Jennison and Turnbull (Jennison, C., Turnbull, B. W. (1993). Group sequential tests for bivariate response: interim analyses of clinical trials with both efficacy and safety endpoints. Biometrics 49:741-752) and Cook and Farewell (Cook, R. J., Farewell, V. T. (1994). Guidelines for monitoring efficacy and toxicity responses in clinical trials. Biometrics 50:1146-1152) in the context of one efficacy and one safety response. These expositions are in terms of normally distributed data with known covariance. The methods proposed require specification of the correlation, ρ between test statistics monitored as part of the sequential test. It can be difficult to quantify ρ and previous authors have suggested simply taking the lowest plausible value, as this will guarantee power. This paper begins with an illustration of the effect that inappropriate specification of ρ can have on the preservation of trial error rates. It is shown that both the type I error and the power can be adversely affected. As a possible solution to this problem, formulas are provided for the calculation of correlation from data collected as part of the trial. An adaptive approach is proposed and evaluated that makes use of these formulas and an example is provided to illustrate the method. Attention is restricted to the bivariate case for ease of computation, although the formulas derived are applicable in the general multivariate case.

Sequential designs for phase III clinical trials incorporating treatment selection

Most statistical methodology for phase III clinical trials focuses on the comparison of a single experimental treatment with a control. An increasing desire to reduce the time before regulatory approval of a new drug is sought has led to development of two-stage or sequential designs for trials that combine the definitive analysis associated with phase III with the treatment selection element of a phase II study. In this paper we consider a trial in which the most promising of a number of experimental treatments is selected at the first interim analysis. This considerably reduces the computational load associated with the construction of stopping boundaries compared to the approach proposed by Follman, Proschan and Geller (Biometrics 1994; 50: 325-336). The computational requirement does not exceed that for the sequential comparison of a single experimental treatment with a control. Existing methods are extended in two ways. First, the use of the efficient score as a test statistic makes the analysis of binary, normal or failure-time data, as well as adjustment for covariates or stratification straightforward. Second, the question of trial power is also considered, enabling the determination of sample size required to give specified power. Copyright © 2003 John Wiley & Sons, Ltd.

Decision-theoretic designs for phase II clinical trials allowing for competing studies

This article describes an approach to optimal design of phase II clinical trials using Bayesian decision theory. The method proposed extends that suggested by Stallard (1998, Biometrics54, 279–294) in which designs were obtained to maximize a gain function including the cost of drug development and the benefit from a successful therapy. Here, the approach is extended by the consideration of other potential therapies, the development of which is competing for the same limited resources. The resulting optimal designs are shown to have frequentist properties much more similar to those traditionally used in phase II trials.

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.

A practical comparison of group-sequential and adaptive designs

Sequential methods provide a formal framework by which clinical trial data can be monitored as they accumulate. The results from interim analyses can be used either to modify the design of the remainder of the trial or to stop the trial as soon as sufficient evidence of either the presence or absence of a treatment effect is available. The circumstances under which the trial will be stopped with a claim of superiority for the experimental treatment, must, however, be determined in advance so as to control the overall type I error rate. One approach to calculating the stopping rule is the group-sequential method. A relatively recent alternative to group-sequential approaches is the adaptive design method. This latter approach provides considerable flexibility in changes to the design of a clinical trial at an interim point. However, a criticism is that the method by which evidence from different parts of the trial is combined means that a final comparison of treatments is not based on a sufficient statistic for the treatment difference, suggesting that the method may lack power. The aim of this paper is to compare two adaptive design approaches with the group-sequential approach. We first compare the form of the stopping boundaries obtained using the different methods. We then focus on a comparison of the power of the different trials when they are designed so as to be as similar as possible. We conclude that all methods acceptably control type I error rate and power when the sample size is modified based on a variance estimate, provided no interim analysis is so small that the asymptotic properties of the test statistic no longer hold. In the latter case, the group-sequential approach is to be preferred. Provided that asymptotic assumptions hold, the adaptive design approaches control the type I error rate even if the sample size is adjusted on the basis of an estimate of the treatment effect, showing that the adaptive designs allow more modifications than the group-sequential method.

Variance estimation with highly stratified sampling designs with unequal probabilities

It is common practice to design a survey with a large number of strata. However, in this case the usual techniques for variance estimation can be inaccurate. This paper proposes a variance estimator for estimators of totals. The method proposed can be implemented with standard statistical packages without any specific programming, as it involves simple techniques of estimation, such as regression fitting.

Epidemiological risk assessment using linked network and grid based modelling: Phytophthora ramorum and Phytophthora kernoviae in the UK

We developed a stochastic simulation model incorporating most processes likely to be important in the spread of Phytophthora ramorum and similar diseases across the British landscape (covering Rhododendron ponticum in woodland and nurseries, and Vaccinium myrtillus in heathland). The simulation allows for movements of diseased plants within a realistically modelled trade network and long-distance natural dispersal. A series of simulation experiments were run with the model, representing an experiment varying the epidemic pressure and linkage between natural vegetation and horticultural trade, with or without disease spread in commercial trade, and with or without inspections-with-eradication, to give a 2 x 2 x 2 x 2 factorial started at 10 arbitrary locations spread across England. Fifty replicate simulations were made at each set of parameter values. Individual epidemics varied dramatically in size due to stochastic effects throughout the model. Across a range of epidemic pressures, the size of the epidemic was 5-13 times larger when commercial movement of plants was included. A key unknown factor in the system is the area of susceptible habitat outside the nursery system. Inspections, with a probability of detection and efficiency of infected-plant removal of 80% and made at 90-day intervals, reduced the size of epidemics by about 60% across the three sectors with a density of 1% susceptible plants in broadleaf woodland and heathland. Reducing this density to 0.1% largely isolated the trade network, so that inspections reduced the final epidemic size by over 90%, and most epidemics ended without escape into nature. Even in this case, however, major wild epidemics developed in a few percent of cases. Provided the number of new introductions remains low, the current inspection policy will control most epidemics. However, as the rate of introduction increases, it can overwhelm any reasonable inspection regime, largely due to spread prior to detection. (C) 2009 Elsevier B.V. All rights reserved.