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.

Kinematic and non-linear analysis of foldable barrel vaults

Kinematic analysis is conducted to derive the geometric constraints for the geometric design of foldable barrel vaults (FBV) composed of polar or angulated scissor units. Non-linear structural analysis is followed to determine the structural response of FBVs in the fully deployed configuration under static loading. Two load cases are considered: cross wind and longitudinal wind. The effect of varying member sizes, depth-to-span ratio and geometric imperfections is examined. (C) 2000 Elsevier Science Ltd. All rights reserved.

Analysis and implementation of a new constant acceleration subcycling algorithm

An algorithm for explicit integration of structural dynamics problems with multiple time steps is proposed that averages accelerations to obtain subcycle states at a nodal interface between regions integrated with different time steps. With integer time step ratios, the resulting subcycle updates at the interface sum to give the same effect as a central difference update over a major cycle. The algorithm is shown to have good accuracy, and stability properties in linear elastic analysis similar to those of constant velocity subcycling algorithms. The implementation of a generalised form of the algorithm with non-integer time step ratios is presented. (C) 1997 by John Wiley & Sons, Ltd.

Monitoring the composition and form of urban environments based on the vegetation-impervious surface-soil (VIS) model by sub-pixel analysis techniques

The majority of the world's population now resides in urban environments and information on the internal composition and dynamics of these environments is essential to enable preservation of certain standards of living. Remotely sensed data, especially the global coverage of moderate spatial resolution satellites such as Landsat, Indian Resource Satellite and Systeme Pour I'Observation de la Terre (SPOT), offer a highly useful data source for mapping the composition of these cities and examining their changes over time. The utility and range of applications for remotely sensed data in urban environments could be improved with a more appropriate conceptual model relating urban environments to the sampling resolutions of imaging sensors and processing routines. Hence, the aim of this work was to take the Vegetation-Impervious surface-Soil (VIS) model of urban composition and match it with the most appropriate image processing methodology to deliver information on VIS composition for urban environments. Several approaches were evaluated for mapping the urban composition of Brisbane city (south-cast Queensland, Australia) using Landsat 5 Thematic Mapper data and 1:5000 aerial photographs. The methods evaluated were: image classification; interpretation of aerial photographs; and constrained linear mixture analysis. Over 900 reference sample points on four transects were extracted from the aerial photographs and used as a basis to check output of the classification and mixture analysis. Distinctive zonations of VIS related to urban composition were found in the per-pixel classification and aggregated air-photo interpretation; however, significant spectral confusion also resulted between classes. In contrast, the VIS fraction images produced from the mixture analysis enabled distinctive densities of commercial, industrial and residential zones within the city to be clearly defined, based on their relative amount of vegetation cover. The soil fraction image served as an index for areas being (re)developed. The logical match of a low (L)-resolution, spectral mixture analysis approach with the moderate spatial resolution image data, ensured the processing model matched the spectrally heterogeneous nature of the urban environments at the scale of Landsat Thematic Mapper data.

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.

Computational properties of min/max autocorrelation factors

Minimum/maximum autocorrelation factor (MAF) is a suitable algorithm for orthogonalization of a vector random field. Orthogonalization avoids the use of multivariate geostatistics during joint stochastic modeling of geological attributes. This manuscript demonstrates in a practical way that computation of MAF is the same as discriminant analysis of the nested structures. Mathematica software is used to illustrate MAF calculations from a linear model of coregionalization (LMC) model. The limitation of two nested structures in the LMC for MAF is also discussed and linked to the effects of anisotropy and support. The analysis elucidates the matrix properties behind the approach and clarifies relationships that may be useful for model-based approaches. (C) 2003 Elsevier Science Ltd. All rights reserved.

Overlapping species boundaries and hybridization within the Monastraea "annularis" reef coral complex in the Pleistocene of the Bahama Islands

Recent molecular analyses indicate that many reef coral species belong to hybridizing species complexes or "syngameons." Such complexes consist of numerous genetically distinct-species or lineages, which periodically split and/or fuse as they extend through time. During splitting and fusion, morphologic intermediates form and species overlap. Here we focus on processes associated with lineage fusion, specifically introgressive hybridization, and the recognition of such hybridization in the fossil record. Our approach involves comparing patterns of ecologic and morphologic overlap in genetically characterized modern species with fossil representatives of the same or closely related species. We similarly consider the long-term consequences of past hybridization on the structure of modern-day species boundaries. Our study involves the species complex Montastraea annularis s.l. and is based in the Bahamas, where, unlike other Caribbean locations, two of the three members of the complex today are not genetically distinct. We measured and collected colonies along linear transects across Pleistocene reef terraces of last interglacial age (approximately 125 Ka) on the islands of San Salvador, Andros, and Great Inagua. We performed quantitative ecologic and morphologic analyses of the fossil data, and compared patterns of overlap among species with data from modern localities where species are and are not genetically distinct. Ecologic and morphologic analyses reveal "moderate" overlap (>10%, but statistically significant differences) and sometimes "high" overlap (no statistically significant differences) among Pleistocene growth forms (= "species"). Ecologic analyses show that three species (massive, column, organ-pipe) co-occurred. Although organ-pipes had higher abundances in patch reef environments, columnar and massive species exhibited broad, completely overlapping distributions and had abundances that were not related to reef environment. For morphometric analyses, we used multivariate discriminant analysis on landmark data and linear measurements. The results show that columnar species overlap "moderately" with organ-pipe and massive species. Comparisons with genetically characterized colonies from Panama show that the Pleistocene Bahamas species have intermediate morphologies, and that the observed "moderate" overlap differs from the morphologic separation among the three modern species. In contrast, massive and columnar species from the Pleistocene of the Dominican Republic comprise distinct morphologic clusters, similar to the modern species; organ-pipe species exhibit "low" overlap (

Dependence of public support for survival of wildlife species on their likeability

We surveyed a sample of 204 individuals selected from the public in Brisbane, Australia, to ascertain the extent to which they like or dislike 24 species of wildlife present in tropical Australia. The species belong to three classes: mammals, birds and reptiles. We calculated likeability indices for each of these species. We also asked respondents if they favoured the survival of each of these species and so the percentage of respondents favouring survival of each of these species could be calculated. Thus, using linear regression analysis, the percentage of respondents favouring survival of each of the species was related to their indices of likeability. In addition, the data enables the average likeability of species in the three classes (mammals, birds and reptiles) to be compared with the average support for survival of species in each of these three classes. As a result, we are able to assess how important stated likeability seems to be for preferences for survival of species, and to reconsider the hypothesis in the literature that there is likely to be more public support for the survival of mammals than for birds than for reptiles.

Increased hepatic iron concentration in nonalcoholic steatohepatitis is associated with increased fibrosis

Background & Aims: Nonalcoholic steatohepatitis (NASH) is a chronic liver disease that occasionally progresses to cirrhosis but usually has a benign course. The aim of this study was to investigate the role of the hemochromatosis mutation Cys282Tyr in development of the mild hepatic iron overload found in some patients with NASH and its association with hepatic damage in these patients. Methods: Fifty-one patients with NASH were studied. The presence of the Cys282Tyr mutation was tested in all patients, and the data were analyzed with respect to the histological grade of steatosis, inflammation, Perls' staining, hepatic iron concentration (HIC), and serum iron indices. Results: Thirty-one percent of patients with NASH were either homozygous or heterozygous for the Cys282Tyr mutation. This mutation was significantly associated with Perls' stain grade (P < 0.005), HIC (P < 0.005), and transferrin saturation percentage (P < 0.005) but not with serum ferritin levels. Linear regression analysis showed that increased hepatic iron (Perls' stain or HIC) had the greatest association with the severity of fibrosis (P < 0.0001). Conclusions: The Cys282Tyr mutation is responsible for most of the mild iron overload found in NASH and thus has a significant association with hepatic damage in these patients. Heterozygosity for the hemochromatosis gene mutation therefore cannot always be considered benign.

Referral to post-acute care following traumatic brain injury (TBI) in the Australian context

The study aimed to describe the types of care allocated at the end of acute care to people diagnosed with TBI and to identify the factors associated with variations in referral to care. A retrospective analysis of medical records of 61 patients was conducted based on a sample from two hospitals. While 60.7% of the study sample were referred to formal rehabilitation care, this was primarily non-inpatient rehabilitation care (32.8%). Discriminant analysis was used to determine medical and non-medical predictors of referral. Results indicated that place of treatment and age contribute to group differences and were significant in separating the inpatient rehabilitation group from the non-inpatient and no rehabilitation groups. Review by a rehabilitation physician was associated with referral to inpatient rehabilitation but was not adequate to explain referral to non-inpatient rehabilitation. An in-depth exploration of post-acute referral is warranted to improve policy and practice in relation to continuity of care following TBI.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

Weight loss during prolonged lactation in rural Bangladeshi mothers

To determine the duration of lactation which is associated with weight loss in rural Bangladeshi mothers and also to determine the relationship with consumption patterns of principal food items, a cross-sectional study was carried out among 791 lactating rural Bangladeshi mothers aged 18-40 years. Results were compared with 333 non-pregnant and non-lactating mothers of a similar age group. The duration of lactation was up to 60 months. The mean difference in body-weight and body mass index (BMI) of lactating mothers who breastfed their children up to 24 months was significantly lower compared to non-lactating mothers of the same age group, but no differences were observed for those who breastfed beyond 24 months. The frequency of consumption of principal food items was comparable between the non-lactating and the lactating mothers who breastfed beyond 24 months. Results of multiple linear regression analysis showed that body-weight of mothers was negatively correlated with 1-12 month(s) and 13-24 months of lactation after controlling for height, education, and food consumption (slope -1.04, p < 0.05 and slope -1.23, p < 0.05 respectively). Height and consumption of meat and milk were significantly positively correlated with body-weight (slope 0.53, p < 0.001; slope 1.44, p < 0.001; and slope 0.75, p < 0.05 respectively). The study concluded that Bangladeshi women who breastfed up to 24 months were of lower weight than non-lactating mothers, most likely due to the effect of lactation. These mothers were not taking any additional foods during their lactating period. Based on the findings of the study, it is recommended that mothers consume additional energy-rich foods during the first 24 months of lactation to prevent weight loss.

Evaluation of bioelectrical impedance for prospective nutritional assessment in cystic fibrosis

We have compared the use of bioelectrical impedance analysis (BIA) with anthropometry for the prediction of changes in total body potassium (TBK) in a group (n = 31) of children with cystic fibrosis. Linear regression analysis showed that TBK was highly correlated (r > 0.93) with height(2)/impedance, weight, height, and fat-free mass (FFM) estimated from skin-fold measurements. Changes in TBK were also correlated, but less well, with changes in height(2)/impedance, weight, height, and FFM (r = 0.69, 0.59, 0.44, and 0.40, respectively). The children were divided into two groups: those who had normal accretion of TBK (> 5%/y) and those who had suboptimal accretion of TBK (< 5%/y). Analysis of variance showed that the significant difference in the change in TBK between the groups was detectable by concomitant changes in impedance and weight but not by changes in height, FFM, or weight and height Z scores. The results of this study suggest that serial BIA measures may be useful as a predictor of progressive undernutrition and poor growth in children with cystic fibrosis. (C) Elsevier Science Inc. 1997.