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.

Low-density koala (Phascolarctos cinereus) populations in the mulgalands of south-west Queensland. I. Faecal pellet sampling protocol

Koala (Phascolarctos cinereus) populations in eastern Australia are threatened by land clearing for agricultural and urban development. At the same time, conservation efforts are hindered by a dearth of information about inland populations. Faecal deposits offer a source of information that is readily available and easily collected non-invasively. We detail a faecal pellet sampling protocol that was developed for use in a large rangeland biogeographic region. The method samples trees in belt transects, uses a thorough search at the tree base to quickly identify trees with koala pellets under them, then estimates the abundance of faecal pellets under those trees using 1-m(2) quadrats. There was a strong linear relationship between these estimates and a complete enumeration of pellet abundance under the same trees. We evaluated the accuracy of our method in detecting trees where pellets were present by means of a misclassification index that was weighed more heavily for missed trees that had high numbers of pellets under them. This showed acceptable accuracy in all landforms except riverine, where some trees with large numbers of pellets were missed. Here, accuracy in detecting pellet presence was improved by sampling with quadrats, rather than basal searches. Finally, we developed a method to reliably age pellets and demonstrate how this protocol could be used with the faecal-standing-crop method to derive a regional estimate of absolute koala abundance.

Wolbachia density and virulence attenuation after transfer into a novel host

The factors that control replication rate of the intracellular bacterium Wolbachia pipientis in its insect hosts are unknown and difficult to explore, given the complex interaction of symbiont and host genotypes. Using a strain of Wolbachia that is known to over-replicate and shorten the lifespan of its Drosophila melanogaster host, we have tracked the evolution of replication control in both somatic and reproductive tissues in a novel host/Wolbachia association. After transinfection (the transfer of a Wolbachia strain into a different species) of the over-replicating Wolbachia popcorn strain from D. metanogaster to Drosophila simulans, we demonstrated that initial high densities in the ovaries were in excess of what was required for perfect maternal transmission, and were likely causing reductions in reproductive fitness. Both densities and fitness costs associated with ovary infection rapidly declined in the generations after transinfection. The early death effect in D. simulans attenuated only slightly and was comparable to that induced in D. metanogaster. This study reveals a strong host involvement in Wolbachia replication rates, the independence of density control responses in different tissues, and the strength of natural selection acting on reproductive fitness.

Numerical investigations of flow and energy fields near a thermoacoustic couple

The flow field and the energy transport near thermoacoustic couples are simulated using a 2D full Navier-Stokes solver. The thermoacoustic couple plate is maintained at a constant temperature; plate lengths, which are short and long compared with the particle displacement lengths of the acoustic standing waves, are tested. Also investigated are the effects of plate spacing and the amplitude of the standing wave. Results are examined in the form of energy vectors, particle paths, and overall entropy generation rates. These show that a net heat-pumping effect appears only near the edges of thermoacoustic couple plates, within about a particle displacement distance from the ends. A heat-pumping effect can be seen even on the shortest plates tested when the plate spacing exceeds the thermal penetration depth. It is observed that energy dissipation near the plate increases quadratically as the plate spacing is reduced. The results also indicate that there may be a larger scale vortical motion outside the plates which disappears as the plate spacing is reduced. (C) 2002 Acoustical Society of America.

The relationship of natural killer cell counts, perforin mRNA and CD2 expression to post-exercise natural killer cell activity in humans

The aim of this study was to further investigate the mechanism of suppression of natural killer (NK) cell cytotoxic activity In peripheral blood following strenuous exercise. Blood was collected for analysis of NK cell concentration, cytotoxic activity, CD2 surface expression and perforin gene expression from runners (RUN, n = 6) and resting controls (CONTROL, n = 4) pre-exercise, 0, 1.5, 5, and 24 h following a 60-min treadmill run at 80% of VO2 peak. Natural killer cytotoxic activity, measured using a whole blood chromium release assay, fluctuated minimally in the CONTROL group and increased by 63% and decreased by 43% 0 and 1.5 h post-exercise, respectively, in the RUN group (group x time, P < 0.001). Lytic index (cytotoxic activity per cell) did not change. Perforin mRNA, measured using quantitative real-time polymerase chain reaction (ORT-PCR) decreased from pre- to post-exercise and remained decreased through 24 h, The decrease from pre- to 0 In post-exercise was seen predominately in the RUN group and was inversely correlated r = - 0.95) to pre-exercise perform mRNA. The NK cell surface expression of CD2 (lymphocyte function-associated antigen-2) was determined using fluorescent antibodies and flow cytometry, There was no change in the proportion of NK cells expressing CD2 or CD2 density, We conclude that (1) numerical redistribution accounted for most of the change in NK cytotoxic activity following a strenuous run, (2) decrease in perforin gene expression during the run was inversely related to pre-exercise levels but did not parallel changes in cytotoxic activity, and (3) CD2 surface expression was not affected by exercise.

On the effect of tracer density on axial dispersion in a batch oscillatory baffled column

In this paper, we report our modelling evaluation on the effect of tracer density on axial dispersion in a batch oscillatory baffled column (OBC). Tracer solution of potassium nitrite, its specific density ranged from 1.0 to 1.5, was used in the study, and was injected to the vertical column from either the top or bottom. Local concentration profiles are measured using conductivity probes at two locations along the height of the column. Using the experimental measured concentration profiles together with both 'Tank-in-Series' and 'Plug Flow with Axial Dispersion' models, axial dispersion coefficients were determined and used to describe the effect of specific tracer density on mixing in the OBC. The results showed that the axial dispersion coefficients evaluated by the two models are very similar in both magnitudes and trends, and the range of variations in such coefficients is generally larger for the bottom injection than for the top one. Empirical correlations linking the mechanical energy for mixing, the specific density of tracer and axial dispersion coefficient were established. Using these correlations, we identified the enhancements of up to 269% on axial dispersion for various specific tracer densities. (C) 2002 Elsevier Science B.V. All rights reserved.

Close packed transitions in slit-shaped pores: Density functional theory study of methane adsorption capacity in carbon

An important feature of improving lattice gas models and classical isotherms is the incorporation of a pore size dependent capacity, which has hitherto been overlooked. In this paper, we develop a model for predicting the temperature dependent variation in capacity with pore size. The model is based on the analysis of a lattice gas model using a density functional theory approach at the close packed limit. Fluid-fluid and solid-fluid interactions are modeled by the Lennard-Jones 12-6 potential and Steele's 10-4-3, potential respectively. The capacity of methane in a slit-shaped carbon pore is calculated from the characteristic parameters of the unit cell, which are extracted by minimizing the grand potential of the unit cell. The capacities predicted by the proposed model are in good agreement with those obtained from grand canonical Monte Carlo simulation, for pores that can accommodate up to three adsorbed layers. Single particle and pair distributions exhibit characteristic features that correspond to the sequence of buckling and rhombic transitions that occur as the slit pore width is increased. The model provides a useful tool to model continuous variation in the microstructure of an adsorbed phase, namely buckling and rhombic transitions, with increasing pore width. (C) 2002 American Institute of Physics.

Density functional theory analysis of the influence of pore wall heterogeneity on adsorption in carbons

Density functional theory for adsorption in carbons is adapted here to incorporate a random distribution of pore wall thickness in the solid, and it is shown that the mean pore wall thickness is intimately related to the pore size distribution characteristics. For typical carbons the pore walls are estimated to comprise only about two graphene layers, and application of the modified density functional theory approach shows that the commonly used assumption of infinitely thick walls can severely affect the results for adsorption in small pores under both supercritical and subcritical conditions. Under supercritical conditions the Henry's law coefficient is overpredicted by as much as a factor of 2, while under subcritical conditions pore wall heterogeneity appears to modify transitions in small pores into a sequence of smaller ones corresponding to pores with different wall thicknesses. The results suggest the need to improve current pore size distrubution analysis methods to allow for pore wall heterogeneity. The density functional theory is further extended here to allow for interpore adsorbate interactions, and it appears that these interaction are negligible for small molecules such as nitrogen but significant for more strongly interacting heavier molecules such as butane, for which the traditional independent pore model may not be adequate.

Experimental investigation of contaminant transport in coastal groundwater

Contaminant transport in coastal aquifers is of increasing interest since, with the development of coastal areas, contaminants from surface sources may enter coastal aquifers and pollute the groundwater flow. Coastal groundwater flow is complicated because of the presence of a freshwater-saltwater diffusion zone and the tidal variation of sea level at the seaward end. This paper investigates experimentally the behaviour of contaminant plumes with different densities in an unconfined coastal aquifer. Experiments were performed in a flow tank filled with glass beads as the porous medium. Results show that the dense contaminant has a more diffusive front than the less dense one in the seaward direction towards the coastline. The plume becomes more diffusive when it travels closer to the saltwater interface. On the contrary, the less dense contaminant presents a relatively sharp outline. It tends to migrate in the upper portion of the aquifer and exits in a concentrated manner over a small discharge area at the coastline, not further seaward under the sea. Non-dimensional parameters show that instabilities occur in our experiments for a density difference of 1.2% or larger between the contaminant and the ambient water. The experimental results provide guidance for field monitoring and numerical modelling. (C) 2002 Elsevier Science Ltd. All rights reserved.

Polyethylene flow prediction with a differential multi-mode Pom-Pom model

We report the first steps of a collaborative project between the University of Queensland, Polyflow, Michelin, SK Chemicals, and RMIT University; on simulation, validation and application of a recently introduced constitutive model designed to describe branched polymers. Whereas much progress has been made on predicting the complex flow behaviour of many - in particular linear - polymers, it sometimes appears difficult to predict simultaneously shear thinning and extensional strain hardening behaviour using traditional constitutive models. Recently a new viscoelastic model based on molecular topology, was proposed by McLeish and Larson (1998). We explore the predictive power of a differential multi-mode version of the pom-pom model for the flow behaviour of two commercial polymer melts: a (long-chain branched) low-density polyethylene (LDPE) and a (linear) high-density polyethylene (HDPE). The model responses are compared to elongational recovery experiments published by Langouche and Debbaut (1999), and start-up of simple shear flow, stress relaxation after simple and reverse step strain experiments carried out in our laboratory.

Repeated three-dimensional magnetic resonance imaging of atherosclerosis development in innominate arteries of low-density lipoprotein receptor-knockout mice

Background-In vivo methods to evaluate the size and composition of atherosclerotic lesions in animal models of atherosclerosis would assist in the testing of antiatherosclerotic drugs. We have developed an MRI method of detecting atherosclerotic plaque in the major vessels at the base of the heart in low-density lipoprotein (LDL) receptor-knockout (LDLR-/-) mice on a high-fat diet. Methods and Results-Three-dimensional fast spin-echo magnetic resonance images were acquired at 7 T by use of cardiac and respiratory triggering, with approximate to140-mum isotropic resolution, over 30 minutes. Comparison of normal and fat-suppressed images from female LDLR-/- mice I week before and 8 and 12 weeks after the transfer to a high-fat diet allowed visualization and quantification of plaque development in the innominate artery in vivo. Plaque mean cross-sectional area was significantly greater at week 12 in the LDLR-/- mice (0.14+/-0.086 mm(2) [mean+/-SD]) than in wild-type control mice on a normal diet (0.017+/-0.031 mm(2), p

A comparative study of carbon gasification with O-2 and CO2 by density functional theory calculations

A comparative study of carbon gasification with O-2 and CO2 was conducted by using density functional theory calculations. It was found that the activation energy and the number of active sites in carbon gasification reactions are significantly affected by both the capacity and manner of gas chemisorption. O-2 has a strong adsorption capacity and the dissociative chemisorption of O-2 is thermodynamically favorable on either bare carbon surface or even isolated edge sites. As a result, a large number of semiquinone and o-quinone oxygen can be formed indicating a significant increase in the number of active sites. Moreover, the weaker o-quinone C-C bonds can also drive the reaction forward at (ca. 30%) lower activation energy. Epoxy oxygen forms under relatively high O-2 pressure, and it can only increase the number of active sites, not further reduce the activation energy. CO2 has a lower adsorption capacity. Dissociative chemisorption of CO2 can only occur on two consecutive edge sites and o-quinone oxygen formed from CO2 chemisorption is negligible, let alone epoxy oxygen. Therefore, CO2-carbon reaction needs (ca 30%) higher activation energy. Furthermore, the effective active sites are also reduced by the manner Of CO2 chemisorption. A combination of the higher activation energy and the fewer active sites leads to the much lower reaction rate Of CO2-carbon.

Composition and energy density of eggs from two species of freshwater turtle with twofold ranges in egg size

Lipid, protein, ash, carbohydrate and water content and energy density of eggs were measured from different clutches over a range of egg size in two species of freshwater turtle. Dry egg contents consisted of protein (54-60%), lipid (25-31%) and ash (5-6%) while carbohydrate was found to be negligible (