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.

Combining the analyses of introgressive hybridisation and linkage mapping to investigate the genetic architecture of population divergence in the lake whitefish (Coregonus clupeaformis, Mitchill)

Adaptation and reproductive isolation, the engines of biological diversity, are still elusive when discussing the genetic bases of speciation. Namely, the number of genes and magnitude of selection acting positively or negatively on genomic traits implicated in speciation is contentious. Here, we describe the first steps of an ongoing research program aimed at understanding the genetic bases of population divergence and reproductive isolation in the lake whitefish (Coregonus clupeaformis). A preliminary linkage map originating from a hybrid cross between dwarf and normal ecotypes is presented, whereby some of the segregating AFLP markers were found to be conserved among natural populations. Maximum-likelihood was used to estimate hybrid indices from non-diagnostic markers at 998 AFLP loci. This allowed identification of the most likely candidate loci that have been under the influence of selection during the natural hybridisation of whitefish originating from different glacial races. As some of these loci could be identified on the linkage map, the possibility that selection of traits in natural populations may eventually be correlated to specific chromosomal regions was demonstrated. The future prospects and potential of these approaches to elucidate the genetic bases of adaptation and reproductive isolation among sympatric ecotypes of lake whitefish is discussed.

Nonverbal narratives: Listening to people with severe intellectual disability

This article describes an exploratory study that examined the perspectives of practitioners who spend much of their working day listening to and in some ways interpreting for people with severe intellectual disabilities. On the basis of focus group interviews with 23 professional disability-sector workers, including speech therapists, psychologists, and human service workers, the article reports on the importance of a practitioner's values and experience in successful interactions with individuals who rely on self-developed nonsymbolic communication repertoires. The article includes a discussion of the likelihood of including individuals with severe intellectual disabilities in narrative research.

Outcome of stroke patients admitted to intensive care: Experience from an Australian teaching hospital

The objective of this study was to determine the mortality rate and the functional outcomes of stroke patients admitted to the intensive care unit (ICU) and to identify predictors of poor outcome in this population. The records of all patients admitted to the ICU with the diagnosis of stroke between January 1994 and December 1999 were reviewed. Patients with subarachnoid haemorrhage were excluded. Data were collected on clinical and biological variables, risk factors for stroke and the presence of comorbidities. Mortality (ICU, in-hospital and three-month) and functional outcome were used as end-points. In the six-year-period, 61 patients were admitted to the ICU with either haemorrhagic or ischaemic stroke. Medical records were available for only 58 patients. There were 23 ischaemic and 35 haemorrhagic strokes. The ICU, in-hospital and three-month mortality rates were 36%, 47% and 52% respectively. There were no significant differences in the prevalence of premorbid risk factors between survivors and non-survivors. The mean Barthel score was significantly different between the independent and dependent survivors (94 +/- 6 vs 45 +/- 26, P < 0.001). A substantial number of patients with good functional outcomes had lower Rankin scores (92% vs 11%, P < 0.001). Only 46% of those who were alive at three months were functionally independent. Intensive care admission was associated with a high mortality rate and a high likelihood of dependent lifestyle after hospital discharge. Haemorrhagic stroke, fixed dilated pupil(s) and GCS < 10 during assessment were associated with increased mortality and poor functional outcome.

Longitudinal study of veterinarians from entry to the veterinary course to ten years after graduation: career paths

Objective To map out the career paths of veterinarians during their first 10 years after graduation, and to determine if this could have been predicted at entry to the veterinary course. 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 in their second, sixth and eleventh year as veterinarians. Methods Data from 129 (96%) questionnaires completed during the eleventh year after graduation were coded numerically then analysed, together with data from previous questionnaires, with SAS System 7 for Windows 95. Results Ten years after they graduated, 80% were doing veterinary work, 60% were in private practice, 40% in small animal practice and 18% in mixed practice. The equivalent of 25% of the working time of all females was taken up by family duties. When part-time work was taken into account, veterinary work constituted the equivalent of 66% of the group working full-time. That 66% consisted of 52% on small animals, 7% on horses, 6% on cattle/sheep and 1% on pigs/poultry. Those who had grown up on farms with animals were twice as likely to be working with farm animals as were those from other backgrounds. Forecasts made on entry to the veterinary course were of no value in predicting who would remain in mixed practice. Conclusions Fewer than one-fifth of graduates were in mixed practice after 10 years, but the number was higher for those who grew up on farms with animals. Forecasts that may be made at interview before entry to the course were of little value in predicting the likelihood of remaining in mixed veterinary practice.

Simulation of insect movement with respect to plant architecture and morphogenesis

Developments in computer and three dimensional (3D) digitiser technologies have made it possible to keep track of the broad range of data required to simulate an insect moving around or over the highly heterogeneous habitat of a plant's surface. Properties of plant parts vary within a complex canopy architecture, and insect damage can induce further changes that affect an animal's movements, development and likelihood of survival. Models of plant architectural development based on Lindenmayer systems (L-systems) serve as dynamic platforms for simulation of insect movement, providing ail explicit model of the developing 3D structure of a plant as well as allowing physiological processes associated with plant growth and responses to damage to be described and Simulated. Simple examples of the use of the L-system formalism to model insect movement, operating Lit different spatial scales-from insects foraging on an individual plant to insects flying around plants in a field-are presented. Such models can be used to explore questions about the consequences of changes in environmental architecture and configuration on host finding, exploitation and its population consequences. In effect this model is a 'virtual ecosystem' laboratory to address local as well as landscape-level questions pertinent to plant-insect interactions, taking plant architecture into account. (C) 2002 Elsevier Science B.V. All rights reserved.

Remitting behaviour of Nepalese rural-to-urban migrants: Implications for theory and policy

This article uses data for Nepal to test contemporary hypotheses about the remitting behaviour and associated motives of rural-to-urban migrants and to consider the likely impact of such remittances on rural development. Possibilities for inheritance, degree of family attachment, likelihood of eventual return to place of origin and family investment in the education of the migrants are found to be significant influences on levels of remittances by Nepalese migrants. However, in Nepal, remittances do not seem to result in long-term capital investment in rural areas and so may not promote long-term development of these areas.

Using crop simulation to generate genotype by environment interaction effects for sorghum in water-limited environments

Multi-environment trials (METs) used to evaluate breeding lines vary in the number of years that they sample. We used a cropping systems model to simulate the target population of environments (TPE) for 6 locations over 108 years for 54 'near-isolines' of sorghum in north-eastern Australia. For a single reference genotype, each of 547 trials was clustered into 1 of 3 'drought environment types' (DETs) based on a seasonal water stress index. Within sequential METs of 2 years duration, the frequencies of these drought patterns often differed substantially from those derived for the entire TPE. This was reflected in variation in the mean yield of the reference genotype. For the TPE and for 2-year METs, restricted maximum likelihood methods were used to estimate components of genotypic and genotype by environment variance. These also varied substantially, although not in direct correlation with frequency of occurrence of different DETs over a 2-year period. Combined analysis over different numbers of seasons demonstrated the expected improvement in the correlation between MET estimates of genotype performance and the overall genotype averages as the number of seasons in the MET was increased.

An EM-based Semi-Parametric Mixture Model Approach to the Regression Analysis of Competing-Risks Data

We consider a mixture model approach to the regression analysis of competing-risks data. Attention is focused on inference concerning the effects of factors on both the probability of occurrence and the hazard rate conditional on each of the failure types. These two quantities are specified in the mixture model using the logistic model and the proportional hazards model, respectively. We propose a semi-parametric mixture method to estimate the logistic and regression coefficients jointly, whereby the component-baseline hazard functions are completely unspecified. Estimation is based on maximum likelihood on the basis of the full likelihood, implemented via an expectation-conditional maximization (ECM) algorithm. Simulation studies are performed to compare the performance of the proposed semi-parametric method with a fully parametric mixture approach. The results show that when the component-baseline hazard is monotonic increasing, the semi-parametric and fully parametric mixture approaches are comparable for mildly and moderately censored samples. When the component-baseline hazard is not monotonic increasing, the semi-parametric method consistently provides less biased estimates than a fully parametric approach and is comparable in efficiency in the estimation of the parameters for all levels of censoring. The methods are illustrated using a real data set of prostate cancer patients treated with different dosages of the drug diethylstilbestrol. Copyright (C) 2003 John Wiley Sons, Ltd.

The validity of publication and citation counts for Sociology and other selected disciplines

The use of bibliometric data is a means of comparing. research productivity and scholarly. impact for individuals, work groups, institutions and nations within and between disciplines. Central to this debate is the notion that disciplines differ in the ways in which,they exchange ideas and disseminate information and therefore have diverse publishing and citation patterns. In this article we use two different approaches to compiling bibliometric data to compare publishing patterns of five different disciplines that encompass Molecular Biology; Administration/Political Science, Psychology,. Philosophy and Sociology/Anthropology. We find that the social sciences differ from each other as well as from the physical sciences in their publication and citation patterns. Further, while the different ways of organizing the data produce somewhat different results, the substantive findings for the general patterning of publications and citations of disciplines are consistent for both data sets. Sociology/Anthropology, when compared with the other disciplines, shows substantial differences across universities.

Centile reference for total energy expenditure in infants from 1 to 12 months

Background: A knowledge of energy expenditure in infancy is required for the estimation of recommended daily amounts of food energy, for designing artificial infant feeds, and as a reference standard for studies of energy metabolism in disease states. Objectives: The objectives of this study were to construct centile reference charts for total energy expenditure (TEE) in infants across the first year of life. Methods: Repeated measures of TEE using the doubly labeled water technique were made in 162 infants at 1.5, 3, 6, 9 and 12 months. In total, 322 TEE measurements were obtained. The LMS method with maximum penalized likelihood was used to construct the centile reference charts. Centiles were constructed for TEE expressed as MJ/day and also expressed relative to body weight (BW) and fat-free mass (FFM). Results: TEE increased with age and was 1.40,1.86, 2.64, 3.07 and 3.65 MJ/day at 1.5, 3, 6, 9 and 12 months, respectively. The standard deviations were 0.43, 0.47, 0.52, 0.66 and 0.88, respectively. TEE in MJ/kg increased from 0.29 to 0.36 and in MJ/day/kg FFM from 0.36 to 0.48. Conclusions: We have presented centile reference charts for TEE expressed as MJ/day and expressed relative to BW and FFM in infants across the first year of life. There was a wide variation or biological scatter in TEE values seen at all ages. We suggest that these centile charts may be used to assess and possibly quantify abnormal energy metabolism in disease states in infants.

An integrated system to support electricity tariff contract definition

This paper presents an integrated system that helps both retail companies and electricity consumers on the definition of the best retail contracts and tariffs. This integrated system is composed by a Decision Support System (DSS) based on a Consumer Characterization Framework (CCF). The CCF is based on data mining techniques, applied to obtain useful knowledge about electricity consumers from large amounts of consumption data. This knowledge is acquired following an innovative and systematic approach able to identify different consumers’ classes, represented by a load profile, and its characterization using decision trees. The framework generates inputs to use in the knowledge base and in the database of the DSS. The rule sets derived from the decision trees are integrated in the knowledge base of the DSS. The load profiles together with the information about contracts and electricity prices form the database of the DSS. This DSS is able to perform the classification of different consumers, present its load profile and test different electricity tariffs and contracts. The final outputs of the DSS are a comparative economic analysis between different contracts and advice about the most economic contract to each consumer class. The presentation of the DSS is completed with an application example using a real data base of consumers from the Portuguese distribution company.

Introduction to the special issue on real world applications of intelligent tutoring systems

The idea behind creating this special issue on real world applications of intelligent tutoring systems was to bring together in a single publication some of the most important examples of success in the use of ITS technology. This will serve as a reference to all researchers working in the area. It will also be an important resource for the industry, showing the maturity of ITS technology and creating an atmosphere for funding new ITS projects. Simultaneously, it will be valuable to academic groups, motivating students for new ideas of ITS and promoting new academic research work in the area.

Ethical dilemmas in scientific publication: pitfalls and solutions for editors

Editors of scientific journals need to be conversant with the mechanisms by which scientific misconduct is amplified by publication practices. This paper provides definitions, ways to document the extent of the problem, and examples of editorial attempts to counter fraud. Fabrication, falsification, duplication, ghost authorship, gift authorship, lack of ethics approval, non-disclosure, 'salami' publication, conflicts of interest, auto-citation, duplicate submission, duplicate publications, and plagiarism are common problems. Editorial misconduct includes failure to observe due process, undue delay in reaching decisions and communicating these to authors, inappropriate review procedures, and confounding a journal's content with its advertising or promotional potential. Editors also can be admonished by their peers for failure to investigate suspected misconduct, failure to retract when indicated, and failure to abide voluntarily by the six main sources of relevant international guidelines on research, its reporting and editorial practice. Editors are in a good position to promulgate reasonable standards of practice, and can start by using consensus guidelines on publication ethics to state explicitly how their journals function. Reviewers, editors, authors and readers all then have a better chance to understand, and abide by, the rules of publishing.

The gambin model provides a superior fit to species abundance distributions with a single free parameter : evidence, implementation and interpretation

The species abundance distribution (SAD) has been a central focus of community ecology for over fifty years, and is currently the subject of widespread renewed interest. The gambin model has recently been proposed as a model that provides a superior fit to commonly preferred SAD models. It has also been argued that the model's single parameter (α) presents a potentially informative ecological diversity metric, because it summarises the shape of the SAD in a single number. Despite this potential, few empirical tests of the model have been undertaken, perhaps because the necessary methods and software for fitting the model have not existed. Here, we derive a maximum likelihood method to fit the model, and use it to undertake a comprehensive comparative analysis of the fit of the gambin model. The functions and computational code to fit the model are incorporated in a newly developed free-to-download R package (gambin). We test the gambin model using a variety of datasets and compare the fit of the gambin model to fits obtained using the Poisson lognormal, logseries and zero-sum multinomial distributions. We found that gambin almost universally provided a better fit to the data and that the fit was consistent for a variety of sample grain sizes. We demonstrate how α can be used to differentiate intelligibly between community structures of Azorean arthropods sampled in different land use types. We conclude that gambin presents a flexible model capable of fitting a wide variety of observed SAD data, while providing a useful index of SAD form in its single fitted parameter. As such, gambin has wide potential applicability in the study of SADs, and ecology more generally.