Parametric analysis of discrepant sets of data

In this paper a new parametric method to deal with discrepant experimental results is developed. The method is based on the fit of a probability density function to the data. This paper also compares the characteristics of different methods used to deduce recommended values and uncertainties from a discrepant set of experimental data. The methods are applied to the (137)Cs and (90)Sr published half-lives and special emphasis is given to the deduced confidence intervals. The obtained results are analyzed considering two fundamental properties expected from an experimental result: the probability content of confidence intervals and the statistical consistency between different recommended values. The recommended values and uncertainties for the (137)Cs and (90)Sr half-lives are 10,984 (24) days and 10,523 (70) days, respectively. (C) 2009 Elsevier B.V. All rights reserved.

Assessing food dietary intakes in Japanese-Brazilians using factor analysis

We describe the use of factor analysis for assessing food habits in Japanese-Brazilians. Dietary data from 1,283 participants of a cross-sectional study were used. Besides statistical criteria, we also used the conceptual meaning of identified profiles to obtain scores for dietary patterns (Japanese or Western profile). Paired Student t test, linear regression and Poisson models were used to verify the existence of relationship between these scores and generation, body mass index (BMI), waist circumference and presence of metabolic syndrome, respectively. First generation subjects had higher mean Japanese profile scores and lower Western profile scores than those of second generation. The Western dietary pattern was associated with BMI (p = 0.001), waist circumference (p = 0.023) and metabolic syndrome (p < 0.05). We concluded that these scores were able to discriminate subjects who maintained their traditional Japanese lifestyle or otherwise, and that the incorporation of a Western lifestyle is associated to high values of BMI, waist circumference and presence of metabolic syndrome.

Trends in dermatomyositis- and polymyositis-related mortality in the state of Sao Paulo, Brazil, 1985-2007: multiple cause-of-death analysis

Background: Dermatomyositis (DM) and polymyositis (PM) are rare systemic autoimmune rheumatic diseases with high fatality rates. There have been few population-based mortality studies of dermatomyositis and polymyositis in the world, and none have been conducted in Brazil. The objective of the present study was to employ multiple-cause of-death methodology in the analysis of trends in mortality related to dermatomyositis and polymyositis in the state of Sao Paulo, Brazil, between 1985 and 2007. Methods: We analyzed mortality data from the Sao Paulo State Data Analysis System, selecting all death certificates on which DM or PM was listed as a cause of death. The variables sex, age and underlying, associated or total mentions of causes of death were studied using mortality rates, proportions and historical trends. Statistical analysis were performed by chi-square and H Kruskal-Wallis tests, variance analysis and linear regression. A p value less than 0.05 was regarded as significant. Results: Over a 23-year period, there were 318 DM-related deaths and 316 PM-related deaths. Overall, DM/PM was designated as an underlying cause in 55.2% and as an associated cause in 44.8%; among 634 total deaths females accounted for 71.5%. During the study period, age-and gender-adjusted DM mortality rates did not change significantly, although PM as an underlying cause and total mentions of PM trended lower (p < 0.05). The mean ages at death were 47.76 +/- 20.81 years for DM and 54.24 +/- 17.94 years for PM (p = 0.0003). For DM/PM, respectively, as underlying causes, the principal associated causes of death were as follows: pneumonia (in 43.8%/33.5%); respiratory failure (in 34.4%/32.3%); interstitial pulmonary diseases and other pulmonary conditions (in 28.9%/17.6%); and septicemia (in 22.8%/15.9%). For DM/PM, respectively, as associated causes, the following were the principal underlying causes of death: respiratory disorders (in 28.3%/26.0%); circulatory disorders (in 17.4%/20.5%); neoplasms (in 16.7%/13.7%); infectious and parasitic diseases (in 11.6%/9.6%); and gastrointestinal disorders (in 8.0%/4.8%). Of the 318 DM-related deaths, 36 involved neoplasms, compared with 20 of the 316 PM-related deaths (p = 0.03). Conclusions: Our study using multiple cause of deaths found that DM/PM were identified as the underlying cause of death in only 55.2% of the deaths, indicating that both diseases were underestimated in the primary mortality statistics. We observed a predominance of deaths in women and in older individuals, as well as a trend toward stability in the mortality rates. We have confirmed that the risk of death is greater when either disease is accompanied by neoplasm, albeit to lesser degree in individuals with PM. The investigation of the underlying and associated causes of death related to DM/PM broaden the knowledge of the natural history of both diseases and could help integrate mortality data for use in the evaluation of control measures for DM/PM.

Log-modified Weibull regression models with censored data: Sensitivity and residual analysis

This paper proposes a regression model considering the modified Weibull distribution. This distribution can be used to model bathtub-shaped failure rate functions. Assuming censored data, we consider maximum likelihood and Jackknife estimators for the parameters of the model. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and we also present some ways to perform global influence. Besides, for different parameter settings, sample sizes and censoring percentages, various simulations are performed and the empirical distribution of the modified deviance residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for a martingale-type residual in log-modified Weibull regression models with censored data. Finally, we analyze a real data set under log-modified Weibull regression models. A diagnostic analysis and a model checking based on the modified deviance residual are performed to select appropriate models. (c) 2008 Elsevier B.V. All rights reserved.

Regression models for grouped survival data: Estimation and sensitivity analysis

In this study, regression models are evaluated for grouped survival data when the effect of censoring time is considered in the model and the regression structure is modeled through four link functions. The methodology for grouped survival data is based on life tables, and the times are grouped in k intervals so that ties are eliminated. Thus, the data modeling is performed by considering the discrete models of lifetime regression. The model parameters are estimated by using the maximum likelihood and jackknife methods. To detect influential observations in the proposed models, diagnostic measures based on case deletion, which are denominated global influence, and influence measures based on small perturbations in the data or in the model, referred to as local influence, are used. In addition to those measures, the local influence and the total influential estimate are also employed. Various simulation studies are performed and compared to the performance of the four link functions of the regression models for grouped survival data for different parameter settings, sample sizes and numbers of intervals. Finally, a data set is analyzed by using the proposed regression models. (C) 2010 Elsevier B.V. All rights reserved.

The exponentiated generalized gamma distribution with application to lifetime data

A four-parameter extension of the generalized gamma distribution capable of modelling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma and generalized Rayleigh, among others. We derive two infinite sum representations for its moments. We calculate the density of the order statistics and two expansions for their moments. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is obtained. Finally, a real data set from the medical area is analysed.

PCA Tomography: how to extract information from data cubes

Astronomy has evolved almost exclusively by the use of spectroscopic and imaging techniques, operated separately. With the development of modern technologies, it is possible to obtain data cubes in which one combines both techniques simultaneously, producing images with spectral resolution. To extract information from them can be quite complex, and hence the development of new methods of data analysis is desirable. We present a method of analysis of data cube (data from single field observations, containing two spatial and one spectral dimension) that uses Principal Component Analysis (PCA) to express the data in the form of reduced dimensionality, facilitating efficient information extraction from very large data sets. PCA transforms the system of correlated coordinates into a system of uncorrelated coordinates ordered by principal components of decreasing variance. The new coordinates are referred to as eigenvectors, and the projections of the data on to these coordinates produce images we will call tomograms. The association of the tomograms (images) to eigenvectors (spectra) is important for the interpretation of both. The eigenvectors are mutually orthogonal, and this information is fundamental for their handling and interpretation. When the data cube shows objects that present uncorrelated physical phenomena, the eigenvector`s orthogonality may be instrumental in separating and identifying them. By handling eigenvectors and tomograms, one can enhance features, extract noise, compress data, extract spectra, etc. We applied the method, for illustration purpose only, to the central region of the low ionization nuclear emission region (LINER) galaxy NGC 4736, and demonstrate that it has a type 1 active nucleus, not known before. Furthermore, we show that it is displaced from the centre of its stellar bulge.

Comparing diagnostic tests with missing data

When missing data occur in studies designed to compare the accuracy of diagnostic tests, a common, though naive, practice is to base the comparison of sensitivity, specificity, as well as of positive and negative predictive values on some subset of the data that fits into methods implemented in standard statistical packages. Such methods are usually valid only under the strong missing completely at random (MCAR) assumption and may generate biased and less precise estimates. We review some models that use the dependence structure of the completely observed cases to incorporate the information of the partially categorized observations into the analysis and show how they may be fitted via a two-stage hybrid process involving maximum likelihood in the first stage and weighted least squares in the second. We indicate how computational subroutines written in R may be used to fit the proposed models and illustrate the different analysis strategies with observational data collected to compare the accuracy of three distinct non-invasive diagnostic methods for endometriosis. The results indicate that even when the MCAR assumption is plausible, the naive partial analyses should be avoided.

Consumo alimentar e ecologia de populações ribeirinhas em dois ecossistemas amazônicos: um estudo comparativo

OBJETIVO: Este artigo analisa e compara os dados de consumo alimentar de duas populações ribeirinhas da Amazônia vivendo em ecossistemas contrastantes de floresta tropical: a várzea estacional e a floresta de terra firme. MÉTODOS: Foi estudado o consumo alimentar de 11 unidades domésticas na várzea (Ilha de Ituqui, Município de Santarém) e 17 na terra firme (Floresta Nacional de Caxiuanã, Municípios de Melgaço e Portel). O método utilizado foi o recordatório de 24 horas. As análises estatísticas foram executadas com o auxílio do programa Statistical Package for Social Sciences 12.0. RESULTADOS: Em ambos os ecossistemas, os resultados confirmam a centralidade do pescado e da mandioca na dieta local. Porém, a contribuição de outros itens alimentares secundários, tais como o açaí (em Caxiuanã) e o leite in natura (em Ituqui), também foi significante. Além disso, o açúcar revelou ser uma fonte de energia confiável para enfrentar as flutuações sazonais dos recursos naturais. Parece haver ainda uma maior contribuição energética dos peixes para a dieta de Ituqui, provavelmente em função da maior produtividade dos rios e lagos da várzea em relação à terra firme. Por fim, Ituqui revelou uma maior dependência de itens alimentares comprados, enquanto Caxiuanã mostrou estar ainda bastante vinculada à agricultura e às redes locais de troca. CONCLUSÃO: Além dos resultados confirmarem a importância do pescado e da mandioca, também mostraram que produtos industrializados, como o açúcar, têm um papel importante nas dietas, podendo apontar para tendências no consumo alimentar relacionadas com a atual transição nutricional e com a erosão, em diferentes níveis, dos sistemas de subsistência locais.

Pre-processing for noise detection in gene expression classification data

Due to the imprecise nature of biological experiments, biological data is often characterized by the presence of redundant and noisy data. This may be due to errors that occurred during data collection, such as contaminations in laboratorial samples. It is the case of gene expression data, where the equipments and tools currently used frequently produce noisy biological data. Machine Learning algorithms have been successfully used in gene expression data analysis. Although many Machine Learning algorithms can deal with noise, detecting and removing noisy instances from the training data set can help the induction of the target hypothesis. This paper evaluates the use of distance-based pre-processing techniques for noise detection in gene expression data classification problems. This evaluation analyzes the effectiveness of the techniques investigated in removing noisy data, measured by the accuracy obtained by different Machine Learning classifiers over the pre-processed data.

Variância intrapessoal para ajuste da distribuição de nutrientes em estudos epidemiológicos

O objetivo do estudo foi apresentar a fração da variância intrapessoal para ajuste da distribuição de nutrientes de adultos e idosos. Utilizaram-se dados de inquérito populacional com amostra representativa (n = 511) de indivíduos com 19 anos ou mais do município de São Paulo, SP, em 2007. A fração da variância intrapessoal foi obtida pelo método proposto pela Iowa State University. Observaram-se diferenças nas frações das variâncias intrapessoais de nutrientes segundo sexo. Esses valores devem ser utilizados para ajustar a distribuição da ingestão de nutrientes, pois sua não utilização pode resultar em viés na análise e interpretação de dados.

Avaliação da ingestão de alimentos por nipo-brasileiros utilizando análise fatorial

Descreve-se o uso da análise fatorial na avaliação dos hábitos alimentares de nipo-brasileiros. Utilizaram-se dados dietéticos de 1.283 participantes de estudo transversal. A partir de critérios estatísticos e do significado conceitual dos padrões identificados, geraram-se escores que definiram os perfis dietéticos (japonês ou ocidental). Empregou-se o teste t de Student pareado, os modelos de regressão linear e de Poisson para examinar as relações desses escores com, respectivamente, a geração, índice de massa corporal (IMC), perímetro abdominal e a presença de síndrome metabólica. Aqueles de primeira geração, em relação aos de segunda, apresentaram escores maiores para o perfil japonês e menores para o ocidental. O perfil ocidental relacionou-se com o IMC (p = 0,001), perímetro abdominal (p = 0,023) e a síndrome metabólica (p < 0,05). Conclui-se que os escores discriminaram sujeitos que mantêm ou não estilo de vida tradicional japonês e que a incorporação de hábitos ocidentais associou-se a maiores valores de IMC, perímetro abdominal e a presença de síndrome metabólica

Kinematic parameters and membership probabilities of open clusters in the Bordeaux PM2000 catalogue

Aims. We derive lists of proper-motions and kinematic membership probabilities for 49 open clusters and possible open clusters in the zone of the Bordeaux PM2000 proper motion catalogue (+ 11 degrees <= delta <= + 18 degrees). We test different parametrisations of the proper motion and position distribution functions and select the most successful one. In the light of those results, we analyse some objects individually. Methods. We differenciate between cluster and field member stars, and assign membership probabilities, by applying a new and fully automated method based on both parametrisations of the proper motion and position distribution functions, and genetic algorithm optimization heuristics associated with a derivative-based hill climbing algorithm for the likelihood optimization. Results. We present a catalogue comprising kinematic parameters and associated membership probability lists for 49 open clusters and possible open clusters in the Bordeaux PM2000 catalogue region. We note that this is the first determination of proper motions for five open clusters. We confirm the non-existence of two kinematic populations in the region of 15 previously suspected non-existent objects.

Automated supervised classification of variable stars in the CoRoT programme Method and application to the first four exoplanet fields

Aims. In this work, we describe the pipeline for the fast supervised classification of light curves observed by the CoRoT exoplanet CCDs. We present the classification results obtained for the first four measured fields, which represent a one-year in-orbit operation. Methods. The basis of the adopted supervised classification methodology has been described in detail in a previous paper, as is its application to the OGLE database. Here, we present the modifications of the algorithms and of the training set to optimize the performance when applied to the CoRoT data. Results. Classification results are presented for the observed fields IRa01, SRc01, LRc01, and LRa01 of the CoRoT mission. Statistics on the number of variables and the number of objects per class are given and typical light curves of high-probability candidates are shown. We also report on new stellar variability types discovered in the CoRoT data. The full classification results are publicly available.

Recurrence quantification analysis of turbulent fluctuations in the plasma edge of Tokamak Chauffage Alfveacuten Breacutesilien tokamak

Recurrences are close returns of a given state in a time series, and can be used to identify different dynamical regimes and other related phenomena, being particularly suited for analyzing experimental data. In this work, we use recurrence quantification analysis to investigate dynamical patterns in scalar data series obtained from measurements of floating potential and ion saturation current at the plasma edge of the Tokamak Chauffage Alfveacuten Breacutesilien [R. M. O. Galva approximate to o , Plasma Phys. Controlled Fusion 43, 1181 (2001)]. We consider plasma discharges with and without the application of radial electric bias, and also with two different regimes of current ramp. Our results indicate that biasing improves confinement through destroying highly recurrent regions within the plasma column that enhance particle and heat transport.