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.

Automatic construction of explicit R matrices for the one-parameter families of irreducible typical highest weight (0(m)vertical bar alpha(n)) representations of U-q[gl(m vertical bar n)]

We detail the automatic construction of R matrices corresponding to (the tensor products of) the (O-m\alpha(n)) families of highest-weight representations of the quantum superalgebras Uq[gl(m\n)]. These representations are irreducible, contain a free complex parameter a, and are 2(mn)-dimensional. Our R matrices are actually (sparse) rank 4 tensors, containing a total of 2(4mn) components, each of which is in general an algebraic expression in the two complex variables q and a. Although the constructions are straightforward, we describe them in full here, to fill a perceived gap in the literature. As the algorithms are generally impracticable for manual calculation, we have implemented the entire process in MATHEMATICA; illustrating our results with U-q [gl(3\1)]. (C) 2002 Published by Elsevier Science B.V.

Growth and compaction behaviour of copper concentrate granules in a rotating drum

In order to understand the growth and compaction behaviour of chalcopyrite (copper concentrate), batch granulation tests were carried out using a rotating drum. The granule growth exhibited induction-type behaviour, as defined by Iveson and Litster [AIChE J. 44 (1998) 15 10]. There were two consecutive stages during granulation: the induction stage, during which the granules are gradually being compacted and little or no growth occurs, and the rapid growth stage, which starts when the granules have become surface wet and are rapidly growing. In agreement with earlier findings. an increased amount of binder liquid shortened the induction time. The compaction behaviour was also investigated. A displaced volume method was adopted to determine the porosity of the granules. It was shown that this technique had a limitation as it was unable to detect the reduction of the volumes of the granule pores after the granules had become surface wet. Due to this, some of the measurements were not suited for fitting a three-parameter empirical model. Attempts were made to determine whether the rapid growth stage started with the pore saturation exceeding a certain critical value, but due to the scatter in the porosity measurements and the fact that some of the measurements could not be used, it was not possible to determine a critical pore saturation, However, the porosity measurements clearly demonstrated that the porosity of the granules decreased during the induction stage of an experiment and that when rapid growth occurred, the granules had a pore saturation was around 0.85. This value was slightly lower than unity, which is most likely due to trapped air bubbles. (C) 2002 Published by Elsevier Science B.V.

Vacancy solution theory for binary adsorption equilibria in heterogeneous carbon

A heterogeneous modified vacancy solution model of adsorption developed is evaluated. The new model considers the adsorption process through a mass-action law and is thermodynamically consistent, while maintaining the simplicity in calculation of multicomponent adsorption equilibria, as in the original vacancy solution theory. It incorporates the adsorbent heterogeneity through a pore-width-related potential energy, represented by Steele's 10-4-3 potential expression. The experimental data of various hydrocarbons, CO2 and SO2 on four different activated carbons - Ajax, Norit, Nuxit, and BPL - at multiple temperatures over a wide range of pressures were studied by the heterogeneous modified VST model to obtain the isotherm parameters and micropore-size distribution of carbons. The model successfully correlates the single-component adsorption equilibrium data for all compounds studied on various carbons. The fitting results for the vacancy occupancy parameter are consistent with the pressure change on different carbons, and the effect of pore heterogeneity is important in adsorption at elevated pressure. It predicts binary adsorption equilibria better than the IAST scheme, reflecting the significance of molecular size nonideality.

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.

Beyond hearing aid fitting: improving communication for older adults

Many older adults with hearing impairment continue to have substantial communication difficulties after being fitted with hearing aids, and many do not choose to wear hearing aids. Two group communication education programs aimed at such older people are described. The 'Keep on Talking' program has a health promotion focus, and is aimed at maintaining communication for older adults living in the community. An experimental group (n=120) attended the program, and a control group (n=130) received a communication assessment but no intervention. Significant improvements were found in the experimental participants in terms of knowledge about communication changes with age and about strategies to maintain communication skills. At the follow-up evaluation at 1 year, 45% of the experimental group, compared to 10% of the control group, had acted to improve their communication skills. The 'Active Communication Education' program focuses on the development of problem-solving strategies to improve communication in everyday life situations. Preliminary outcomes have been assessed on a small scale (n=14) to date. It is concluded that communication programs represent an important adjunct to, or supplement for, the traditional approach that focuses on hearing aid fitting.

Effect of Co addition on the lattice parameter, electrical conductivity and sintering of gadolinia-doped ceria

Co-sintering aid has been added to Ce1.9Gd0.1O1.95 (CGO) by treating a commercial powder with Co(NO3)(2) (COCGO), X-ray diffraction (XRD) measurements of lattice parameter indicated that the Co was located on the CGO particle surface after calcination at 650 degreesC. After heat treatment at temperatures above 650 degreesC, the room temperature lattice parameter of CGO was found to increase, indicating redistribution of the Gd. Compared to CGO, the lattice parameter of CGO + 2 cation% Co (2CoCGO) was lower for a given temperature (650-1100 degreesC), A.C. impedance revealed that the lattice conductivity of 2CoCGO was enhanced when densified at lower temperatures, Transmission electron microscopy (TEM) showed that, even after sintering for 4 h at 980 degreesC, most of the Co was located at grain boundaries. (C) 2002 Published by Elsevier Science B.V.

Link invariants associated with gauge equivalent solutions of the Yang-Baxter equation: The one-parameter family of minimal typical representations of Uq[gl(2|1)]

In this paper we investigate the construction of state models for link invariants using representations of the braid group obtained from various gauge choices for a solution of the trigonometric Yang-Baxter equation. Our results show that it is possible to obtain invariants of regular isotopy (as defined by Kauffman) which may not be ambient isotopic. We illustrate our results with explicit computations using solutions of the trigonometric Yang-Baxter equation associated with the one-parameter family of minimal typical representations of the quantum superalgebra U-q,[gl(2/1)]. We have implemented MATHEMATICA code to evaluate the invariants for all prime knots up to 10 crossings.

Sistema de Informação Geográfica: uma arma estratégica para o Database Marketing

Os Sistemas de Informações Geográficas (SIG) tornaram-se disponíveis para uso em computadores de mesa, sendo extremamente adaptáveis para as atividades de marketing de relacionamento. Uma característica desses sistemas é que permitem uma análise de dados de forma praticamente inviável com outras ferramentas. A sua utilização depende da existência de mapas e dados possivelmente caros e difíceis de obter. No Brasil, seu uso vai popularizar-se nos próximos dois anos.

Dynamical behaviour on the parameter space: new populational growth models proportional to beta densities

We present new populational growth models, generalized logistic models which are proportional to beta densities with shape parameters p and 2, where p > 1, with Malthusian parameter r. The complex dynamical behaviour of these models is investigated in the parameter space (r, p), in terms of topological entropy, using explicit methods, when the Malthusian parameter r increases. This parameter space is split into different regions, according to the chaotic behaviour of the models.

Quantificação por imagem ponderada em difusão (DWI) das lesões tumorais da mama

Introdução – A técnica de Difusão por Ressonância Magnética (RM), ao avaliar o movimento das moléculas de água nos tecidos, permite inferir sobre a arquitetura dos mesmos, em particular relativamente à celularidade, volume celular e permeabilidade das membranas. O Coeficiente de Difusão Aparente (ADC) é um parâmetro quantificável da imagem ponderada em difusão (DWI). A sua análise poderá fornecer informação clínica adicional sobre estas lesões, sobretudo em relação à sua caracterização histológica. Objetivos – Caracterizar e diferenciar tipos e alguns subtipos de lesões mamárias através da análise do ADC. Metodologia – 20 Mulheres com 23 lesões mamárias foram submetidas a RM mamária: 3 lesões benignas (3 Fibroadenomas-FA) e 20 malignas (16 Carcinomas Ductais Invasivos-CDI, 2 Carcinomas Ductais In Situ-CDIS e 2 Carcinomas Invasivos Lobulares-CLI). Num equipamento 1.5T aplicou-se uma sequência de Difusão (b=0,50,250,500,750,1000 s/mm2). Obteve-se o ADC através do ajuste exponencial da intensidade de sinal das lesões em função do valor de b, fazendo-se corresponder os valores de ADC à respetiva caracterização histológica e posterior comparação com a literatura. Resultados e Discussão – As lesões malignas apresentaram ADCs significativamente (p=0,014) inferiores [(0,94±0,22)x10-3 mm2/s] aos das benignas [(1,43±0,25)x10-3 mm2/s]. A justificação pode residir no aumento da celularidade e consequente restrição da Difusão que se observa nas lesões malignas. Os CDI apresentaram ADCs baixos [(0,88±0,17)x10-3 mm2/s], enquanto que os CDIS apresentaram ADCs mais elevados [(1,33±0,29)x10-3 mm2/s]. Estes resultados estão de acordo com o facto dos CDIS estarem limitados aos ductos mamários, mantendo-se menos alterada a estrutura do tecido adjacente e resultando numa menor restrição à difusão que nos CDI. Verificaram-se diferenças significativas entre FA e CDI (p=0,010) e entre CDI e CDIS (p=0,049). Conclusões – O ADC possibilita a diferenciação entre lesões mamárias benignas e malignas, bem como entre alguns tipos histológicos. O desenvolvimento deste conceito pode representar um avanço no papel da RM na avaliação destas neoplasias. ABSTRACT - Introduction – The Magnetic Resonance (MR) diffusion technique measures the movement of water molecules in tissues. Therefore, it provides useful information about tissue architecture, specially regarding tissue cellularity, cell volume and membrane permeability. The quantification of diffusion weighted imaging (DWI) data is done by measuring the so-called. Apparent Diffusion Coefficient (ADC). This parameter provides additional clinical information about breast lesions, and potentially allows for in-vivo histological characterization. Purpose – To characterize and differentiate breast lesions through ADC analysis. Methodology – The study comprised 20 women, with 23 breast lesions: 3 benign lesions - 3 Fibroadenomas (FA); and 20 malignant - 16 Invasive Ductal Carcinomas (CDI), 2 Ductal Carcinomas In Situ (CDIS), 2 Invasive Lobular Carcinoma (CLI). On a 1.5T equipment a diffusion-weighted sequence with 6 b-values (b=0,50,250,500,750,1000 s/mm2) was used to examine the patients. ADC was obtained by fitting an exponential to data of lesion signal intensity vs. b values. A correspondence of ADC values to histological lesion characterization was done and finally, the results were comparison with the literature. Results and Discussion – Malignant lesions showed inferior ADCs significantly (p=0.014) lower ((0.94±0.22)x10-3 mm2/s) than the benign lesions ((1.43±0.25)x10-3 mm2/s). This may be associated to increasead cellularity in malignant lesions, which result in higher tissue restriction to diffusion. CDI showed low ADC values ((0.88±0.17)x10-3 mm2/s), while the CDIS showed higher ADCs ((1.33±0.29)x10-3 mm2/s). These results agree with the fact that CDIS are limited to mammary ducts, maintaining a less altered neighboring tissue structure, which results in a lower restriction to diffusion than observed in CDI. Significant differences between FA and CDI (p=0.010) and between CDI and CDIS (p=0.049) were observed. Conclusion – The ADC parameter is able to differentiate between malignant and benign breast lesions, as well as between some histological types.

Multi-apprentice learning for meta-heuristics parameter tuning in a multi agent scheduling system

The scheduling problem is considered in complexity theory as a NP-hard combinatorial optimization problem. Meta-heuristics proved to be very useful in the resolution of this class of problems. However, these techniques require parameter tuning which is a very hard task to perform. A Case-based Reasoning module is proposed in order to solve the parameter tuning problem in a Multi-Agent Scheduling System. A computational study is performed in order to evaluate the proposed CBR module performance.