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.

Conductivity/microstructure study of a pyrochlore-type ordered phase in the compound Zr0.75Ce0.08Nd0.17O1.92

The compound Zr0.75Ce0.08Nd0.17O1.92 was investigated as part of a much larger electrical conductivity/microstructure study of the systems ZrO2-CeO2-M2O3 (where M=Nd, Sm, ..., Yb) [Solid State Ionics (2002)]. Electrical conductivity measurements performed in air at 800 degreesC showed significant conductivity degradation over a period of 200 h. Investigation of the annealed and as-fired specimens by ATEM revealed the presence of an emerging, ordered pyrochlore-type phase within the Zr0.75Ce0.08Nd0.17O1.92 defect-fluorite solid solution at much lower dopant levels than observed previously for zirconia binary systems. (C) 2002 Elsevier Science B.V. All rights reserved.

Predicting the HER2 status of breast cancer from basic histopathology data: an analysis of 1500 breast cancers as part of the HER2000 International Study

The tests that are currently available for the measurement of overexpression of the human epidermal growth factor-2 (HER2) in breast cancer have shown considerable problems in accuracy and interlaboratory reproducibility. Although these problems are partly alleviated by the use of validated, standardised 'kits', there may be considerable cost involved in their use. Prior to testing it may therefore be an advantage to be able to predict from basic pathology data whether a cancer is likely to overexpress HER2. In this study, we have correlated pathology features of cancers with the frequency of HER2 overexpression assessed by immunohistochemistry (IHC) using HercepTest (Dako). In addition, fluorescence in situ hybridisation (FISH) has been used to re-test the equivocal cancers and interobserver variation in assessing HER2 overexpression has been examined by a slide circulation scheme. Of the 1536 cancers, 1144 (74.5%) did not overexpress HER2. Unequivocal overexpression (3+ by IHC) was seen in 186 cancers (12%) and an equivocal result (2+ by IHC) was seen in 206 cancers (13%). Of the 156 IHC 3+ cancers for which complete data was available, 149 (95.5%) were ductal NST and 152 (97%) were histological grade 2 or 3. Only 1 of 124 infiltrating lobular carcinomas (0.8%) showed HER2 overexpression. None of the 49 'special types' of carcinoma showed HER2 overexpression. Re-testing by FISH of a proportion of the IHC 2+ cancers showed that only 25 (23%) of those assessable exhibited HER2 gene amplification, but 46 of the 47 IHC 3+ cancers (98%) were confirmed as showing gene amplification. Circulating slides for the assessment of HER2 score showed a moderate level of agreement between pathologists (kappa 0.4). As a result of this study we would advocate consideration of a triage approach to HER-2 testing. Infiltrating lobular and special types of carcinoma may not need to be routinely tested at presentation nor may grade 1 NST carcinomas in which only 1.4% have been shown to overexpress HER2. Testing of these carcinomas may be performed when HER2 status is required to assist in therapeutic or other clinical/prognostic decision-making. The highest yield of HER2 overexpressing carcinomas is seen in the grade 3 NST subgroup in which 24% are positive by IHC. (C) 2003 Elsevier Science Ltd. All rights reserved.

Capillaries within compartments: Microvascular interpretation of dynamic positron emission tomography data

Measurement of exchange of substances between blood and tissue has been a long-lasting challenge to physiologists, and considerable theoretical and experimental accomplishments were achieved before the development of the positron emission tomography (PET). Today, when modeling data from modern PET scanners, little use is made of earlier microvascular research in the compartmental models, which have become the standard model by which the vast majority of dynamic PET data are analysed. However, modern PET scanners provide data with a sufficient temporal resolution and good counting statistics to allow estimation of parameters in models with more physiological realism. We explore the standard compartmental model and find that incorporation of blood flow leads to paradoxes, such as kinetic rate constants being time-dependent, and tracers being cleared from a capillary faster than they can be supplied by blood flow. The inability of the standard model to incorporate blood flow consequently raises a need for models that include more physiology, and we develop microvascular models which remove the inconsistencies. The microvascular models can be regarded as a revision of the input function. Whereas the standard model uses the organ inlet concentration as the concentration throughout the vascular compartment, we consider models that make use of spatial averaging of the concentrations in the capillary volume, which is what the PET scanner actually registers. The microvascular models are developed for both single- and multi-capillary systems and include effects of non-exchanging vessels. They are suitable for analysing dynamic PET data from any capillary bed using either intravascular or diffusible tracers, in terms of physiological parameters which include regional blood flow. (C) 2003 Elsevier Ltd. All rights reserved.

An EM-based Semi-Parametric Mixture Model Approach to the Regression Analysis of Competing-Risks Data

We consider a mixture model approach to the regression analysis of competing-risks data. Attention is focused on inference concerning the effects of factors on both the probability of occurrence and the hazard rate conditional on each of the failure types. These two quantities are specified in the mixture model using the logistic model and the proportional hazards model, respectively. We propose a semi-parametric mixture method to estimate the logistic and regression coefficients jointly, whereby the component-baseline hazard functions are completely unspecified. Estimation is based on maximum likelihood on the basis of the full likelihood, implemented via an expectation-conditional maximization (ECM) algorithm. Simulation studies are performed to compare the performance of the proposed semi-parametric method with a fully parametric mixture approach. The results show that when the component-baseline hazard is monotonic increasing, the semi-parametric and fully parametric mixture approaches are comparable for mildly and moderately censored samples. When the component-baseline hazard is not monotonic increasing, the semi-parametric method consistently provides less biased estimates than a fully parametric approach and is comparable in efficiency in the estimation of the parameters for all levels of censoring. The methods are illustrated using a real data set of prostate cancer patients treated with different dosages of the drug diethylstilbestrol. Copyright (C) 2003 John Wiley Sons, Ltd.

Effect of changing data collection parameters on statistical motor unit number estimates

The effect of number of samples and selection of data for analysis on the calculation of surface motor unit potential (SMUP) size in the statistical method of motor unit number estimates (MUNE) was determined in 10 normal subjects and 10 with amyotrophic lateral sclerosis (ALS). We recorded 500 sequential compound muscle action potentials (CMAPs) at three different stable stimulus intensities (10–50% of maximal CMAP). Estimated mean SMUP sizes were calculated using Poisson statistical assumptions from the variance of 500 sequential CMAP obtained at each stimulus intensity. The results with the 500 data points were compared with smaller subsets from the same data set. The results using a range of 50–80% of the 500 data points were compared with the full 500. The effect of restricting analysis to data between 5–20% of the CMAP and to standard deviation limits was also assessed. No differences in mean SMUP size were found with stimulus intensity or use of different ranges of data. Consistency was improved with a greater sample number. Data within 5% of CMAP size gave both increased consistency and reduced mean SMUP size in many subjects, but excluded valid responses present at that stimulus intensity. These changes were more prominent in ALS patients in whom the presence of isolated SMUP responses was a striking difference from normal subjects. Noise, spurious data, and large SMUP limited the Poisson assumptions. When these factors are considered, consistent statistical MUNE can be calculated from a continuous sequence of data points. A 2 to 2.5 SD or 10% window are reasonable methods of limiting data for analysis. Muscle Nerve 27: 320–331, 2003

Model-based clustering in gene expression microarrays: An application to breast cancer data

In microarray studies, the application of clustering techniques is often used to derive meaningful insights into the data. In the past, hierarchical methods have been the primary clustering tool employed to perform this task. The hierarchical algorithms have been mainly applied heuristically to these cluster analysis problems. Further, a major limitation of these methods is their inability to determine the number of clusters. Thus there is a need for a model-based approach to these. clustering problems. To this end, McLachlan et al. [7] developed a mixture model-based algorithm (EMMIX-GENE) for the clustering of tissue samples. To further investigate the EMMIX-GENE procedure as a model-based -approach, we present a case study involving the application of EMMIX-GENE to the breast cancer data as studied recently in van 't Veer et al. [10]. Our analysis considers the problem of clustering the tissue samples on the basis of the genes which is a non-standard problem because the number of genes greatly exceed the number of tissue samples. We demonstrate how EMMIX-GENE can be useful in reducing the initial set of genes down to a more computationally manageable size. The results from this analysis also emphasise the difficulty associated with the task of separating two tissue groups on the basis of a particular subset of genes. These results also shed light on why supervised methods have such a high misallocation error rate for the breast cancer data.

Localization of a novel melanoma susceptibility locus to 1p22

Over the past 20 years, the incidence of cutaneous malignant melanoma (CMM) has increased dramatically worldwide. A positive family history of the disease is among the most established risk factors for CMM; it is estimated that 10% of CMM cases result from an inherited predisposition. Although mutations in two genes, CDKN2A and CDK4, have been shown to confer an increased risk of CMM, they account for only 20%-25% of families with multiple cases of CMM. Therefore, to localize additional loci involved in melanoma susceptibility, we have performed a genomewide scan for linkage in 49 Australian pedigrees containing at least three CMM cases, in which CDKN2A and CDK4 involvement has been excluded. The highest two-point parametric LOD score (1.82; recombination fraction [theta] 0.2) was obtained at D1S2726, which maps to the short arm of chromosome 1 (1p22). A parametric LOD score of 4.65 (theta = 0) and a nonparametric LOD score of 4.19 were found at D1S2779 in nine families selected for early age at onset. Additional typing yielded seven adjacent markers with LOD scores 13 in this subset, with the highest parametric LOD score, 4.95 (theta = 0) ( nonparametric LOD score 5.37), at D1S2776. Analysis of 33 additional multiplex families with CMM from several continents provided further evidence for linkage to the 1p22 region, again strongest in families with the earliest mean age at diagnosis. A nonparametric ordered sequential analysis was used, based on the average age at diagnosis in each family. The highest LOD score, 6.43, was obtained at D1S2779 and occurred when the 15 families with the earliest ages at onset were included. These data provide significant evidence of a novel susceptibility gene for CMM located within chromosome band 1p22.

Zygosity Diagnosis in the Absence of Genotypic Data: An Approach Using Latent Class Analysis

For zygosity diagnosis in the absence of genotypic data, or in the recruitment phase of a twin study where only single twins from same-sex pairs are being screened, or to provide a test for sample duplication leading to the false identification of a dizygotic pair as monozygotic, the appropriate analysis of respondents' answers to questions about zygosity is critical. Using data from a young adult Australian twin cohort (N = 2094 complete pairs and 519 singleton twins from same-sex pairs with complete responses to all zygosity items), we show that application of latent class analysis (LCA), fitting a 2-class model, yields results that show good concordance with traditional methods of zygosity diagnosis, but with certain important advantages. These include the ability, in many cases, to assign zygosity with specified probability on the basis of responses of a single informant (advantageous when one zygosity type is being oversampled); and the ability to quantify the probability of misassignment of zygosity, allowing prioritization of cases for genotyping as well as identification of cases of probable laboratory error. Out of 242 twins (from 121 like-sex pairs) where genotypic data were available for zygosity confirmation, only a single case was identified of incorrect zygosity assignment by the latent class algorithm. Zygosity assignment for that single case was identified by the LCA as uncertain (probability of being a monozygotic twin only 76%), and the co-twin's responses clearly identified the pair as dizygotic (probability of being dizygotic 100%). In the absence of genotypic data, or as a safeguard against sample duplication, application of LCA for zygosity assignment or confirmation is strongly recommended.

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.

Multiple ordered phases in Langmuir-Blodgett films of cadmium arachidate at elevated temperatures

Using synchrotron X-ray grazing incidence diffraction, superlattice structures have been observed to develop in Langmuir-Blodgett films of cadmium arachidate as the temperature is raised. The previously reported superstructure in the stacked lamellae at room temperature changes at about 70 degreesC and there are further changes at about 90 and 103 degreesC before the major phase transition from stacked lamellae to hexagonally packed rods occurs at 107 degreesC (Langmuir 1997, 13, 1602). Between 70 and 103 degreesC there is a 1 x 10 one-dimensional in-plane superstructure, which is commensurate with the local structure and has an interlayer shift along [01] by a distance of b (of the local structure) at lower temperatures, and a further shift at about 90 degreesC. At lower (

Data Mart para apoio ?? capacita????o em gest??o de pol??ticas p??blicas voltado para uma Escola de Administra????o P??blica

Este trabalho consiste num estudo de caso que se destina ao desenvolvimento de um Data Mart que possibilite a Escola Nacional de Administra????o P??blica ??? ENAP conhecer o perfil e o panorama geral da situa????o funcional dos servidores p??blicos federais que se capacitaram na Escola nos ??ltimos 7 anos. O aplicativo foi desenvolvido cruzando o banco de dados do sistema gerenciador dos cursos ministrados pela ENAP, onde est??o armazenadas informa????es sobre os alunos capacitados, os cursos realizados, os resultados alcan??ados, o perfil dos docentes e demais informa????es relativas ??s atividades da Escola, com os dados gerados pelo Sistema Integrado de Administra????o de Recursos Humanos ??? SIAPE, cuja extra????o de dados foi direcionada para os registros sobre a situa????o funcional, cargos, carreiras, fun????es, ??rg??os e alguns dados pessoais dos alunos, servidores p??blicos federais que se encontram registrados no SIAPE

New species of Bakeriella (Hymenoptera, Bethylidae) from the New World and new data on geographic distribution

Bakeriella lata sp. nov. (Brazil, Rondônia), Bakeriella aurata sp. nov. (Brazil, Amazonas) and Bakeriella sulcaticeps sp. nov. (Brazil, Amazonas) are described and illustrated. New geographic records and variation data for B. cristata Evans, 1964, B. floridana Evans, 1964, B. flavicornis Kieffer, 1910, B. incompleta Azevedo, 1994, B. mira Evans, 1997, B. montivaga (Kieffer, 1910), B. olmeca Evans, 1964 and B. subcarinata Evans, 1965 are provided. The male of B. incompleta is described for the first time.

“Somos ensinados a pensar em sexo”: Representações Sociais de Masculinidades e de Amor em Travestis Homens Gays e Homens Heterossexuais

Estudos realizados no campo das masculinidades mostraram que a adesão e identificação de determinados grupos a um conjunto de regras e comportamentos que definem o ser homem implicava em adoecimento e submissão a um regime que supostamente lhes oferece vantagens e poder. Compreender como diferentes grupos constituem-se em relação à divisão social sexual torna-se relevante principalmente se considerarmos sujeitos que fogem ao padrão que define modos de ser homem e ser mulher. Esta tese objetivou investigar as representações sociais de masculinidades e de amor de sujeitos com identidades de gênero e orientações sexuais diversas, especificamente: identificar e analisar representações sociais de masculinidade e amor de travestis, homens gays e homens heterossexuais; apreender experiências de preconceito e discriminação vividas em função da identidade de gênero e orientações sexuais. A pesquisa foi realizada em duas etapas: entrevistas semiestruturadas com 21 travestis; aplicação de questionários com questões abertas e fechadas em 52 homens gays, 40 homens heterossexuais e 39 travestis. Os dados coletados na primeira etapa foram submetidos aos procedimentos da análise de conteúdo categorial temática. O tratamento dos dados coletados na segunda etapa foi realizado através do software ALCESTE. Os homens heterossexuais representam a masculinidade a partir de uma perspectiva evolutiva, ancorada no discurso biológico e científico que descarta as influências da cultura na constituição dos sexos e identidades de gênero. Nos homens gays encontramos aproximações das representações dos sujeitos ao modelo da masculinidade hegemônica fortemente identificada à imagem do homem heterossexual, bem como elementos que ora apresentam aproximações, ora distanciamentos desse padrão. As travestis também utilizam elementos da ideia de masculinidade hegemônica para construção desta representação, objetivando-a na construção do modelo do homem forte e viril, sobre o qual não se identificam mas direcionam seu desejo. Os três grupos representam o amor como elemento estranho à masculinidade, posto que este é objeto que pertence à feminilidade. O amor distancia-se da construção da representação social de masculinidade de homens heterossexuais, é elemento estranho e conflituoso às representações dos homens gays e valor que corresponde a uma ética para o grupo de travestis. Nesse contexto, as experiências de preconceito e discriminação pouco modificam a elaboração das representações de heterossexuais; influenciam fortemente a construção da imagem de homem e mulher por onde homens gays e travestis alimentam suas identificações. Concluímos apostando na defesa da emergência de práticas e identidades sexo-diversas como artifícios potentes à desestabilização do padrão da dominação masculina para determinação de sexos e identidades.