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.

The efficiency and condition of oysters and macroalgae used as biological filters of shrimp pond effluent

Current shrimp pond management practices generally result in elevated concentrations of nutrients, suspended solids, bacteria and phytoplankton compared with the influent water. Concerns about adverse environmental impacts caused by discharging pond effluent directly into adjacent waterways have prompted the search for cost-effective methods of effluent treatment. One potential method of effluent treatment is the use of ponds or raceways stocked with plants or animals that act as natural biofilters by removing waste nutrients. In addition to improving effluent water quality prior to discharge, the use of natural biofilters provides a method for capturing otherwise wasted nutrients. This study examined the potential of the native oyster, Saccostrea commercialis (Iredale and Roughley) and macroalgae, Gracilaria edulis (Gmelin) Silva to improve effluent water quality from a commercial Penaeus japonicus (Bate) shrimp farm, A system of raceways was constructed to permit recirculation of the effluent through the oysters to maximize the filtration of bacteria, phytoplankton and total suspended solids. A series of experiments was conducted to test the ability of oysters and macroalgae to improve effluent water quality in a flow-through system compared with a recirculating system. In the flow-through system, oysters reduced the concentration of bacteria to 35% of the initial concentration, chlorophyll a to 39%, total particulates (2.28-35.2 mum) to 29%, total nitrogen to 66% and total phosphorus to 56%. Under the recirculating flow regime, the ability of the oysters to improve water quality was significantly enhanced. After four circuits, total bacterial numbers were reduced to 12%, chlorophyll a to 4%, and total suspended solids to 16%. Efforts to increase biofiltration by adding additional layers of oyster trays and macroalgae-filled mesh bags resulted in fouling of the lower layers causing the death of oysters and senescence of macroalgae. Supplementary laboratory experiments were designed to examine the effects of high effluent concentrations of suspended particulates on the growth and condition of oysters and macroalgae. The results demonstrated that high concentrations of particulates inhibited growth and reduced the condition of oysters and macroalgae. Allowing the effluent to settle before biofiltration improved growth and reduced signs of stress in the oysters and macroalgae. A settling time of 6 h reduced particulates to a level that prevented fouling of the oysters and macroalgae.

Rebreathing potential of infant mattresses and bedcovers

Objective : To establish the CO2 dispersion and retention properties of some mattresses and bed coverings commercially available in Australia. Methods : Five mattresses were studied in (i) an in vivo model in which an infant's head was covered by a headbox, rebreathing was allowed to occur, and the final steady state CO2 concentration was measured; and (ii) an in vitro model in which 5% CO2 in a headbox was allowed to disperse, and the time taken for the concentration to reach 1% was measured. Five types of bedcover were studied in (i) an in vivo model in which an infant's head was covered by a bedcover and the final steady state CO2 concentration was measured; and (ii) an in vitro model in which 5% CO2 under a bedcover was allowed to disperse, and the time taken for the concentration to reach 1% was measured. Results : The steady state CO2 concentrations ranged from 0.6% to 3.0% for the mattresses (P < 0.05). The time for CO2 to disperse ranged from 5.5 min to 30.4 min (P < 0.05). Steady state CO2 concentrations ranged from 2.5% to 3.6% for the bedcoverings (P > 0.05). The time for CO2 to disperse ranged from 5.4 min to 7.7 min (P > 0.05). Conclusions : Some commercial cot mattresses and bedcoverings allow high concentrations of CO2 to accumulate in rebreathing environments. Some mattress types studied were more diffusive to CO2 , whereas there was no difference between the bedcovers studied. This may have implications for vulnerable infants at risk of sudden infant death syndrome.

Control of chaotic instability in a dual-spin spacecraft with dissipation using energy methods

Control of chaotic instability in a rotating multibody system in the form of a dual-spin spacecraft with an axial nutational damper is achieved using an algorithm derived using energy methods. The control method is implemented on two realistic spacecraft parameter configurations which have been found to exhibit chaotic instability when a sinusoidally varying torque is applied to the spacecraft for a range of forcing amplitudes and frequencies. Such a torque, in practice, may arise under malfunction of the control system or from an unbalanced rotor. Chaotic instabilities arising from these torques could introduce uncertainties and irregularities into a spacecraft's attitude and consequently impair pointing accuracy. The control method is formulated from nutational stability results derived using an energy sink approximation for a dual-spin spacecraft with an asymmetric platform and axisymmetric rotor. The effectiveness of the control method is shown numerically and the results are studied by means of time history, phase space, Poincare map, Lyapunov characteristic exponents and Bifurcation diagrams.

Improvements to peroxide oxidation methods for analysing sulfur in acid sulfate soils

Improvements to peroxide oxidation methods for analysing acid sulfate soils (ASS) are introduced. The soil solution ratio has been increased to 1 : 40, titrations are performed in suspension, and the duration of the peroxide digest stage is substantially shortened. For 9 acid sulfate soils, the peroxide oxidisable sulfur value obtained using the improved method was compared with the reduced inorganic sulfur result obtained using the chromium reducible sulfur method. Their regression was highly significant, the slope of the regression line was not significantly different (P = 0.05) from unity, and the intercept not significantly different from zero. A complete sulfur budget for the improved method showed there was no loss of sulfur as has been reported for earlier peroxide oxidation techniques. When soils were very finely ground, efficient oxidation of sulfides was achieved, despite the milder digestion conditions. Highly sulfidic and organic soils were shown to be the most difficult to analyse using either the improved method or the chromium method. No single analytical method can be universally applied to all ASS, rather a suite of methods is necessary for a thorough understanding of many ASS. The improved peroxide method, in combination with the chromium method and the 4 M HCl extraction, form a sound platform for informed decision making on the management of acid sulfate soils.

Body mass index-related human adipocyte agouti expression is sex-specific but not depot-specific

Objective: To determine if human adipocyte agouti signal protein (ASIP) mRNA expression is associated with obesity and is gender and/or depot specific. Research Methods and Procedures: Subjects included 8 men (64 +/- 3 years) and 14 women (56 +/- 15 years) undergoing elective abdominal surgery. ASIP mRNA levels in isolated omental and subcutaneous abdominal adipocytes were measured by quantitative reverse transcription polymerase chain reaction. Results: No significant depot difference was observed between genders; ASIP mRNA levels of omental and subcutaneous abdominal adipocytes were pooled for this analysis. BMI and ASIP gene expression were negatively correlated in men (p = -0.70; p < 0.05), whereas a positive relationship was observed in women (p = 0.48; p < 0.05). No significant difference was observed in age, body weight, body mass index (BMI), and waist circumference between groups. Hip circumference was significantly higher in women than in men (p < 0.05). Also, no significant difference in ASIP mRNA expression was observed between men and women, regardless of the fat depot. Discussion: These results show that men and women of similar age and BMI present similar ASIP mRNA levels in omental and subcutaneous abdominal adipocytes. However, a sexual dimorphism exists in the relationship between ASIP expression and BMI. If ASIP is involved in appetite regulation or energy homeostasis in humans, this observation may contribute to the recognized differences in these parameters between men and women.

Inhibition of adventitious rooting in Backhousia citriodora F. Muell cuttings correlate with the concentration of essential oil

Backhousia citriodora is typical of the many commercially valuable woody Australian Myrtaceae species that are recalcitrant in forming adventitious roots from cuttings after maturation. A series of experiments were conducted to identify an endogenous rooting inhibitor in line with established criteria. Endogenous levels of citral were correlated with the rooting capacities of juvenile versus mature, and easy- versus difficult-to-root genotypes of B. citriodora, in both winter and summer. The biological activity of citral was confirmed in bioassays on mung beans and easy-to-root B. citriodora seedlings. Evidence of a common mechanism of root inhibition with other species in the Myrtaceae and the role of action of citral are discussed.

Co-crystallization of sucrose at high concentration in the presence of glucose and fructose

Co-crystallization of sucrose from a highly concentrated sucrose syrup (less than or equal to 7% moisture, w/w) at 131 degreesC with 0, 5, 10, 15, and 20% of fructose, glucose, or a mixture of fructose and glucose was investigated. The crystallization of sucrose was delayed in presence of these lower molecular weight sugars. The DSC melting endotherm of cocrystallized samples exhibited a decrease in crystalline sucrose in the sample as a function of increased level of glucose and fructose. The mechanical strength of co-crystallized granules was found to be related to the moisture content and the amount of glucose or fructose content in the sample. The samples containing 10, 15, and 20% glucose in co-crystallized product demonstrated crystallization of glucose in its monohydrate form during 1 mo of storage.

Predictor-corrector methods of Runge-Kutta type for stochastic differential equations

In this paper we construct predictor-corrector (PC) methods based on the trivial predictor and stochastic implicit Runge-Kutta (RK) correctors for solving stochastic differential equations. Using the colored rooted tree theory and stochastic B-series, the order condition theorem is derived for constructing stochastic RK methods based on PC implementations. We also present detailed order conditions of the PC methods using stochastic implicit RK correctors with strong global order 1.0 and 1.5. A two-stage implicit RK method with strong global order 1.0 and a four-stage implicit RK method with strong global order 1.5 used as the correctors are constructed in this paper. The mean-square stability properties and numerical results of the PC methods based on these two implicit RK correctors are reported.

An atomistic simulation of solid state sintering using Monte Carlo methods

This paper presents results on the simulation of the solid state sintering of copper wires using Monte Carlo techniques based on elements of lattice theory and cellular automata. The initial structure is superimposed onto a triangular, two-dimensional lattice, where each lattice site corresponds to either an atom or vacancy. The number of vacancies varies with the simulation temperature, while a cluster of vacancies is a pore. To simulate sintering, lattice sites are picked at random and reoriented in terms of an atomistic model governing mass transport. The probability that an atom has sufficient energy to jump to a vacant lattice site is related to the jump frequency, and hence the diffusion coefficient, while the probability that an atomic jump will be accepted is related to the change in energy of the system as a result of the jump, as determined by the change in the number of nearest neighbours. The jump frequency is also used to relate model time, measured in Monte Carlo Steps, to the actual sintering time. The model incorporates bulk, grain boundary and surface diffusion terms and includes vacancy annihilation on the grain boundaries. The predictions of the model were found to be consistent with experimental data, both in terms of the microstructural evolution and in terms of the sintering time. (C) 2002 Elsevier Science B.V. All rights reserved.

Comparison of an ELISA and an LC/MS/MS method for measuring tacrolimus concentrations and making dosage decisions in transplant recipients

This study compared an enzyme-linked immunosorbent assay (ELISA) to a liquid chromatography-tandem mass spectrometry (LC/MS/MS) technique for measurement of tacrolimus concentrations in adult kidney and liver transplant recipients, and investigated how assay choice influenced pharmacokinetic parameter estimates and drug dosage decisions. Tacrolimus concentrations measured by both ELISA and LC/MS/MS from 29 kidney (n = 98 samples) and 27 liver (n = 97 samples) transplant recipients were used to evaluate the performance of these methods in the clinical setting. Tacrolimus concentrations measured by the two techniques were compared via regression analysis. Population pharmacokinetic models were developed independently using ELISA and LC/MS/MS data from 76 kidney recipients. Derived kinetic parameters were used to formulate typical dosing regimens for concentration targeting. Dosage recommendations for the two assays were compared. The relation between LC/MS/MS and ELISA measurements was best described by the regression equation ELISA = 1.02 . (LC/MS/MS) + 0.14 in kidney recipients, and ELISA = 1.12 . (LC/MS/MS) - 0.87 in liver recipients. ELISA displayed less accuracy than LC/MS/MS at lower tacrolimus concentrations. Population pharmacokinetic models based on ELISA and LC/MS/MS data were similar with residual random errors of 4.1 ng/mL and 3.7 ng/mL, respectively. Assay choice gave rise to dosage prediction differences ranging from 0% to 30%. ELISA measurements of tacrolimus are not automatically interchangeable with LC/MS/MS values. Assay differences were greatest in adult liver recipients, probably reflecting periods of liver dysfunction and impaired biliary secretion of metabolites. While the majority of data collected in this study suggested assay differences in adult kidney recipients were minimal, findings of ELISA dosage underpredictions of up to 25% in the long term must be investigated further.