Modelling High-Dimensional Data by Mixtures of Factor Analyzers

We focus on mixtures of factor analyzers from the perspective of a method for model-based density estimation from high-dimensional data, and hence for the clustering of such data. This approach enables a normal mixture model to be fitted to a sample of n data points of dimension p, where p is large relative to n. The number of free parameters is controlled through the dimension of the latent factor space. By working in this reduced space, it allows a model for each component-covariance matrix with complexity lying between that of the isotropic and full covariance structure models. We shall illustrate the use of mixtures of factor analyzers in a practical example that considers the clustering of cell lines on the basis of gene expressions from microarray experiments. (C) 2002 Elsevier Science B.V. All rights reserved.

Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

A note on quasi-stationary distributions of birth-death processes and the SIS logistic epidemic

For Markov processes on the positive integers with the origin as an absorbing state, Ferrari, Kesten, Martinez and Picco studied the existence of quasi-stationary and limiting conditional distributions by characterizing quasi-stationary distributions as fixed points of a transformation Phi on the space of probability distributions on {1, 2,.. }. In the case of a birth-death process, the components of Phi(nu) can be written down explicitly for any given distribution nu. Using this explicit representation, we will show that Phi preserves likelihood ratio ordering between distributions. A conjecture of Kryscio and Lefevre concerning the quasi-stationary distribution of the SIS logistic epidemic follows as a corollary.

Conceptual modelling and ontology: Possibilities and pitfalls

Conceptual modelling is an activity undertaken during information systems development work to build a representation of selected semantics about some real-world domain. Ontological theories have been developed to account for the structure and behavior of the real world in general. In this paper, I discuss why ontological theories can be used to inform conceptual modelling research, practice, and pedagogy. I provide examples from my research to illustrate how a particular ontological theory has enabled me to improve my understanding of certain conceptual modelling practices and grammars. I describe, also, how some colleagues and I have used this theory to generate several counter-intuitive, sometimes surprising predictions about widely advocated conceptual modelling practices - predictions that subsequently were supported in empirical research we undertook. Finally, I discuss several possibilities and pitfalls I perceived to be associated with our using ontological theories to underpin research on conceptual modelling.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

Conspicuous males suffer higher predation risk: visual modelling and experimental evidence from lizards

Colour pattern variation is a striking and widespread phenomenon. Differential predation risk between individuals is often invoked to explain colour variation, but empirical support for this hypothesis is equivocal. We investigated differential conspicuousness and predation risk in two species of Australian rock dragons, Ctenophorus decresii and C. vadnappa. To humans, the coloration of males of these species varies between 'bright' and 'dull'. Visual modelling based on objective colour measurements and the spectral sensitivities of avian visual pigments showed that dragon colour variants are differentially conspicuous to the visual system of avian predators when viewed against the natural background. We conducted field experiments to test for differential predation risk, using plaster models of 'bright' and 'dull' males. 'Bright' models were attacked significantly more often than 'dull' models suggesting that differential conspicuousness translates to differential predation risk in the wild. We also examined the influence of natural geographical range on predation risk. Results from 22 localities suggest that predation rates vary according to whether predators are familiar with the prey species. This study is among the first to demonstrate both differential conspicuousness and differential predation risk in the wild using an experimental protocol. (C) 2003 Published by Elsevier Ltd on behalf of The Association for the Study of Animal Behaviour.

Discrete element modelling of the influence of lifters on power draw of tumbling mills

Crushing and grinding are the most energy intensive part of the mineral recovery process. A major part of rock size reduction occurs in tumbling mills. Empirical models for the power draw of tumbling mills do not consider the effect of lifters. Discrete element modelling was used to investigate the effect of lifter condition on the power draw of tumbling mill. Results obtained with PFC3D code show that lifter condition will have a significant influence on the power draw and on the mode of energy consumption in the mill. Relatively high lifters will consume less power than low lifters, under otherwise identical conditions. The fraction of the power that will be consumed as friction will increase as the height of the lifters decreases. This will result in less power being used for high intensity comminution caused by the impacts. The fraction of the power that will be used to overcome frictional resistance is determined by the material's coefficient of friction. Based on the modelled results, it appears that the effective coefficient of friction for in situ mill is close to 0.1. (C) 2003 Elsevier Science Ltd. All rights reserved.

Applying discrete element modelling to vertical and horizontal shaft impact crushers

The PFC3D (particle flow code) that models the movement and interaction of particles by the DEM techniques was employed to simulate the particle movement and to calculate the velocity and energy distribution of collision in two types of impact crusher: the Canica vertical shaft crusher and the BJD horizontal shaft swing hammer mill. The distribution of collision energies was then converted into a product size distribution for a particular ore type using JKMRC impact breakage test data. Experimental data of the Canica VSI crusher treating quarry and the BJD hammer mill treating coal were used to verify the DEM simulation results. Upon the DEM procedures being validated, a detailed simulation study was conducted to investigate the effects of the machine design and operational conditions on velocity and energy distributions of collision inside the milling chamber and on the particle breakage behaviour. (C) 2003 Elsevier Ltd. All rights reserved.

Modelling the size of the crushed zone around a blasthole

A new model to predict the extent of crushing around a blasthole is presented. The model is based on the back-analysis of a comprehensive experimental program that included the direct measurement of the zone of crushing from 92 blasting tests on concrete blocks using two commercial explosives. The concrete blocks varied from low, medium to high strength and measured 1.5 in in length, 1.0 m in width and 1.1 m in height. A dimensionless parameter called the crushing zone index (CZI) is introduced. This index measures the crushing potential of a charged blasthole and is a function of the borehole pressure, the unconfined compressive strength of the rock material, dynamic Young's modulus and Poisson's ratio. It is shown that the radius of crushing is a function of the CZI and the blasthole radius. A good correlation between the new model and measured results was obtained. A number of previously proposed models could not approximate the conditions measured in the experimental work and there are noted discrepancies between the different approaches reviewed, particularly for smaller diameter holes and low strength rock conditions. The new model has been verified with full scale tests reported in the literature. Results from this validation and model evaluations show its applicability to production blasting. (C) 2003 Elsevier Science Ltd. All rights reserved.

Discrete element modelling of power draw of tumbling mills

The power required to operate large mills is typically 5-10 MW. Hence, optimisation of power consumption will have a significant impact on overall economic performance and environmental impact. Power draw modelling results using the discrete element code PFC3D have been compared with results derived from the widely used empirical Model of Morrell. This is achieved by calculating the power draw for a range of operating conditions for constant mill size and fill factor using two modelling approaches. fThe discrete element modelling results show that, apart from density, selection of the appropriate material damping ratio is critical for the accuracy of modelling of the mill power draw. The relative insensitivity of the power draw to the material stiffness allows selection of moderate stiffness values, which result in acceptable computation time. The results obtained confirm that modelling of the power draw for a vertical slice of the mill, of thickness 20% of the mill length, is a reliable substitute for modelling the full mill. The power draw predictions from PFC3D show good agreement with those obtained using the empirical model. Due to its inherent flexibility, power draw modelling using PFC3D appears to be a viable and attractive alternative to empirical models where necessary code and computer power are available.

Polymerase chain reaction tests for the identification of Ross River, Kunjin and Murray Valley encephalitis virus infections in horses

Objective To develop and validate specific, sensitive and rapid diagnostic tests using RT-PCR for the detection of Ross River virus (RRV), Kunjin virus (KV) and Murray Valley encephalitis virus (MVEV) infections in horses. Methods Primer sets based on nucleotide sequence encoding the envelope glycoprotein E2 of RRV and on the nonstructural protein 5 (NS5) of KV and MVEV were designed and used in single round PCRs to test for the respective viruses in infected cell cultures and, in the case of RRV, in samples of horse blood and synovial fluid. Results The primer pairs designed for each of the three viruses amplified a product of expected size from prototype viruses that were grown in cell culture. The identity of each of the products was confirmed by nucleotide sequencing indicating that in the context used the RT-PCRs were specific. RRV was detected in serums from 8 horses for which there were clinical signs consistent with RRV infection such that an acute-phase serum sample was taken and submitted for RRV serology testing. The RRV RT-PCR was analytically sensitive in that it was estimated to detect as little as 50 TCID50 of RRV per mL of serum and was specific in that the primer pairs did not amplify other products from the 8 serum samples. The RRV primers also detected virus in three independent mosquito pools known to contain RRV by virus isolation in cell culture. Samples from horses suspected to be infected with KV and MVEV were not available. Conclusion Despite much anecdotal and serological evidence for infection of horses with RRV actual infection and associated clinical disease are infrequently confirmed. The availability of a specific and analytically sensitive RT-PCR for the detection of RRV provides additional opportunities to confirm the presence of this virus in clinical samples. The RTPCR primers for the diagnosis of KV and MVEV infections were shown to be specific for cell culture grown viruses but the further validation of these tests requires the availability of appropriate clinical samples from infected horses.

Quadriceps activation in closed and in open kinetic chain exercise

Purpose: For treatment of various knee disorders, muscles are trained in open or closed kinetic chain tasks. Coordination between the heads of the quadriceps muscle is important for stability and optimal joint loading for both the tibiofemoral and the patellofemoral joint. The aim of this study was to examine whether the quadriceps femoris muscles are activated differently in open versus closed kinetic chain tasks. Methods: Ten healthy men and women (mean age 28.5 +/- 0.7) extended the knees isometrically in open and closed kinetic chain tasks in a reaction time paradigm using moderate force. Surface electromyography (EMG) recordings were made from four different parts of the quadriceps muscle. The onset and amplitude of EMG and force data were measured. Results: In closed chain knee extension, the onset of EMG activity of the four different muscle portions of the quadriceps was more simultaneous than in the open chain. In open chain, rectus femoris (RF) had the earliest EMG onset while vastus medialis obliquus was activated last (7 +/- 13 ms after RF EMG onset) and with smaller amplitude (40 +/- 30% of maximal voluntary contraction (MVC)) than in closed chain (46 +/- 43% MVC). Conclusions: Exercise in closed kinetic chain promotes more balanced initial quadriceps activation than does exercise in open kinetic chain. This may be of importance in designing training programs aimed toward control of the patellofemoral joint.

Alternatively spliced products of the human kinesin light chain 1 (KNS2) gene

Conventional kinesin is a microtubule-based molecular motor involved in the transport of membranous and non-membranous cargoes. The kinesin holoenzyme exists as a heterotetramer, consisting of two heavy chain and two light chain subunits. It is thought that one function of the light chains is to interact with the cargo. Alternative splicing of kinesin light chain pre-mRNA has been observed in lower organisms, although evidence for alternative splicing of the human gene has not been reported. We have identified 19 variants of the human KNS2 gene (KLC1) that are generated by alternative splicing of downstream exons, but calculate that KNS2 has the potential to produce 285919 spliceforms. Corresponding spliceforms of the mouse KLC1 gene were also identified. The alternative exons are all located 3' of exon 12 and the novel spliceforms produce both alternative carboxy termini and alternative 3' untranslated regions. The observation of multiple light chain isoforms is consistent with their proposed role in specific cargo attachment.