A Bayesian design criterion for locating the optimum point on a response surface

Most factorial experiments in industrial research form one stage in a sequence of experiments and so considerable prior knowledge is often available from earlier stages. A Bayesian A-optimality criterion is proposed for choosing designs, when each stage in experimentation consists of a small number of runs and the objective is to optimise a response. Simple formulae for the weights are developed, some examples of the use of the design criterion are given and general recommendations are made. (C) 2003 Elsevier B.V. All rights reserved.

Optimisation of high temperature puffing of potato cubes using response surface methodology

One cubic centimetre potato cubes were blanched, sulfited, dried initially for between 40 and 80 min in air at 90 degreesC in a cabinet drier, puffed in a high temperature fluidised bed and then dried for up to 180 min in a cabinet drier. The final moisture content was 0.05 dwb. The resulting product was optimised using response surface methodology, in terms of volume and colour (L-*, a(*) and b(*) values) of the dry product, as well as rehydration ratio and texture of the rehydrated product. The operating conditions resulting in the optimised product were found to be blanching for 6 min in water at 100 degreesC, dipping in 400 ppm sodium metabisulfite solution for 10 min, initially drying for 40 min and puffing in air at 200 degreesC for 40 s, followed by final drying to a moisture content of 0.05 dwb. (C) 2003 Elsevier Ltd. All rights reserved.

Effect of cutting time, temperature and calcium on curd moisture, whey fat loses and curd yield by response surface methodology

Response surface methodology was used to study the effect of temperature, cutting time, and calcium chloride addition level on curd moisture content, whey fat losses, and curd yield. Coagulation and syneresis were continuously monitored using 2 optical sensors detecting light backscatter. The effect of the factors on the sensors’ response was also examined. Retention of fat during cheese making was found to be a function of cutting time and temperature, whereas curd yield was found to be a function of those 2 factors and the level of calcium chloride addition. The main effect of temperature on curd moisture was to increase the rate at which whey was expelled. Temperature and calcium chloride addition level were also found to affect the light backscatter profile during coagulation whereas the light backscatter profile during syneresis was a function of temperature and cutting time. The results of this study suggest that there is an optimum firmness at which the gel should be cut to achieve maximum retention of fat and an optimum curd moisture content to maximize product yield and quality. It was determined that to maximize curd yield and quality, it is necessary to maximize firmness while avoiding rapid coarsening of the gel network and microsyneresis. These results could contribute to the optimization of the cheese-making process.

Blood-clotting response tests for resistance to diphacinone and chlorophacinone in the Norway rat (Rattus norvegicus Berk.)

Resistance baselines were obtained for the first generation anticoagulant rodenticides chlorophacinone and diphacinone using laboratory, caesarian-derived Norway rats (Rattus norvegicus) as the susceptible strain and the blood clotting response test method. The ED99 estimates for a quantal response were: chlorophacinone, males 0.86 mg kg−1, females 1.03 mg kg−1; diphacinone, males 1.26 mg kg−1, females 1.60 mg kg−1. The dose-response data also showed that chlorophacinone was significantly (p<0.0001) more potent than diphacinone for both male and female rats, and that male rats were more susceptible than females to both compounds (p<0.002). The ED99 doses were then given to groups of five male and five female rats of the Welsh and Hampshire warfarin-resistant strains. Twenty-four hours later, prothrombin times were slightly elevated in both strains but all the animals were classified as resistant to the two compounds, indicating cross-resistance from warfarin to diphacinone and chlorophacinone. When rats of the two resistant strains were fed for six consecutive days on baits containing either diphacinone or chlorophacinone, many animals survived, indicating that their resistance might enable them to survive treatments with these compounds in the field.

Excess open solar magnetic flux from satellite data: 2. A survey of kinematic effects

We investigate the “flux excess” effect, whereby open solar flux estimates from spacecraft increase with increasing heliocentric distance. We analyze the kinematic effect on these open solar flux estimates of large-scale longitudinal structure in the solar wind flow, with particular emphasis on correcting estimates made using data from near-Earth satellites. We show that scatter, but no net bias, is introduced by the kinematic “bunching effect” on sampling and that this is true for both compression and rarefaction regions. The observed flux excesses, as a function of heliocentric distance, are shown to be consistent with open solar flux estimates from solar magnetograms made using the potential field source surface method and are well explained by the kinematic effect of solar wind speed variations on the frozen-in heliospheric field. Applying this kinematic correction to the Omni-2 interplanetary data set shows that the open solar flux at solar minimum fell from an annual mean of 3.82 × 1016 Wb in 1987 to close to half that value (1.98 × 1016 Wb) in 2007, making the fall in the minimum value over the last two solar cycles considerably faster than the rise inferred from geomagnetic activity observations over four solar cycles in the first half of the 20th century.

The importance of temperature and moisture to the egg-laying behaviour of a pest slug, Deroceras reticulatum

Oviposition behaviour is important when modelling the population dynamics of many invertebrates. The numbers of eggs laid are frequently used to describe fecundity, but this measure may differ significantly from realised fecundity. Oviposition has been shown to be important when describing the dynamics of slug populations, which are important agricultural pests. The numbers of eggs laid by Deroceras reticulatum and their viability were measured across a range of 16 temperature (4, 10, 15 and 23 degrees C) by moisture (33%, 42%, 53% and 58% by dry soil weight) experimental combinations. A fitted quadratic response surface model was used to estimate how D. reticulatum adjusted its egg laying to the surrounding temperature and moisture conditions, with most eggs being laid at a combination of 53% soil moisture and 18 degrees C. The number and proportion of viable eggs also covaried with temperature and moisture, suggesting that D. reticulatum may alter their investment in reproduction to maximise their fitness. We have shown that the number of viable eggs differs from the total number of eggs laid by D. reticulatum. Changes in egg viability with temperature and moisture may also be seen in other species and should be considered when modelling populations of egg-laying invertebrates.

Fitting models for the joint action of two drugs using SAS(R)

Combinations of drugs are increasingly being used for a wide variety of diseases and conditions. A pre-clinical study may allow the investigation of the response at a large number of dose combinations. In determining the response to a drug combination, interest may lie in seeking evidence of synergism, in which the joint action is greater than the actions of the individual drugs, or of antagonism, in which it is less. Two well-known response surface models representing no interaction are Loewe additivity and Bliss independence, and Loewe or Bliss synergism or antagonism is defined relative to these. We illustrate an approach to fitting these models for the case in which the marginal single drug dose-response relationships are represented by four-parameter logistic curves with common upper and lower limits, and where the response variable is normally distributed with a common variance about the dose-response curve. When the dose-response curves are not parallel, the relative potency of the two drugs varies according to the magnitude of the desired effect and the models for Loewe additivity and synergism/antagonism cannot be explicitly expressed. We present an iterative approach to fitting these models without the assumption of parallel dose-response curves. A goodness-of-fit test based on residuals is also described. Implementation using the SAS NLIN procedure is illustrated using data from a pre-clinical study. Copyright © 2007 John Wiley & Sons, Ltd.

Modelling climate impact on floods with ensemble climate projections

The evidence provided by modelled assessments of future climate impact on flooding is fundamental to water resources and flood risk decision making. Impact models usually rely on climate projections from global and regional climate models (GCM/RCMs). However, challenges in representing precipitation events at catchment-scale resolution mean that decisions must be made on how to appropriately pre-process the meteorological variables from GCM/RCMs. Here the impacts on projected high flows of differing ensemble approaches and application of Model Output Statistics to RCM precipitation are evaluated while assessing climate change impact on flood hazard in the Upper Severn catchment in the UK. Various ensemble projections are used together with the HBV hydrological model with direct forcing and also compared to a response surface technique. We consider an ensemble of single-model RCM projections from the current UK Climate Projections (UKCP09); multi-model ensemble RCM projections from the European Union's FP6 ‘ENSEMBLES’ project; and a joint probability distribution of precipitation and temperature from a GCM-based perturbed physics ensemble. The ensemble distribution of results show that flood hazard in the Upper Severn is likely to increase compared to present conditions, but the study highlights the differences between the results from different ensemble methods and the strong assumptions made in using Model Output Statistics to produce the estimates of future river discharge. The results underline the challenges in using the current generation of RCMs for local climate impact studies on flooding. Copyright © 2012 Royal Meteorological Society

Twenty-three cycles of changing open solar magnetic flux

This paper presents a comparison of various estimates of the open solar flux, deduced from measurements of the interplanetary magnetic field, from the aa geomagnetic index and from photospheric magnetic field observations. The first two of these estimates are made using the Ulysses discovery that the radial heliospheric field is approximately independent of heliographic latitude, the third makes use of the potential-field source surface method to map the total flux through the photosphere to the open flux at the top of the corona. The uncertainties associated with using the Ulysses result are 5%, but the effects of the assumptions of the potential field source surface method are harder to evaluate. Nevertheless, the three methods give similar results for the last three solar cycles when the data sets overlap. In 11-year running means, all three methods reveal that 1987 marked a significant peak in the long-term variation of the open solar flux. This peak is close to the solar minimum between sunspot cycles 21 and 22, and consequently the mean open flux (averaged from minimum to minimum) is similar for these two cycles. However, this similarity between cycles 21 and 22 in no way implies that the open flux is constant. The long-term variation shows that these cycles are fundamentally different in that the average open flux was rising during cycle 21 (from consistently lower values in cycle 20 and toward the peak in 1987) but was falling during cycle 22 (toward consistently lower values in cycle 23). The estimates from the geomagnetic aa index are unique as they extend from 1842 onwards (using the Helsinki extension). This variation gives strong anticorrelations, with very high statistical significance levels, with cosmic ray fluxes and with the abundances of the cosmogenic isotopes that they produce. Thus observations of photospheric magnetic fields, of cosmic ray fluxes, and of cosmogenic isotope abundances all support the long-term drifts in open solar flux reported by Lockwood et al. [1999a, 1999b].

Optimisation of oat milk formulation to obtain fermented derivatives by using probiotic Lactobacillus reuteri microorganisms

Functional advantages of probiotics combined with interesting composition of oat were considered as an alternative to dairy products. In this study, fermentation of oat milk with Lactobacillus reuteri and Streptococcus thermophilus was analysed to develop a new probiotic product. Central composite design with response surface methodology was used to analyse the effect of different factors (glucose, fructose, inulin and starters) on the probiotic population in the product. Optimised formulation was characterised throughout storage time at 4 ℃ in terms of pH, acidity, β-glucan and oligosaccharides contents, colour and rheological behaviour. All formulations studied were adequate to produce fermented foods and minimum dose of each factor was considered as optimum. The selected formulation allowed starters survival above 107/cfu ml to be considered as a functional food and was maintained during the 28 days controlled. β-glucans remained in the final product with a positive effect on viscosity. Therefore, a new probiotic non-dairy milk was successfully developed in which high probiotic survivals were assured throughout the typical yoghurt-like shelf life.

A functional proteomic method for the enrichment of peripheral membrane proteins reveals the collagen binding protein Hsp47 is exposed on the surface of activated human platelets

Platelets are small blood cells vital for hemostasis. Following vascular damage, platelets adhere to collagens and activate, forming a thrombus that plugs the wound and prevents blood loss. Stimulation of the platelet collagen receptor glycoprotein VI (GPVI) allows recruitment of proteins to receptor-proximal signaling complexes on the inner-leaflet of the plasma membrane. These proteins are often present at low concentrations; therefore, signaling-complex characterization using mass spectrometry is limited due to high sample complexity. We describe a method that facilitates detection of signaling proteins concentrated on membranes. Peripheral membrane proteins (reversibly associated with membranes) were eluted from human platelets with alkaline sodium carbonate. Liquid-phase isoelectric focusing and gel electrophoresis were used to identify proteins that changed in levels on membranes from GPVI-stimulated platelets. Immunoblot analysis verified protein recruitment to platelet membranes and subsequent protein phosphorylation was preserved. Hsp47, a collagen binding protein, was among the proteins identified and found to be exposed on the surface of GPVI-activated platelets. Inhibition of Hsp47 abolished platelet aggregation in response to collagen, while only partially reducing aggregation in response to other platelet agonists. We propose that Hsp47 may therefore play a role in hemostasis and thrombosis.

A time-domain probe method for three-dimensional rough surface reconstructions

The task of this paper is to develop a Time-Domain Probe Method for the reconstruction of impenetrable scatterers. The basic idea of the method is to use pulses in the time domain and the time-dependent response of the scatterer to reconstruct its location and shape. The method is based on the basic causality principle of timedependent scattering. The method is independent of the boundary condition and is applicable for limited aperture scattering data. In particular, we discuss the reconstruction of the shape of a rough surface in three dimensions from time-domain measurements of the scattered field. In practise, measurement data is collected where the incident field is given by a pulse. We formulate the time-domain fieeld reconstruction problem equivalently via frequency-domain integral equations or via a retarded boundary integral equation based on results of Bamberger, Ha-Duong, Lubich. In contrast to pure frequency domain methods here we use a time-domain characterization of the unknown shape for its reconstruction. Our paper will describe the Time-Domain Probe Method and relate it to previous frequency-domain approaches on sampling and probe methods by Colton, Kirsch, Ikehata, Potthast, Luke, Sylvester et al. The approach significantly extends recent work of Chandler-Wilde and Lines (2005) and Luke and Potthast (2006) on the timedomain point source method. We provide a complete convergence analysis for the method for the rough surface scattering case and provide numerical simulations and examples.