Survival analysis statistics applied to threshold data obtained from the ascending forced-choice method of limits

Sensory thresholds are often collected through ascending forced-choice methods. Group thresholds are important for comparing stimuli or populations; yet, the method has two problems. An individual may correctly guess the correct answer at any concentration step and might detect correctly at low concentrations but become adapted or fatigued at higher concentrations. The survival analysis method deals with both issues. Individual sequences of incorrect and correct answers are adjusted, taking into account the group performance at each concentration. The technique reduces the chance probability where there are consecutive correct answers. Adjusted sequences are submitted to survival analysis to determine group thresholds. The technique was applied to an aroma threshold and a taste threshold study. It resulted in group thresholds similar to ASTM or logarithmic regression procedures. Significant differences in taste thresholds between younger and older adults were determined. The approach provides a more robust technique over previous estimation methods.

Coupling liquid MALDI MS to liquid chromatography

Matrix-assisted laser desorption/ionisation (MALDI) coupled with time-of-flight (TOF) mass spectrometry (MS) is a powerful tool for the analysis of biological samples, and nanoflow high-performance liquid chromatography (nanoHPLC) is a useful separation technique for the analysis of complex proteomics samples. The off-line combination of MALDI and nanoHPLC has been extensively investigated and straightforward techniques have been developed, focussing particularly on automated MALDI sample preparation that yields sensitive and reproducible spectra. Normally conventional solid MALDI matrices such as α-cyano-4-hydroxycinnamic acid (CHCA) are used for sample preparation. However, they have limited usefulness in quantitative measurements and automated data acquisition because of the formation of heterogeneous crystals, resulting in highly variable ion yields and desorption/ ionization characteristics. Glycerol-based liquid support matrices (LSM) have been proposed as an alternative to the traditional solid matrices as they provide increased shot-to-shot reproducibility, leading to prolonged and stable ion signals and therefore better results. This chapter focuses on the integration of the liquid LSM MALDI matrices into the LC-MALDI MS/MS approach in identifying complex and large proteomes. The interface between LC and MALDI consists of a robotic spotter, which fractionates the eluent from the LC column into nanoliter volumes, and co-spots simultaneously the liquid matrix with the eluent fractions onto a MALDI target plate via sheath flow. The efficiency of this method is demonstrated through the analysis of trypsin digests of both bovine serum albumin (BSA) and Lactobacillus plantarum WCFS1 proteins.

Reduction of arsenic uptake by lettuce with ferrous sulfate applied to contaminated soil

Soil contamination by arsenic (As) presents a hazard in many countries and there is a need for techniques to minimize As uptake by plants. A proposed in situ remediation method was tested by growing lettuce (Lactuca sativa L. cv. Kermit) in a greenhouse pot experiment on soil that contained 577 mg As kg(-1), taken from a former As smelter site. All combinations of iron (Fe) oxides, at concentrations of 0.00, 0.22, 0.54, and 1.09% (w/w), and lime, at concentrations of 0.00, 0.27, 0.68, and 1.36% (w/w), were tested in a factorial design. To create the treatments, field-moist soil, commercial-grade FeSO4, and ground agricultural lime were mixed and stored for one week, allowing Fe oxides to precipitate. Iron oxides gave highly significant (P < 0.001) reductions in lettuce As concentrations, down to 11% of the lettuce As concentration for untreated soil. For the Fe oxides and lime treatment combinations where soil pH was maintained nearly constant, the lettuce As concentration declined in an exponential relationship with increasing FeSO4 application rate and lettuce yield was almost unchanged. Iron oxides applied at a concentration of 1.09% did not give significantly lower lettuce As concentrations than the 0.54% treatment. Simultaneous addition of lime with FeSO4 was essential. Ferrous sulfate with insufficient lime lowered soil pH and caused mobilization of Al, Ba, Co, Cr, Cu, Fe, K, Mg, Mn, Na, Ni, Pb, Sr, and Zn. At the highest Fe oxide to lime ratios, Mn toxicity caused severe yield loss.

Reduction of arsenic uptake by lettuce with ferrous sulfate applied to contaminated soil

Soil contamination by arsenic (As) presents a hazard in many countries and there is a need for techniques to minimize As uptake by plants. A proposed in situ remediation method was tested by growing lettuce (Lactuca sativa L. cv. Kermit) in a greenhouse pot experiment on soil that contained 577 mg As kg(-1), taken from a former As smelter site. All combinations of iron (Fe) oxides, at concentrations of 0.00, 0.22, 0.54, and 1.09% (w/w), and lime, at concentrations of 0.00, 0.27, 0.68, and 1.36% (w/w), were tested in a factorial design. To create the treatments, field-moist soil, commercial-grade FeSO4, and ground agricultural lime were mixed and stored for one week, allowing Fe oxides to precipitate. Iron oxides gave highly significant (P < 0.001) reductions in lettuce As concentrations, down to 11% of the lettuce As concentration for untreated soil. For the Fe oxides and lime treatment combinations where soil pH was maintained nearly constant, the lettuce As concentration declined in an exponential relationship with increasing FeSO4 application rate and lettuce yield was almost unchanged. Iron oxides applied at a concentration of 1.09% did not give significantly lower lettuce As concentrations than the 0.54% treatment. Simultaneous addition of lime with FeSO4 was essential. Ferrous sulfate with insufficient lime lowered soil pH and caused mobilization of Al, Ba, Co, Cr, Cu, Fe, K, Mg, Mn, Na, Ni, Pb, Sr, and Zn. At the highest Fe oxide to lime ratios, Mn toxicity caused severe yield loss.

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.

Clustering techniques applied to multiple-models structures

In this work a method for building multiple-model structures is presented. A clustering algorithm that uses data from the system is employed to define the architecture of the multiple-model, including the size of the region covered by each model, and the number of models. A heating ventilation and air conditioning system is used as a testbed of the proposed method.

New fast SCFT algorithm applied to binary diblock copolymer/homopolymer blends

We present an efficient strategy for mapping out the classical phase behavior of block copolymer systems using self-consistent field theory (SCFT). With our new algorithm, the complete solution of a classical block copolymer phase can be evaluated typically in a fraction of a second on a single-processor computer, even for highly segregated melts. This is accomplished by implementing the standard unit-cell approximation (UCA) for the cylindrical and spherical phases, and solving the resulting equations using a Bessel function expansion. Here the method is used to investigate blends of AB diblock copolymer and A homopolymer, concentrating on the situation where the two molecules are of similar size.

Application of cylindrical distribution functions to wide-angle X-ray scattering from oriented polymers

A systematic approach is presented for obtaining cylindrical distribution functions (CDF's) of noncrystalline polymers which have been oriented by extension. The scattering patterns and CDF's are also sharpened by the method proposed by Deas and by Ruland. Data from atactic poly(methyl methacrylate) and polystyrene are analysed by these techniques. The methods could also be usefully applied to liquid crystals.

Restriction fragment length polymorphism of polymerase chain reaction products applied to the differentiation of poultry campylobacters for epidemiological investigations

A technique for subtyping Camplobacter jejuni isolates has been developed by using the restriction fragment length polymorphism (Rnp) of polymerase chain reaction (PCR) products of the fluA and flaB genes. The technique was validated by using strains representing 28 serotypes of C jejuni and it may also be applied to C coli. From these strains 12 distinct RFLP profiles were observed but there was no direct relationship between the RFLP profile and the serotype. One hundred and thirty-five campylobacter isolates from 15 geographically distinct broiler flocks were investigated. All the isolates could be subtyped by using the RFLP method. Isolates from most of the flocks had a single RFLP profile despite data indicating that several serotypes were involved. Although it is possible that further restriction analysis may have demonstrated profile variations in these strains, it is more likely that antigenic variation can occur within genotypically related campylobacters. As a result, serotyping may give conflicting information for veterinary epidemiological purposes. This RFLP typing scheme appears to provide a suitable tool for the investigation of the sources and routes of transmission of campylobacters in chickens.

Detection of Chlamydia psittaci DNA in avian clinical samples by polymerase chain reaction

A polymerase chain reaction (PCR) assay was developed to detect Chlamydia psittaci DNA in faeces and tissue samples from avian species. Primers were designed to amplify a 264 bp product derived from part of the 5' non-translated region and part of the coding region of the ompA gene which encodes the major outer membrane protein. Amplified sequences were confirmed by Southern hybridization using an internal probe. The sensitivity of the combined assay was found to be between 60 to 600 fg of chlamydial DNA (approximately 6 to 60 genome copies). The specificity of the assay was confirmed since PCR product was not obtained from samples containing several serotypes of C. trachomatis, strains of C. pneumoniae, the type strain of C. pecorum, nor from samples containing microorganisms commonly found in the avian gut flora. In this study, 404 avian faeces and 141 avian tissue samples received by the Central Veterinary Laboratory over a 6 month period were analysed by PCR, antigen detection ELISA and where possible, cell culture isolation. PCR performed favourably compared with ELISA and cell culture, or with ELISA alone. The PCR assay was especially suited to the detection of C. psittaci DNA in avian faeces samples. The test was also useful when applied to tissue samples from small contact birds associated with a case of human psittacosis where ELISA results were negative and chlamydial isolation was a less favourable method due to the need for rapid diagnosis.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

Species distribution models for crop pollination: a modelling framework applied to Great Britain

Insect pollination benefits over three quarters of the world's major crops. There is growing concern that observed declines in pollinators may impact on production and revenues from animal pollinated crops. Knowing the distribution of pollinators is therefore crucial for estimating their availability to pollinate crops; however, in general, we have an incomplete knowledge of where these pollinators occur. We propose a method to predict geographical patterns of pollination service to crops, novel in two elements: the use of pollinator records rather than expert knowledge to predict pollinator occurrence, and the inclusion of the managed pollinator supply. We integrated a maximum entropy species distribution model (SDM) with an existing pollination service model (PSM) to derive the availability of pollinators for crop pollination. We used nation-wide records of wild and managed pollinators (honey bees) as well as agricultural data from Great Britain. We first calibrated the SDM on a representative sample of bee and hoverfly crop pollinator species, evaluating the effects of different settings on model performance and on its capacity to identify the most important predictors. The importance of the different predictors was better resolved by SDM derived from simpler functions, with consistent results for bees and hoverflies. We then used the species distributions from the calibrated model to predict pollination service of wild and managed pollinators, using field beans as a test case. The PSM allowed us to spatially characterize the contribution of wild and managed pollinators and also identify areas potentially vulnerable to low pollination service provision, which can help direct local scale interventions. This approach can be extended to investigate geographical mismatches between crop pollination demand and the availability of pollinators, resulting from environmental change or policy scenarios.

A single-column model ensemble approach applied to the TWP-ICE experiment

Single-column models (SCM) are useful test beds for investigating the parameterization schemes of numerical weather prediction and climate models. The usefulness of SCM simulations are limited, however, by the accuracy of the best estimate large-scale observations prescribed. Errors estimating the observations will result in uncertainty in modeled simulations. One method to address the modeled uncertainty is to simulate an ensemble where the ensemble members span observational uncertainty. This study first derives an ensemble of large-scale data for the Tropical Warm Pool International Cloud Experiment (TWP-ICE) based on an estimate of a possible source of error in the best estimate product. These data are then used to carry out simulations with 11 SCM and two cloud-resolving models (CRM). Best estimate simulations are also performed. All models show that moisture-related variables are close to observations and there are limited differences between the best estimate and ensemble mean values. The models, however, show different sensitivities to changes in the forcing particularly when weakly forced. The ensemble simulations highlight important differences in the surface evaporation term of the moisture budget between the SCM and CRM. Differences are also apparent between the models in the ensemble mean vertical structure of cloud variables, while for each model, cloud properties are relatively insensitive to forcing. The ensemble is further used to investigate cloud variables and precipitation and identifies differences between CRM and SCM particularly for relationships involving ice. This study highlights the additional analysis that can be performed using ensemble simulations and hence enables a more complete model investigation compared to using the more traditional single best estimate simulation only.

Spherical nonparametric estimators applied to the UGAMP model integration for AMIP

The aim of this paper is essentially twofold: first, to describe the use of spherical nonparametric estimators for determining statistical diagnostic fields from ensembles of feature tracks on a global domain, and second, to report the application of these techniques to data derived from a modern general circulation model. New spherical kernel functions are introduced that are more efficiently computed than the traditional exponential kernels. The data-driven techniques of cross-validation to determine the amount elf smoothing objectively, and adaptive smoothing to vary the smoothing locally, are also considered. Also introduced are techniques for combining seasonal statistical distributions to produce longer-term statistical distributions. Although all calculations are performed globally, only the results for the Northern Hemisphere winter (December, January, February) and Southern Hemisphere winter (June, July, August) cyclonic activity are presented, discussed, and compared with previous studies. Overall, results for the two hemispheric winters are in good agreement with previous studies, both for model-based studies and observational studies.

Influence of time on the plant availability of Cd, Ni and Zn after sewage sludge has been applied to soils

A pot experiment was conducted to test the hypothesis that decomposition of organic matter in sewage sludge and the consequent formation of dissolved organic compounds (DOC) would lead to an increase in the bioavailability of the heavy metals. Two Brown Earth soils, one with clayey loam texture (CL) and the other a loamy sand (LS) were mixed with sewage sludge at rates equivalent to 0, 10 and 50 1 dry sludge ha(-1) and the pots were sown with ryegrass (Lolium perenne L.). The organic matter content and heavy metal availability assessed with soil extractions with 0.05 M CaCl2 were monitored over a residual time of two years, while plant uptake over one year, after addition of the sludge. It was found that the concentrations of Cd and Ni in both the ryegrass and the soil extracts increased slightly but significantly during the first year. In most cases, this increase was most evident especially at the higher sludge application rate (50 t ha(-1)). However, in the second year metal availability reached a plateau. Zinc concentrations in the ryegrass did not show an increase but the CaCl2 extracts increased during the first year. In contrast, organic matter content decreased rapidly in the first months of the first year and much more slowly in the second (total decrease of 16%). The concentrations of DOC increased significantly in the more organic rich CL soil in the course of two years. The pattern followed by the decomposition of organic matter with time and the production of DOC may provide at least a partial explanation for trend towards increased metal availability.