Maximum Likelihood Estimation of Mixture Densities for Binned and Truncated Multivariate Data

Binning and truncation of data are common in data analysis and machine learning. This paper addresses the problem of fitting mixture densities to multivariate binned and truncated data. The EM approach proposed by McLachlan and Jones (Biometrics, 44: 2, 571-578, 1988) for the univariate case is generalized to multivariate measurements. The multivariate solution requires the evaluation of multidimensional integrals over each bin at each iteration of the EM procedure. Naive implementation of the procedure can lead to computationally inefficient results. To reduce the computational cost a number of straightforward numerical techniques are proposed. Results on simulated data indicate that the proposed methods can achieve significant computational gains with no loss in the accuracy of the final parameter estimates. Furthermore, experimental results suggest that with a sufficient number of bins and data points it is possible to estimate the true underlying density almost as well as if the data were not binned. The paper concludes with a brief description of an application of this approach to diagnosis of iron deficiency anemia, in the context of binned and truncated bivariate measurements of volume and hemoglobin concentration from an individual's red blood cells.

Treatment outcome in individuals with chronic pain: is the Pain Stages of Change Questionnaire (PSOCQ) a useful tool? Jenny StrongCorresponding Author Contact Information, E-mail The Corresponding Author, a, Kym Westburya, Glen Smithb, Ian McKenziec and William Ryan

The efficacy of psychological treatments emphasising a self-management approach to chronic pain has been demonstrated by substantial empirical research. Nevertheless, high drop-out and relapse rates and low or unsuccessful engagement in self-management pain rehabilitation programs have prompted the suggestion that people vary in their readiness to adopt a self-management approach to their pain. The Pain Stages of Change Questionnaire (PSOCQ) was developed to assess a patient's readiness to adopt a self-management approach to their chronic pain. Preliminary evidence has supported the PSOCQ's psychometric properties. The current study was designed to further examine the psychometric properties of the PSOCQ, including its reliability, factorial structure and predictive validity. A total of 107 patients with an average age of 36.2 years (SD = 10.63) attending a multi-disciplinary pain management program completed the PSOCQ, the Pain Self-Efficacy Questionnaire (PSEQ) and the West Haven-Yale Multidimensional Pain Inventory (WHYMPI) pre-admission and at discharge from the program. Initial data analysis found inadequate internal consistencies of the precontemplation and action scales of the PSOCQ and a high correlation (r = 0.66, P < 0.01) between the action and maintenance scales. Principal component analysis supported a two-factor structure: 'Contemplation' and 'Engagement'. Subsequent analyses revealed that the PSEQ was a better predictor of treatment outcome than the PSOCQ scales. Discussion centres upon the utility of the PSOCQ in a clinical pain setting in light of the above findings, and a need for further research. (C) 2002 International Association for the Study of Pain. Published by Elsevier Science B.V. All rights reserved.

Quantitative evaluation of altered hepatic spaces and membrane transport in fibrotic rat liver

Four animal models were used to quantitatively evaluate hepatic alterations in this study: (1) a carbon tetrachloride control group (phenobarbital treatment only), (2) a CCl4-treated group (phenobarbital with CCl4 treatment), (3) an alcohol-treated group (liquid diet with alcohol treatment), and (4) a pair-fed alcohol control group (liquid diet only). At the end of induction, single-pass perfused livers were used to conduct multiple indicator dilution (MID) studies. Hepatic spaces (vascular space, extravascular albumin space, extravascular sucrose space, and cellular distribution volume) and water hepatocyte permeability/surface area product were estimated from nonlinear regression of outflow concentration versus time profile data. The hepatic extraction ratio of H-3-taurocholate was determined by the nonparametric moments method. Livers were then dissected for histopathologic analyses (e.g., fibrosis index, number of fenestrae). In these 4 models, CCl4-treated rats were found to have the smallest vascular space, extravascular albumin space, H-3-taurocholate extraction, and water hepatocyte permeability/surface area product but the largest extravascular sucrose space and cellular distribution volume. In addition, a linear relationship was found to exist between histopathologic analyses (fibrosis index or number of fenestrae) and hepatic spaces. The hepatic extraction ratio of H-3-taurocholate and water hepatocyte permeability/surface area product also correlated to the severity of fibrosis as defined by the fibrosis index. In conclusion, the multiple indicator dilution data obtained from the in situ perfused rat liver can be directly related to histopathologic analyses.

Rational combinatorial chemistry-based selection, synthesis and evaluation of an affinity adsorbent for recombinant human clotting factor VII

The selection, synthesis and chromatographic evaluation of a synthetic affinity adsorbent for human recombinant factor VIIa is described. The requirement for a metal ion-dependent immunoadsorbent step in the purification of the recombinant human clotting factor, FVIIa, has been obviated by using the X-ray crystallographic structure of the complex of tissue factor (TF) and Factor VIIa and has directed our combinatorial approach to select, synthesise and evaluate a rationally-selected affinity adsorbent from a limited library of putative ligands. The selected and optimised ligand comprises a triazine scaffold bis-substituted with 3-aminobenzoic acid and has been shown to bind selectively to FVIIa in a Ca2+-dependent manner. The adsorbent purifies FVIIa to almost identical purity (>99%), yield (99%), activation/degradation profile and impurity content (∼1000 ppm) as the current immunoadsorption process, while displaying a 10-fold higher static capacity and substantially higher reusability and durability. © 2002 Elsevier Science B.V. All rights reserved.

Air-water flows down stepped chutes: turbulence and flow structure observations

Interactions between turbulent waters and atmosphere may lead to strong air-water mixing. This experimental study is focused on the flow down a staircase channel characterised by very strong flow aeration and turbulence. Interfacial aeration is characterised by strong air-water mixing extending down to the invert. The size of entrained bubbles and droplets extends over several orders of magnitude, and a significant number of bubble/droplet clusters was observed. Velocity and turbulence intensity measurements suggest high levels of turbulence across the entire air-water flow. The increase in turbulence levels, compared to single-phase flow situations, is proportional to the number of entrained particles. (C) 2002 Elsevier Science Ltd. All rights reserved.

Clinical evaluation of the WOMAC 3.0 OA Index in numeric rating scale format using a computerized touch screen version

Background: The Western Ontario and McMaster Universities (WOMAC) Osteoarthritis Index is a previously described self-administered questionnaire covering three domains: pain, stiffness and function. It has been validated in patients with osteoarthritis (OA) of the hip or knee in a paper-based format. Aim: To validate the WOMAC 3.0 using a numerical rating scale in a computerized touch screen format allowing immediate evaluation of the questionnaire. In the computed version cartoons, written and audio instruments were included in order facilitate application. Methods: Fifty patients, demographically balanced, with radiographically proven primary hip or knee OA completed the classical paper and the new computerized WOMAC version. Subjects were randomized either to paper format or computerized format first to balance possible order effects, Results: The intra-class correlation coefficients for pain, stiffness and function values were 0.915, 0.745 and 0.940, respectively. The Spearman correlation coefficients for pain, stiffness and function were 0.88, 0.77 and 0.87, respectively. Conclusion: These data indicate that the computerized WOMAC OA index 3.0 is comparable to the paper WOMAC in all three dimensions. The computerized version would allow physicians to get an immediate result and if present a direct comparison with a previous exam. (C) 2002 OsteoArthritis Research Society International. Published by Elsevier Science Ltd. All rights reserved.

Social structure in migrating humpback whales (Megaptera novaeangliae)

Although largely solitary, humpback whales exhibit a number of behaviours where individuals co-operate with one another, for example during bubble net feeding. Such cases could be due to reciprocal altruism brought on by exceptional circumstances, for example the presence of abundant shoaling fish. An alternative explanation is that these behaviours have evolved through kin selection. With little restriction to either communication or movement, diffuse groups of relatives could maintain some form of social organization without the need to travel in tight-nit units. To try to distinguish between these hypotheses, we took advantage of the fact that migrating humpback whales often swim together in small groups. If kin selection is important in humpback whale biology, these groups should be enriched for relatives. Consequently, we analysed biopsy samples from 57 groups of humpback whales migrating off Eastern Australia in 1992. A total of 142 whales were screened for eight microsatellite markers. Mitochondrial DNA sequences (371 bp) were also used to verify and assist kinship identification. Our data add support to the notion that mothers travel with their offspring for the first year of the calf's life. However, beyond the presence of mother-calf/yearling pairs, no obvious relatedness pattern was found among whales sampled either in the same pod or on the same day. Levels of relatedness did not vary between migratory phases (towards or away from the breeding ground), nor between the two sexes considered either overall or in the north or south migrations separately. These findings suggest that, if any social organization does exist, it is formed transiently when needed rather than being a constant feature of the population, and hence is more likely based on reciprocal altruism than kin selection.

Evaluation of various pre-treatments for the dehydration of banana and selection of suitable drying models

The thin-layer drying behaviour of bananas in a beat pump dehumidifier dryer was examined. Four pre-treatments (blanching, chilling, freezing and combined blanching and freezing) were applied to the bananas, which were dried at 50 degreesC with an air velocity of 3.1 m s(-1) and with the relative humidity of the inlet air of 10-35%. Three drying models, the simple model, the two-term exponential model and the Page model were examined. All models were evaluated using three statistical measures, correlation coefficient, root means square error, and mean absolute percent error. Moisture diffusivity was calculated based on the diffusion equation for an infinite cylindrical shape using the slope method. The rate of drying was higher for the pre-treatments involving freezing. The sample which was blanched only did not show any improvement in drying rate. In fact, a longer drying time resulted due to water absorption during blanching. There was no change in the rate for the chilled sample compared with the control. While all models closely fitted the drying data, the simple model showed greatest deviation from the experimental results. The two-term exponential model was found to be the best model for describing the drying curves of bananas because its parameters represent better the physical characteristics of the drying process. Moisture diffusivities of bananas were in the range 4.3-13.2 x 10(-10) m(2)s(-1). (C) 2002 Published by Elsevier Science Ltd.

Preparing students for the virtual organisation: an evaluation of learning with virtual learning technologies

The paper explores the development of learning behaviours in a virtual management course and the factors that impacted on this development. Data suggest that most teams experienced three kinds of learning behaviours – social, operational and content learning. We propose that the need for technical expertise and team participation will vary during these different stages of learning. Addressing the characteristics of these stages, we comment on the development of a ‘completion phase’ of team development. We argue that the extent to which teams demonstrate different learning stages has a significant impact on the development of on-line learning behaviours. Discussing these results, we suggest why different teams develop distinct learning behaviours, with accordant emphasis on teaching as a moderating and co ordinating role, despite current virtual team pedagogical expectations.

The Three-dimensional Solution Structure of NaD1, a New Floral Defensin from Nicotiana alata and its Application to a Homology Model of the Crop Defense Protein alfAFP

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of NaD1, a novel antifungal and insecticidal protein isolated from the flowers of Nicotiana alata. NaD1 is a basic, cysteine-rich protein of 47 residues and is the first example of a plant defensin from flowers to be characterized structurally. Its three-dimensional structure consists of an a-helix and a triple-stranded anti-parallel beta-sheet that are stabilized by four intramolecular disulfide bonds. NaD1 features all the characteristics of the cysteine-stabilized up motif that has been described for a variety of proteins of differing functions ranging from antibacterial insect defensins and ion channel-perturbing scorpion toxins to an elicitor of the sweet taste response. The protein is biologically active against insect pests, which makes it a potential candidate for use in crop protection. NaD1 shares 31% sequence identity with alfAFP, an antifungal protein from alfalfa that confers resistance to a fungal pathogen in transgenic potatoes. The structure of NaD1 was used to obtain a homology model of alfAFP, since NaD1 has the highest level of sequence identity with alfAFP of any structurally characterized antifungal defensin. The structures of NaD1 and alfAFP were used in conjunction with structure - activity data for the radish defensin Rs-AFP2 to provide an insight into structure-function relationships. In particular, a putative effector site was identified in the structure of NaD1 and in the corresponding homology model of alfAFP. (C) 2002 Elsevier Science Ltd. All rights reserved.

Secondary structure of the third extracellular loop responsible for ligand selectivity of a mammalian gonadotropin-releasing hormone receptor

The extracellular loop 3 (ECL3) of the mammalian gonadotropin-releasing hormone receptor (GnRH-R) contains an acidic amino acid (Glu(301) in the mouse GnRH-R,) that confers agonist selectivity for Are in mammalian GnRH. It is proposed that a specific conformation of ECL3 is necessary to orientate the carboxyl side chain of the acidic residue for interaction with Arg(8) of GnRH, which is supported by decreased affinity for Arg(8) GnRH but not Gln(8) GnRH when an adjacent Pro is mutated to Ala. To probe the structural contribution of the loop domain to the proposed presentation of the carboxyl side chain, we synthesized a model peptide (CGPEMLNRVSEPGC) representing residues 293-302 of mouse ECL3, where Cys and Gly residues are added symmetrically at the N and C termini, respectively, allowing the introduction of a disulfide bridge to simulate the distances at which the ECL3 is tethered to the transmembrane domains 6 and 7 of the receptor. The ability of the ECL3 peptide to bind GnRH with low affinity was demonstrated by its inhibition of GnRH stimulation of inositol phosphate production in cells expressing the GnRH-R. The CD bands of the ECL3 peptides exhibited a superposition of predominantly unordered structure and partial contributions from beta-sheet structure. Likewise, the analysis of the amide I and amide III bands from micro-Raman and FT Raman experiments revealed mainly unordered conformations of the cyclic and of the linear peptide. NMR data demonstrated the presence of a beta-hairpin among an ensemble of largely disordered structures in the cyclic peptide. The location of the turn linking the two strands of the hairpin was assigned to the three central residues L-296, N-297, and R-298. A small population of structured species among an ensemble of predominantly random coil conformation suggests that the unliganded receptor represents a variety of structural conformers, some of which have the potential to make contacts with the ligand. We propose a mechanism of receptor activation whereby binding of the agonist to the inactive receptor state induces and stabilizes a particular structural state of the loop domain, leading to further conformational rearrangements across the transmembrane domain and signal propagating interaction with G proteins. Interaction of the Glu(301) of the receptor with Arg(8) of GnRH induces a folded configuration of the ligand. Our proposal thus suggests that conformational changes of both ligand and receptor result from this interaction.

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.

Systematic characterization of mutations in yeast acetohydroxyacid synthase: Interpretation of herbicide-resistance data

Acetohydroxyacid synthase (AHAS, EC 4.1.3.18) catalyses the first step in branched-chain amino acid biosynthesis and is the target for sulfonylurea and imidazolinone herbicides, which act as potent and specific inhibitors. Mutants of the enzyme have been identified that are resistant to particular herbicides. However, the selectivity of these mutants towards various sulfonylureas and imidazolinones has not been determined systematically. Now that the structure of the yeast enzyme is known, both in the absence and presence of a bound herbicide, a detailed understanding of the molecular interactions between the enzyme and its inhibitors becomes possible. Here we construct 10 active mutants of yeast AHAS, purify the enzymes and determine their sensitivity to six sulfonylureas and three imidazolinones. An additional three active mutants were constructed with a view to increasing imidazolinone sensitivity. These three variants were purified and tested for their sensitivity to the imidazolinones only. Substantial differences are observed in the sensitivity of the 13 mutants to the various inhibitors and these differences are interpreted in terms of the structure of the herbicide-binding site on the enzyme.

Numerical simulations of wind and temperature structure within an Alpine lake basin, Lake Tekapo, New Zealand

High-resolution numerical model simulations have been used to study the local and mesoscale thermal circulations in an Alpine lake basin. The lake (87 km(2)) is situated in the Southern Alps, New Zealand and is located in a glacially excavated rock basin surrounded by mountain ranges that reach 3000 m in height. The mesoscale model used (RAMS) is a three-dimensional non-hydrostatic model with a level 2.5 turbulence closure scheme. The model demonstrates that thermal forcing at local (within the basin) and regional (coast-to-basin inflow) scales drive the observed boundary-layer airflow in the lake basin during clear anticyclonic summertime conditions. The results show that the lake can modify (perturb) both the local and regional wind systems. Following sunrise, local thermal circulations dominate, including a lake breeze component that becomes embedded within the background valley wind system. This results in a more divergent flow in the basin extending across the lake shoreline. However, a closed lake breeze circulation is neither observed nor modelled. Modelling results indicate that in the latter part of the day when the mesoscale (coast-to-basin) inflow occurs, the relatively cold pool of lake air in the basin can cause the intrusion to decouple from the surface. Measured data provide qualitative and quantitative support for the model results.

Evaluation of cyclooxygenase 2 inhibitor use in patients admitted to a large teaching hospital

Background: The heavy usage of coxibs in Australia far outstrips the predicted usage that was based on the treatment of patients with risk factors for upper gastro-intestinal adverse events from conventional anti--inflammatory agents. This raises questions regarding the appropriateness of prescribing. Aims: To determine: (i) the relationship between prescriptions for cyclooxygenase 2 (COX-2) inhibitors and objective evidence of inflammatory arthritis, (ii) prior experience with paracetamol and/or conventional non-steroidal anti-inflammatory drugs (NSAIDs), and (iii) contraindications to the use of NSAIDs. Methods: Drug utilization evaluation and rheumato-logical assessment was conducted on 70 consecutive patients admitted on COX-2 inhibitors to a 480-bed metropolitan hospital. The main outcome measures were: the indication for COX-2 inhibitor; objective -evidence of inflammatory arthritis; previous trial of -paracetamol or conventional NSAIDs; and patient -satisfaction. Results: Only 11 patients (16%) had symptoms or signs of an inflammatory arthropathy, and met Pharmaceut-ical Benefits Schedule criteria for prescribing a COX-2 inhibitor. Fifty-nine patients (84%) had chronic osteo-arthritis, degenerative spinal disease, injury or malignancy, without overt active inflammation. Fourteen patients (20%) had trialled regular paracetamol prior to using any NSAID treatment. Conventional NSAIDs had been previously used by 51 patients (73%). Eleven patients (16%) reported previous adverse gastrointestinal effects from conventional NSAIDs. On the basis of significant renal impairment (creatinine clearance 5/10). Conclusions: Drug utilization data indicate that COX-2 inhibitors are frequently used first line for degenerative osteoarthritis in the absence of overt inflammation, without prior adequate trial of paracetamol and with disregard for the cautions and contraindications of these agents. These findings may explain the unprecedented Pharmaceutical Benefits Schedule expenditure on COX-2 inhibitors in Australia.