A mixture model-based approach to the clustering of microarray expression data

Motivation: This paper introduces the software EMMIX-GENE that has been developed for the specific purpose of a model-based approach to the clustering of microarray expression data, in particular, of tissue samples on a very large number of genes. The latter is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. A feasible approach is provided by first selecting a subset of the genes relevant for the clustering of the tissue samples by fitting mixtures of t distributions to rank the genes in order of increasing size of the likelihood ratio statistic for the test of one versus two components in the mixture model. The imposition of a threshold on the likelihood ratio statistic used in conjunction with a threshold on the size of a cluster allows the selection of a relevant set of genes. However, even this reduced set of genes will usually be too large for a normal mixture model to be fitted directly to the tissues, and so the use of mixtures of factor analyzers is exploited to reduce effectively the dimension of the feature space of genes. Results: The usefulness of the EMMIX-GENE approach for the clustering of tissue samples is demonstrated on two well-known data sets on colon and leukaemia tissues. For both data sets, relevant subsets of the genes are able to be selected that reveal interesting clusterings of the tissues that are either consistent with the external classification of the tissues or with background and biological knowledge of these sets.

Estimating Two-Stage Models for Genetic Influences on Alcohol, Tobacco or Drug Use Initiation and Dependence Vulnerability in Twin and Family Data

Genetic research on risk of alcohol, tobacco or drug dependence must make allowance for the partial overlap of risk-factors for initiation of use, and risk-factors for dependence or other outcomes in users. Except in the extreme cases where genetic and environmental risk-factors for initiation and dependence overlap completely or are uncorrelated, there is no consensus about how best to estimate the magnitude of genetic or environmental correlations between Initiation and Dependence in twin and family data. We explore by computer simulation the biases to estimates of genetic and environmental parameters caused by model misspecification when Initiation can only be defined as a binary variable. For plausible simulated parameter values, the two-stage genetic models that we consider yield estimates of genetic and environmental variances for Dependence that, although biased, are not very discrepant from the true values. However, estimates of genetic (or environmental) correlations between Initiation and Dependence may be seriously biased, and may differ markedly under different two-stage models. Such estimates may have little credibility unless external data favor selection of one particular model. These problems can be avoided if Initiation can be assessed as a multiple-category variable (e.g. never versus early-onset versus later onset user), with at least two categories measurable in users at risk for dependence. Under these conditions, under certain distributional assumptions., recovery of simulated genetic and environmental correlations becomes possible, Illustrative application of the model to Australian twin data on smoking confirmed substantial heritability of smoking persistence (42%) with minimal overlap with genetic influences on initiation.

Stream-bank shade and larval distribution of the Philippine malaria vector Anopheles flavirostris

The principal malaria vector in the Philippines, Anopheles flavirostris (Ludlow) (Diptera: Culicidae), is regarded as 'shade-loving' for its breeding sites, i.e. larval habitats. This long-standing belief, based on circumstantial observations rather than ecological analysis, has guided larval control methods such as 'stream-clearing' or the removal of riparian vegetation, to reduce the local abundance of An. flavirostris . We measured the distribution and abundance of An. flavirostris larvae in relation to canopy vegetation cover along a stream in Quezon Province, the Philippines. Estimates of canopy openness and light measurements were obtained by an approximation method that used simplified assumptions about the sun, and by hemispherical photographs analysed using the program hemiphot(C) . The location of larvae, shade and other landscape features was incorporated into a geographical information system (GIS) analysis. Early larval instars of An. flavirostris were found to be clustered and more often present in shadier sites, whereas abundance was higher in sunnier sites. For later instars, distribution was more evenly dispersed and only weakly related to shade. The best predictor of late-instar larvae was the density of early instars. Distribution and abundance of larvae were related over time (24 days). This pattern indicates favoured areas for oviposition and adult emergence, and may be predictable. Canopy measurements by the approximation method correlated better with larval abundance than hemispherical photography, being economical and practical for field use. Whereas shade or shade-related factors apparently have effects on larval distribution of An. flavirostris , they do not explain it completely. Until more is known about the bionomics of this vector and the efficacy and environmental effects of stream-clearing, we recommend caution in the use of this larval control method.

Testing insect movement theory: empirical analysis of pest data routinely collected from agricultural crops

Observations of an insect's movement lead to theory on the insect's flight behaviour and the role of movement in the species' population dynamics. This theory leads to predictions of the way the population changes in time under different conditions. If a hypothesis on movement predicts a specific change in the population, then the hypothesis can be tested against observations of population change. Routine pest monitoring of agricultural crops provides a convenient source of data for studying movement into a region and among fields within a region. Examples of the use of statistical and computational methods for testing hypotheses with such data are presented. The types of questions that can be addressed with these methods and the limitations of pest monitoring data when used for this purpose are discussed. (C) 2002 Elsevier Science B.V. All rights reserved.

A hedging point strategy - balancing effluent quality, economy and robustness in the control of wastewater treatment plants

An operational space map is an efficient tool to compare a large number of operational strategies to find an optimal choice of setpoints based on a multicriterion. Typically, such a multicriterion includes a weighted sum of cost of operation and effluent quality. Due to the relative high cost of aeration such a definition of optimality result in a relatively high fraction of the effluent total nitrogen in the form of ammonium. Such a strategy may however introduce a risk into operation because a low degree of ammonium removal leads to a low amount of nitrifiers. This in turn leads to a reduced ability to reject event disturbances, such as large variations in the ammonium load, drop in temperature, the presence of toxic/inhibitory compounds in the influent etc. Hedging is a risk minimisation tool, with the aim to "reduce one's risk of loss on a bet or speculation by compensating transactions on the other side" (The Concise Oxford Dictionary (1995)). In wastewater treatment plant operation hedging can be applied by choosing a higher level of ammonium removal to increase the amount of nitrifiers. This is a sensible way to introduce disturbance rejection ability into the multi criterion. In practice, this is done by deciding upon an internal effluent ammonium criterion. In some countries such as Germany, a separate criterion already applies to the level of ammonium in the effluent. However, in most countries the effluent criterion applies to total nitrogen only. In these cases, an internal effluent ammonium criterion should be selected in order to secure proper disturbance rejection ability.

Vertical fluxes of sediment in oscillatory sheet flow

Time series of vertical sediment fluxes are derived from concentration time series in sheet flow under waves. While the concentrations C(z,t) vary very little with time for \z\ < 10d(50), the measured vertical sediment fluxes Q(zs)(z,t) vary strongly with time in this vertical band and their time variation follows, to some extent, the variation of the grain roughness Shields parameter 02,5(t). Thus, sediment distribution models based on the pickup function boundary condition are in some qualitative agreement with the measurements. However, the pickup function models are only able to model the upward bursts of sediment during the accelerating phases of the flow. They are, so far, unable to model the following strong downward sediment fluxes, which are observed during the periods of flow deceleration. Classical pickup functions, which essentially depend on the Shields parameter, are also incapable of modelling the secondary entrainment fluxes, which sometimes occur at free stream velocity reversal. The measured vertical fluxes indicate that the effective sediment settling velocity in the high [(0.3 < C(z,t) < 0.4] concentration area is typically only a few percent of the clear water settling velocity, while the measurements of Richardson and Jeronimo [Chem. Eng. Sci. 34 (1979) 1419], from a different physical setting, lead to estimates of the order 20%. The data does not support gradient diffusion as a model for sediment entrainment from the bed. That is, detailed modelling of the observed near-bed fluxes would require diffusivities that go negative during periods of flow deceleration. An observed general trend for concentration variability to increase with elevation close to the bed is also irreconcilable with diffusion models driven by a bottom boundary condition. (C) 2002 Published by Elsevier Science B.V.

Longitudinal study of veterinarians from entry to the veterinary course to 10 years after graduation: attitudes to work, career and profession

Objective To describe the attitudes of veterinarians to their work, career and profession during the 10 years after graduation. Design Longitudinal study of students who started their course at The University of Queensland in 1985 and 1986, and who completed questionnaires in their first and fifth year as students, and after one, five and 10 years as veterinarians. Methods Data from 129 (96%) questionnaires completed after 10 years as a veterinarian were coded numerically then analysed, together with data from previous questionnaires, with SAS System 7 for Windows 95. Results After 10 years, almost all respondents were either very glad they had done the veterinary course (57%) or generally glad, though with some misgivings (37%). Despite this, only 55% would definitely become a veterinarian if they 'had to do it over again'. The responses for about one-third were different from those given five years earlier. The views of many were related to the level of support and encouragement received in their first job after graduation. There were 42% who were working less than half-time as veterinarians, and their main reasons were, in order, raising children, long hours of work, attitudes of bosses and clients, and poor pay. A majority was concerned about the ethics and competence of some colleagues, and almost all believed that consideration of costs must influence the type of treatment animals receive. Conclusions Most veterinarians were glad to have done the veterinary course, but for about one-quarter their career had not lived up to expectations and almost half would not do it again in another incarnation. Stress, hours of work, difficulties in balancing personal life with career and low income were important concerns for many. Low income may contribute to the low number of males entering the veterinary profession.

Phase equilibrium data and liquidus for the system "MnO"-CaO-(Al2O3+SiO2) at Al2O3/SiO2 of 0.41

Phase-equilibrium data and the liquidus for the system. "MnO"-CaO-(Al2O3-SiO2) at a manganese-rich alloy saturation have been determined in the temperature range from 1423 to 1723 K. The results are presented in the form of a pseudoternary section "MnO"-CaO-(Al2O3 + SiO2) with an Al2O3/SiO2 weight ratio of 0.41. The following primary phases are present in the range of conditions investigated:, 3Al(2)O(3).2SiO(2); SiO2; MnO.Al2O3-2SiO(2); (Mn,Ca)O.SiO2; 2(Mn,Ca)O.SiO2; MnO.Al2O3; (Mn,Ca)O; alpha-2CaO.SiO2; alpha'-2CaO.SiO2; 2CaO.Al2O3.SiO2; CaO.SiO2, and CaO.Al2O3.2SiO(2). The presence of alumina in this system is shown to have a significant effect on the liquidus compared to the system "MnO"-CaO-SiO2, leading to, the stabilization of the anorthite and gehlenite phases.

Balancing acts: Respiratory sensations, motor control and human posture

1. The present brief review covers some novel aspects of integration between respiration and movement of the body. 2. There are potent viscerosomatic reflexes in animals involving small-diameter pulmonary afferents that, when excited, would limit exercise. However, recent studies using lobeline injections to excite pulmonary afferents in awake humans suggest that there is no evoked reflex motoneuronal inhibition. Instead, the noxious respiratory sensations generated by the vagal afferents may be crucial in the decision to stop exercise. 3. While respiratory movements may affect limb movements, the control of the trunk and limbs can involve interaction (and even interference) with key respiratory muscles, such as the diaphragm. Recent studies have revealed that not only does the diaphragm receive feed-forward drive prior to some limb movements, but that it also contracts both phasically and tonically during repetitive limb movements. 4. Thus, challenges to posture can indirectly challenge ventilation, while coordinated diaphragm contraction may contribute to control of the trunk.

Mineral chemistry, whole-rock compositions, and petrogenesis of leg 176 gabbros: Data and discussion

We report mineral chemistry, whole-rock major element compositions, and trace element analyses on Hole 735B samples drilled and selected during Leg 176. We discuss these data, together with Leg 176 shipboard data and Leg 118 sample data from the literature, in terms of primary igneous petrogenesis. Despite mineral compositional variation in a given sample, major constituent minerals in Hole 735B gabbroic rocks display good chemical equilibrium as shown by significant correlations among Mg# (= Mg/[Mg+Fe2+]) of olivine, clinopyroxene, and orthopyroxene and An (=Ca/[Ca+Na]) of plagioclase. This indicates that the mineral assemblages olivine + plagioclase in troctolite, plagioclase + clinopyroxene in gabbro, plagioclases + clinopyroxene + olivine in olivine gabbro, and plagioclase + clinopyroxene + olivine + orthopyroxene in gabbronorite, and so on, have all coprecipitated from their respective parental melts. Fe-Ti oxides (ilmenite and titanomagnetite), which are ubiquitous in most of these rocks, are not in chemical equilibrium with olivine, clinopyroxene, and plagioclase, but precipitated later at lower temperatures. Disseminated oxides in some samples may have precipitated from trapped Fe-Ti–rich melts. Oxides that concentrate along shear bands/zones may mark zones of melt coalescence/transport expelled from the cumulate sequence as a result of compaction or filter pressing. Bulk Hole 735B is of cumulate composition. The most primitive olivine, with Fo = 0.842, in Hole 735B suggests that the most primitive melt parental to Hole 735B lithologies must have Mg# ≤ 0.637, which is significantly less than Mg# = 0.714 of bulk Hole 735B.