High-breakdown linear discriminant analysis

The classification rules of linear discriminant analysis are defined by the true mean vectors and the common covariance matrix of the populations from which the data come. Because these true parameters are generally unknown, they are commonly estimated by the sample mean vector and covariance matrix of the data in a training sample randomly drawn from each population. However, these sample statistics are notoriously susceptible to contamination by outliers, a problem compounded by the fact that the outliers may be invisible to conventional diagnostics. High-breakdown estimation is a procedure designed to remove this cause for concern by producing estimates that are immune to serious distortion by a minority of outliers, regardless of their severity. In this article we motivate and develop a high-breakdown criterion for linear discriminant analysis and give an algorithm for its implementation. The procedure is intended to supplement rather than replace the usual sample-moment methodology of discriminant analysis either by providing indications that the dataset is not seriously affected by outliers (supporting the usual analysis) or by identifying apparently aberrant points and giving resistant estimators that are not affected by them.

Prediction of Golgi Type II membrane proteins based on their transmembrane domains

Motivation: A major issue in cell biology today is how distinct intracellular regions of the cell, like the Golgi Apparatus, maintain their unique composition of proteins and lipids. The cell differentially separates Golgi resident proteins from proteins that move through the organelle to other subcellular destinations. We set out to determine if we could distinguish these two types of transmembrane proteins using computational approaches. Results: A new method has been developed to predict Golgi membrane proteins based on their transmembrane domains. To establish the prediction procedure, we took the hydrophobicity values and frequencies of different residues within the transmembrane domains into consideration. A simple linear discriminant function was developed with a small number of parameters derived from a dataset of Type II transmembrane proteins of known localization. This can discriminate between proteins destined for Golgi apparatus or other locations (post-Golgi) with a success rate of 89.3% or 85.2%, respectively on our redundancy-reduced data sets.

Computational differentiation of N-terminal signal peptides and transmembrane helices

Signal peptides and transmembrane helices both contain a stretch of hydrophobic amino acids. This common feature makes it difficult for signal peptide and transmembrane helix predictors to correctly assign identity to stretches of hydrophobic residues near the N-terminal methionine of a protein sequence. The inability to reliably distinguish between N-terminal transmembrane helix and signal peptide is an error with serious consequences for the prediction of protein secretory status or transmembrane topology. In this study, we report a new method for differentiating protein N-terminal signal peptides and transmembrane helices. Based on the sequence features extracted from hydrophobic regions (amino acid frequency, hydrophobicity, and the start position), we set up discriminant functions and examined them on non-redundant datasets with jackknife tests. This method can incorporate other signal peptide prediction methods and achieve higher prediction accuracy. For Gram-negative bacterial proteins, 95.7% of N-terminal signal peptides and transmembrane helices can be correctly predicted (coefficient 0.90). Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 99% (coefficient 0.92). For eukaryotic proteins, 94.2% of N-terminal signal peptides and transmembrane helices can be correctly predicted with coefficient 0.83. Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 87% (coefficient 0.85). The method can be used to complement current transmembrane protein prediction and signal peptide prediction methods to improve their prediction accuracies. (C) 2003 Elsevier Inc. All rights reserved.

Detecting sugarcane 'orange rust' disease using EO-1 Hyperion hyperspectral imagery

This Letter evaluates several narrow-band indices from EO-1 Hyperion imagery in discriminating sugarcane areas affected by 'orange rust' ( Puccinia kuehnii ) disease. Forty spectral vegetation indices (SVIs), focusing on bands related to leaf pigments, leaf internal structure, and leaf water content, were generated from an image acquired over Mackay, Queensland, Australia. Discriminant function analysis was used to select an optimum set of indices based on their correlations with the discriminant function. The predictive ability of each index was also assessed based on the accuracy of classification. Results demonstrated that Hyperion imagery can be used to detect orange rust disease in sugarcane crops. While some indices that only used visible near-infrared (VNIR) bands (e.g. SIPI and R800/R680) offer separability, the combination of VNIR bands with the moisture-sensitive band (1660 nm) yielded increased separability of rust-affected areas. The newly formulated 'Disease-Water Stress Indices' (DWSI-1=R800/R1660; DSWI-2=R1660/R550; DWSI-5=(R800+R550)/(R1660+R680)) produced the largest correlations, indicating their superior ability to discriminate sugarcane areas affected by orange rust disease.

A phylogeny for the genus Isoodon and a range extension for I-obesulus peninsulae based on mtDNA control region and morphology

Short-nosed bandicoots, Isoodon, have undergone marked range contractions since European colonisation of Australia and are currently divided into many subspecies, the validity of which is debated. Discriminant function analysis of morphology and a phylogeny of Isoodon based on mtDNA control region sequences indicate a clear split between two of the three recognised species, I. macrourus and I. obesulus/auratus. However, while all previously recognised taxa within the I. obesulus/auratus group are morphologically distinct, I. auratus and I. obesulus are not phylogenetically distinct for mtDNA. The genetic divergence between I. obesulus and I. auratus (2.6%) is similar to that found among geographic isolates of the former (I. o. obesulus and I. o. peninsulae: 2.7%). Further, the divergence between geographically close populations of two different species (I. o. obesulus from Western Australia and I. a. barrowensis: 1.2%) is smaller than that among subspecies within I. auratus (I. a. barrowensis and I. auratus from northern Western Australia: 1.7%). A newly discovered population of Isoodon in the Lamb Range, far north Queensland, sympatric with a population of I. m. torosus, is shown to represent a range extension of I. o. peninsulae (350 km). It seems plausible that what is currently considered as two species, I. obesulus and I. auratus, was once one continuous species now represented by isolated populations that have diverged morphologically as a consequence of adaptation to the diverse environments that occur throughout their range. The taxonomy of these populations is discussed in relation to their morphological distinctiveness and genetic similarity.

Stock discrimination and movements of narrow-barred Spanish mackerel across northern Australia as indicated by parasites

The parasite fauna of Spanish mackerel Scomberomorus commerson from 10 sites across northern Australia and one site in Indonesia, was examined to evaluate the degree of movement and subsequent stock structure of the fish. Kupang fish (Indonesia) had very few Terranova spp.. Grillotia branchi, Otobothrium cysticum or Pterobothrium sp. compared to Australian fish, indicating that no Australian fish enter the Kupang fishery. Univariate and discriminant function analysis of four 'temporary' parasite species, the copepod Pseudocyenoides armatus and the monogeneans Gotocotyla bivaginalis, Pricea multae and Pseudothoracocotyla ovalis, demonstrated little similarity between areas of northern Australia, indicating minimal short-term exchange between neighbouring groups of S. commerson. Analyses of five 'permanent' parasite species, the larval helminths G. branchi, O. cysticum, Pterobothrium sp., Callitetrarhynchus gracilis and Paranybelinia balli, also revealed large differences between areas thus indicating long-term separation. There are at least six parasitological stocks across northern Australia: Fog Bay/Bathurst Island, Cape Wessel. Groote/Sir Edward Pellew. Mornington Island, Weipa. and the Torres Strait. The occurrence of a few irregular fish in the samples suggested that LIP to 5% of fish moved between stocks during their lifetime. The similarity of within-school variability to that between schools showed that the fish do not form long-term school associations. (C) 2003 The Fisheries Society of the British Isles.

Factors associated with the wish to hasten death: a study of patients with terminal illness

Background. There is considerable debate regarding the clinical issues surrounding the wish to hasten death (WTHD) in the terminally ill. The clinical factors contributing to the WTHD need further investigation among the terminally ill in order to enhance understanding of the clinical assessment and treatment needs that underlie this problem. A more detailed understanding may assist with the development of appropriate therapeutic interventions. Method. A sample of terminally ill cancer patients (N=256) recruited from an in-patient hospice unit, home palliative care service and a general hospital palliative care consulting service from Brisbane Australia between 1998-2001 completed a questionnaire assessing psychological (depression and anxiety), social (family relationship, social support, level of burden on others) and the impact of physical symptoms. The association between these factors and the WTHD was investigated. Results. A high WTHD was reported by 14% of patients. A discriminant function analysis revealed that the following variables were associated with a high WTHD (P

Comparative habitat use in Juniper woodland bird community

We compared vegetation structure used by 14 bird species during the 1998 and 1999 breeding seasons to determine what habitat features best accounted for habitat division and community organization in Utah juniper (Juniperus osteosperma) woodlands of southwestern Wyoming. Habitat use was quantified by measuring 24 habitat variables in 461 bird-centered quadrats, each 0.04 ha in size. Using discriminant function analysis, we differentiated between habitat used by 14 bird species along 3 habitat dimensions: (1) variation in shrub cover, overstory juniper cover, mature tree density, understory height, and decadent tree density; (2) a gradient composed of elevation and forb cover; and (3) variation in grass cover, tree height, seedling/sapling cover, and bare ground/rock cover. Of 14 species considered, 9 exhibited substantial habitat partitioning: Mourning Dove (Zenaida macroura), Bewick's Wren (Thryomanes bewickii), Blue-gray Gnatcatcher (Polioptila caerulea), Mountain Bluebird (Sialia currucoides), Plumbeous Vireo (Vireo plumbeus), Green-tailed Towhee (Pipilo chlorurus), Brewer's Sparrow (Spizella breweri), Dark-eyed Junco (Junco hyemalis), and Cassin's Finch (Carpodacus cassinii). Our results indicate juniper bird communities of southwestern Wyoming are organized along a 3-dimensional habitat gradient composed of woodland maturity, elevation, and juniper recruitment. Because juniper birds partition habitat along successional and altitudinal gradients, indiscriminate woodland clearing as well as continued fire suppression will alter species composition. Restoration efforts should ensure that all successional stages of juniper woodland are present on the landscape.

Use of X-radiographs to distinguish members of the Montastraea annularis reef-coral species complex

Recent work suggests the Montastraea annularis species complex consists of at least three species, which can be distinguished qualitatively in the field using features related to colony growth (e.g. overall growth form. bumpiness, growth along the colony edge). However, when whole colonies are not available and surfaces are eroded, identification becomes problematic when relying on such characteristics. Characters based on internal skeletal structures are less prone to loss due to taphonomic processes. Previous work has shown that internal corallite architectural features measured in transverse thin sections can be used to distinguish species. To determine whether internal colony-level features measured on X-radiographs can be used. eight characters related to corallite budding and accretionary growth were measured on specimens representing three modern members of the M. annularis species complex (M. annularis, M. flaveolata and M. franksi), as well as two fossil forms (columnar and organ-pipe). All eight characters showed significant differences among species. Discriminant function analysis using seven of these characters resulted in distinct species groupings In canonical scores plots and a 100% classification success for specimens from Panama. These results suggest that measurements made on X-radiographs provide a useful tool for quantitatively distinguishing members of the M. annularis complex as well as between other massive reef corals.

Geographic variation in vocalisations of Satin Bowerbirds, Ptilonorynchus violaceus, in south-eastern Queensland

Geographic variation in the advertisement call of the male Satin Bowerbird, Ptilonorhynchus violaceus, was investigated in three populations in south-eastern Queensland. The call was found to differ significantly among the three geographically distinct populations. A discriminant function analysis using five measurements of call frequency and duration provided 100% classification success of the 25 individuals. The observed geographic variation in this call may result from adaptation to the local acoustic environment in these populations, or from genetic or cultural divergence among populations. Further research involving the acoustic properties of the habitats, population genetics and a larger number of populations is required to fully understand this pattern of call variation.

Multiparametric evaluation of dysphonic severity

In recent years, the multiparametric approach for evaluating perceptual rating of voice quality has been advocated. This study evaluates the accuracy of predicting perceived overall severity of voice quality with a minimal set of aerodynamic, voice range profile (phonetogram), and acoustic perturbation measures. One hundred and twelve dysphonic persons (93 women and 19 men) with laryngeal pathologies and 41 normal controls (35 women and six men) with normal voices participated in this study. Perceptual severity judgement was carried out by four listeners rating the G (overall grade) parameter of the GRBAS scale.(1) The minimal set of instrumental measures was selected based on the ability of the measure to discriminate between dysphonic and normal voices, and to attain at least a moderate correlation with perceived overall severity. Results indicated that perceived overall severity was best described by maximum phonation time of sustained /a/, peak intraoral pressure of the consonant-vowel /pi/ strings production, voice range profile area, and acoustic jitter. Direct-entry discriminant function analysis revealed that these four voice measures in combination correctly predicted 67.3% of perceived overall severity levels.

Psychometric assessment of recovery from mild traumatic brain injury

The Rapid Screen of Concussion (RSC) is a brief psychometric test battery, designed to provide a functional criterion to aid clinical diagnosis of mild traumatic brain injury (mTBI). The present research aimed to examine the utility of this instrument for assessing recovery after mTBI. Three studies were conducted. In Study 1, Discriminant Function Analysis was performed to determine how well the RSC differentiated uninjured controls (N¼16), from mTBI patients (N¼22) and moderate to severe TBI patients (N¼14), several months post-injury. As predicted, moderate to severe TBI patients achieved lower scores than the mTBI and control groups. The RSC also successfully differentiated between each of the diagnostic groups, yielding an overall correct classification rate of 75%. Study 2 examined the predictive utility of the RSC in the mTBI sample (N¼22). Acute injury performance on the RSC was correlated with post-injury scores at an average of 5.5 months post-injury. Statistically significant partial correlation coefficients (r¼0.53r¼0.80) were found for each of the subtests, showing that low acute RSC scores were predictive of poor recovery scores on the RSC after mTBI. In the third study, Reliable Change Indices were calculated on the RSC subtests to examine individual patterns of recovery from mTBI. While 17 of the 23 participants made a significant improvement on their acute injury DSST scores (74%), only 13 of 25 made a significant improvement on the Rapid Sentence Judgement Test (52%), highlighting differential recovery of function, and challenging the notion of full recovery from mTBI within 3 months. These overall results offer support for the construct and predictive validity of the RSC and demonstrate that inexpensive tests of brain function may be useful for managing mTBI acutely for prognosis.

The critical power function is dependent on the duration of the predictive exercise tests chosen

The linear relationship between work accomplished (W-lim) and time to exhaustion (t(lim)) can be described by the equation: W-lim = a + CP.t(lim). Critical power (CP) is the slope of this line and is thought to represent a maximum rate of ATP synthesis without exhaustion, presumably an inherent characteristic of the aerobic energy system. The present investigation determined whether the choice of predictive tests would elicit significant differences in the estimated CP. Ten female physical education students completed, in random order and on consecutive days, five art-out predictive tests at preselected constant-power outputs. Predictive tests were performed on an electrically-braked cycle ergometer and power loadings were individually chosen so as to induce fatigue within approximately 1-10 mins. CP was derived by fitting the linear W-lim-t(lim) regression and calculated three ways: 1) using the first, third and fifth W-lim-t(lim) coordinates (I-135), 2) using coordinates from the three highest power outputs (I-123; mean t(lim) = 68-193 s) and 3) using coordinates from the lowest power outputs (I-345; mean t(lim) = 193-485 s). Repeated measures ANOVA revealed that CPI123 (201.0 +/- 37.9W) > CPI135 (176.1 +/- 27.6W) > CPI345 (164.0 +/- 22.8W) (P < 0.05). When the three sets of data were used to fit the hyperbolic Power-t(lim) regression, statistically significant differences between each CP were also found (P < 0.05). The shorter the predictive trials, the greater the slope of the W-lim-t(lim) regression; possibly because of the greater influence of 'aerobic inertia' on these trials. This may explain why CP has failed to represent a maximal, sustainable work rate. The present findings suggest that if CP is to represent the highest power output that an individual can maintain for a very long time without fatigue then CP should be calculated over a range of predictive tests in which the influence of aerobic inertia is minimised.

Dormancy release in Lolium rigidum seeds is a function of thermal after-ripening time and seed water content

Dormancy release in seeds of Lolium rigidum Gaud. (annual ryegrass) was investigated in relation to temperature and seed water content. Freshly matured seeds were collected from cropping fields at Wongan Hills and Merredin, Western Australia. Seeds from Wongan Hills were equilibrated to water contents between 6 and 18% dry weight and after-ripened at constant temperatures between 9 and 50degreesC for up to 23 weeks. Wongan Hills and Merredin seeds at water contents between 7 and 17% were also after-ripened in full sun or shade conditions. Dormancy was tested at regular intervals during after-ripening by germinating seeds on agar at 12-h alternating 15degreesC (dark) and 25degreesC (light) periods. Rate of dormancy release for Wongan Hills seeds was a positive linear function of after-ripening temperature above a base temperature (T-b) of 5.4degreesC. A thermal after-ripening time model for dormancy loss accounting for seed moisture in the range 6-18% was developed using germination data for Wongan Hills seeds after-ripened at constant temperatures. The model accurately predicted dormancy release for Wongan Hills seeds after-ripened under naturally fluctuating temperatures. Seeds from Merredin responded similarly but had lower dormancy at collection and a faster rate of dormancy release in seeds below 9% water content.

Genetic association and cellular function of MC1R variant alleles in human pigmentation

We have examined MC1R variant allele frequencies in the general population of South East Queensland and in a collection of adolescent dizygotic and monozygotic twins and family members to define statistical associations with hair and skin color, freckling, and mole count. Results of these studies are consistent with a linear recessive allelic model with multiplicative penetrance in the inheritance of red hair. Four alleles, D84E, R151C, R160W, and D294H, are strongly associated with red hair and fair skin with multinomial regression analysis showing odds ratios of 63, 118, 50, and 94, respectively. An additional three low-penetrance alleles V60L, V92M, and R163Q have odds ratios 6, 5, and 2 relative to the wild-type allele. To address the cellular effects of MC1R variant alleles in signal transduction, we expressed these receptors in permanently transfected HEK293 cells. Measurement of receptor activity via induction of a cAMP-responsive luciferase reporter gene found that the R151C and R160W receptors were active in the presence of NDP-MSH ligand, but at much reduced levels compared with that seen with the wild-type receptor. The ability to stimulate phosphorylation of the cAMP response element binding protein (CREB) transcription factor was also apparent in all stimulated MC1R variant allele-expressing HEK293 cell extracts as assessed by immunoblotting. In contrast, human melanoma cell lines showed wide variation in the their ability to undergo cAMP-mediated CREB phosphorylation. Culture of human melanocytes of known MC1R genotype may provide the best experimental approach to examine the functional consequences for each MC1R variant allele. With this objective, we have established more than 300 melanocyte cell strains of defined MC1R genotype.