Using mixture models to detect differentially expressed genes

An important and common problem in microarray experiments is the detection of genes that are differentially expressed in a given number of classes. As this problem concerns the selection of significant genes from a large pool of candidate genes, it needs to be carried out within the framework of multiple hypothesis testing. In this paper, we focus on the use of mixture models to handle the multiplicity issue. With this approach, a measure of the local false discovery rate is provided for each gene, and it can be implemented so that the implied global false discovery rate is bounded as with the Benjamini-Hochberg methodology based on tail areas. The latter procedure is too conservative, unless it is modified according to the prior probability that a gene is not differentially expressed. An attractive feature of the mixture model approach is that it provides a framework for the estimation of this probability and its subsequent use in forming a decision rule. The rule can also be formed to take the false negative rate into account.

Assessing the effectiveness of a mammography screening service

Background: Trials have shown that mammography screening reduces mortality and probably decreases morbidity related to breast cancer. Methods: We assessed whether the major mammography service in Western Australia (BreastScreen WA) is likely to reduce mortality by comparing prognostic variables between screen-detected and other cases of breast cancer diagnosed in 1999. We assessed likely reductions in morbidity by comparing treatments received by these two groups. To confirm mortality and morbidity reduction, we also compared prognostic variables and treatments with targets. Information on demographic variables, tumour characteristics at presentation and treatments were collected from medical records for all incident cases of breast cancer in Western Australia in 1999. We matched cases with the Western Australian Cancer Registry records to determine which cases had been detected by BreastScreen WA. Results: BreastScreen WA achieved the targets for mortality reduction. Tumours detected by BreastScreen WA were smaller in size, less likely to have vascular invasion, of lower histological grade and were more likely to be ductal carcinoma in situ alone without invasive carcinoma. Oestrogen receptor status was more likely to be positive, the difference in progesterone status was not significant, and lymph node involvement tended to be lower. BreastScreen WA patients were treated more often with local therapy and less often with systemic therapy, and the proportion of patients treated with breast-conserving surgery was close to the target for minimizing morbidity in breast cancer. Conclusion: Mammographic detection of breast cancer by BreastScreen WA is associated with reduced breast cancer morbidity and a more favourable prognosis.

An Incremental EM-based Learning Approach for On-Line Prediction of Hospital Resource Utilization

Objective: Inpatient length of stay (LOS) is an important measure of hospital activity, health care resource consumption, and patient acuity. This research work aims at developing an incremental expectation maximization (EM) based learning approach on mixture of experts (ME) system for on-line prediction of LOS. The use of a batchmode learning process in most existing artificial neural networks to predict LOS is unrealistic, as the data become available over time and their pattern change dynamically. In contrast, an on-line process is capable of providing an output whenever a new datum becomes available. This on-the-spot information is therefore more useful and practical for making decisions, especially when one deals with a tremendous amount of data. Methods and material: The proposed approach is illustrated using a real example of gastroenteritis LOS data. The data set was extracted from a retrospective cohort study on all infants born in 1995-1997 and their subsequent admissions for gastroenteritis. The total number of admissions in this data set was n = 692. Linked hospitalization records of the cohort were retrieved retrospectively to derive the outcome measure, patient demographics, and associated co-morbidities information. A comparative study of the incremental learning and the batch-mode learning algorithms is considered. The performances of the learning algorithms are compared based on the mean absolute difference (MAD) between the predictions and the actual LOS, and the proportion of predictions with MAD < 1 day (Prop(MAD < 1)). The significance of the comparison is assessed through a regression analysis. Results: The incremental learning algorithm provides better on-line prediction of LOS when the system has gained sufficient training from more examples (MAD = 1.77 days and Prop(MAD < 1) = 54.3%), compared to that using the batch-mode learning. The regression analysis indicates a significant decrease of MAD (p-value = 0.063) and a significant (p-value = 0.044) increase of Prop(MAD

An analysis of the images attached to referral messages in an e-mail based telemedicine system for developing countries

Little is known about the quality of the images transmitted in email telemedicine systems. The present study was designed to survey the quality of images transmitted in the Swinfen Charitable Trust email referral system. Telemedicine cases were examined for a 3 month period in 2002 and a 3 month period in 2006. The number of cases with images attached increased from 8 (38%) to 37 (53%). There were four types of images (clinical photographs, microscope pictures, notes and X-ray images) and the proportion of radiology images increased from 27 to 48%. The cases in 2002 came from four different hospitals and were associated with seven different clinical specialties. In 2006, the cases came from 19 different hospitals and 20 different specialties. The 46 cases (from both study periods) had a total of 159 attached images. The quality of the images was assessed by awarding each image a score in four categories: focus, anatomical perspective, composition and lighting. The images were scored on a five-point scale (1 = very poor to 5 =very good) by a qualified medical photographer. In comparing image quality between the two study periods, there was some evidence that the quality had reduced, although the average size of the attached images had increased. The median score for all images in 2002 was 16 (interquartile range 14-19) and the median score in 2006 was 15 (13-16). The difference was significant (P < 0.001, Mann-Whitney test).

Simulation of air traffic controllers' behaviour using the operator choice model

The Operator Choice Model (OCM) was developed to model the behaviour of operators attending to complex tasks involving interdependent concurrent activities, such as in Air Traffic Control (ATC). The purpose of the OCM is to provide a flexible framework for modelling and simulation that can be used for quantitative analyses in human reliability assessment, comparison between human computer interaction (HCI) designs, and analysis of operator workload. The OCM virtual operator is essentially a cycle of four processes: Scan Classify Decide Action Perform Action. Once a cycle is complete, the operator will return to the Scan process. It is also possible to truncate a cycle and return to Scan after each of the processes. These processes are described using Continuous Time Probabilistic Automata (CTPA). The details of the probability and timing models are specific to the domain of application, and need to be specified using domain experts. We are building an application of the OCM for use in ATC. In order to develop a realistic model we are calibrating the probability and timing models that comprise each process using experimental data from a series of experiments conducted with student subjects. These experiments have identified the factors that influence perception and decision making in simplified conflict detection and resolution tasks. This paper presents an application of the OCM approach to a simple ATC conflict detection experiment. The aim is to calibrate the OCM so that its behaviour resembles that of the experimental subjects when it is challenged with the same task. Its behaviour should also interpolate when challenged with scenarios similar to those used to calibrate it. The approach illustrated here uses logistic regression to model the classifications made by the subjects. This model is fitted to the calibration data, and provides an extrapolation to classifications in scenarios outside of the calibration data. A simple strategy is used to calibrate the timing component of the model, and the results for reaction times are compared between the OCM and the student subjects. While this approach to timing does not capture the full complexity of the reaction time distribution seen in the data from the student subjects, the mean and the tail of the distributions are similar.

Sound source localisation through active audition

This paper presents a novel method for enabling a robot to determine the direction to a sound source through interacting with its environment. The method uses a new neural network, the Parameter-Less Self-Organizing Map algorithm, and reinforcement learning to achieve rapid and accurate response.

Expertise and the perception of kinematic and situational probability information

Two experiments were conducted on the nature of expert perception in the sport of squash. In the first experiment, ten expert and fifteen novice players attempted to predict the direction and force of squash strokes from either a film display (occluded at variable time periods before and after the opposing player had struck the ball) or a matched point-light display (containing only the basic kinematic features of the opponent's movement pattern). Experts outperformed the novices under both display conditions, and the same basic time windows that characterised expert and novice pick-up of information in the film task also persisted in the point-light task. This suggests that the experts' perceptual advantage is directly related to their superior pick-up of essential kinematic information. In the second experiment, the vision of six expert and six less skilled players was occluded by remotely triggered liquid-crystal spectacles at quasi-random intervals during simulated match play. Players were required to complete their current stroke even when the display was occluded and their prediction performance was assessed with respect to whether they moved to the correct half of the court to match the direction and depth of the opponent's stroke. Consistent with experiment 1, experts were found to be superior in their advance pick-up of both directional and depth information when the display was occluded during the opponent's hitting action. However, experts also remained better than chance, and clearly superior to less skilled players, in their prediction performance under conditions where occlusion occurred before any significant pre-contact preparatory movement by the opposing player was visible. This additional source of expert superiority is attributable to their superior attunement to the information contained in the situational probabilities and sequential dependences within their opponent's pattern of play.

The proportion of individual lucerne plants resistant to Phytophthora medicaginis and Colletotrichum trifolii in Australian cultivars

Phytophthora root rot (Phytophthora medicaginis) and colletotrichum crown rot (Colletotrichum trifoli) are the 2 most serious pathogens of lucerne in eastern Australia. Work reported in this paper shows that in glasshouse tests of the 11 most commonly grown Australian lucerne cultivars, the proportion of individual plants with resistance to both pathogens ranges from 0 (Hunter River and Aurora) through to a maximum of 19.8% (Sequel HR). Within 9 of the cultivars, the proportion of individual plants resistant to the 2 pathogens was <7%. Since these 2 diseases are known to cause serious losses in eastern Australia, the results indicate further improvement in lucerne production can be obtained by increasing the proportion of individual plants in a cultivar resistant to both pathogens. This would be best achieved by identifying dominant sources of resistance and incorporating this into on-going lucerne breeding programs.

Influence of patient age and implantation technique on the probability of re-replacement of the homograft aortic valve

Background and aim of the study: Results of valve re-replacement (reoperation) in 898 patients undergoing aortic valve replacement with cryopreserved homograft valves between 1975 and 1998 are reported. The study aim was to provide estimates of unconditional probability of valve reoperation and cumulative incidence function (actual risk) of reoperation. Methods: Valves were implanted by subcoronary insertion (n = 500), inclusion cylinder (n = 46), and aortic root replacement (n = 352). Probability of reoperation was estimated by adopting a mixture model framework within which estimates were adjusted for two risk factors: patient age at initial replacement, and implantation technique. Results: For a patient aged 50 years, the probability of reoperation in his/her lifetime was estimated as 44% and 56% for non-root and root replacement techniques, respectively. For a patient aged 70 years, estimated probability of reoperation was 16% and 25%, respectively. Given that a reoperation is required, patients with non-root replacement have a higher hazard rate than those with root replacement (hazards ratio = 1.4), indicating that non-root replacement patients tend to undergo reoperation earlier before death than root replacement patients. Conclusion: Younger patient age and root versus non-root replacement are risk factors for reoperation. Valve durability is much less in younger patients, while root replacement patients appear more likely to live longer and hence are more likely to require reoperation.

Derivation of the probability distribution function for the local density of states of a disordered quantum wire via the replica trick and supersymmetry

We consider the statistical properties of the local density of states of a one-dimensional Dirac equation in the presence of various types of disorder with Gaussian white-noise distribution. It is shown how either the replica trick or supersymmetry can be used to calculate exactly all the moments of the local density of states.' Careful attention is paid to how the results change if the local density of states is averaged over atomic length scales. For both the replica trick and supersymmetry the problem is reduced to finding the ground state of a zero-dimensional Hamiltonian which is written solely in terms of a pair of coupled spins which are elements of u(1, 1). This ground state is explicitly found for the particular case of the Dirac equation corresponding to an infinite metallic quantum wire with a single conduction channel. The calculated moments of the local density of states agree with those found previously by Al'tshuler and Prigodin [Sov. Phys. JETP 68 (1989) 198] using a technique based on recursion relations for Feynman diagrams. (C) 2001 Elsevier Science B.V. All rights reserved.

Cavity QED analog of the harmonic-oscillator probability distribution function and quantum collapses

We establish a connection between the simple harmonic oscillator and a two-level atom interacting with resonant, quantized cavity and strong driving fields, which suggests an experiment to measure the harmonic-oscillator's probability distribution function. To achieve this, we calculate the Autler-Townes spectrum by coupling the system to a third level. We find that there are two different regions of the atomic dynamics depending on the ratio of the: Rabi frequency Omega (c) of the cavity field to that of the Rabi frequency Omega of the driving field. For Omega (c)

Laying date and laying sequence influence the sex ratio of crimson rosella broods

We examine the patterns of sex allocation in crimson rosellas Platycercus elegans, a socially monogamous Australian parrot. Overall, 41.8% of nestlings were male, a significant female bias. However underlying this population-level bias were non-random patterns of sex allocation within broods. Broods produced early in the season were female-biased, but the proportion of males in a brood increased as the breeding season progressed. Female rosellas may obtain greater fitness benefits from early-fledging daughters than sons because daughters can breed as 1-year-olds whereas sons do not breed until they are at least 2 years old. Laying date and laying sequence also interacted to influence the sex ratio of eggs. The sex of early-laid eggs strongly followed the brood level pattern, whereas the sex of middle- and late-laid eggs did not change significantly as the season progressed. Nevertheless, late-laid eggs were very unlikely to be male at the end of the season. We argue these differing seasonal patterns reflect the relative costs and benefits to producing early-hatched males and females at different times of the season. Female rosellas appear to maximise the probability that daughters are able to breed early but to minimise competitive asymmetries within the brood. In particular, late-hatched male chicks are disadvantaged if their oldest sibling is male, explaining the dearth of broods containing late-hatched males at the end of the breeding season.