Decomposition of $\rm R\times C$ contingency table chi -square: Application to binned DNA fragment size data and population structure analysis

The purpose of this research is to develop a new statistical method to determine the minimum set of rows (R) in a R x C contingency table of discrete data that explains the dependence of observations. The statistical power of the method will be empirically determined by computer simulation to judge its efficiency over the presently existing methods. The method will be applied to data on DNA fragment length variation at six VNTR loci in over 72 populations from five major racial groups of human (total sample size is over 15,000 individuals; each sample having at least 50 individuals). DNA fragment lengths grouped in bins will form the basis of studying inter-population DNA variation within the racial groups are significant, will provide a rigorous re-binning procedure for forensic computation of DNA profile frequencies that takes into account intra-racial DNA variation among populations. ^

A unified approach for representing rows and columns in contingency tables

By using suitable parameters, we present a uni¯ed aproach for describing four methods for representing categorical data in a contingency table. These methods include:correspondence analysis (CA), the alternative approach using Hellinger distance (HD),the log-ratio (LR) alternative, which is appropriate for compositional data, and theso-called non-symmetrical correspondence analysis (NSCA). We then make an appropriate comparison among these four methods and some illustrative examples are given.Some approaches based on cumulative frequencies are also linked and studied usingmatrices.Key words: Correspondence analysis, Hellinger distance, Non-symmetrical correspondence analysis, log-ratio analysis, Taguchi inertia

A unified approach for representing rows and columns in contingency tables

By using suitable parameters, we present a uni¯ed aproach for describing four methods for representing categorical data in a contingency table. These methods include: correspondence analysis (CA), the alternative approach using Hellinger distance (HD), the log-ratio (LR) alternative, which is appropriate for compositional data, and the so-called non-symmetrical correspondence analysis (NSCA). We then make an appropriate comparison among these four methods and some illustrative examples are given. Some approaches based on cumulative frequencies are also linked and studied using matrices. Key words: Correspondence analysis, Hellinger distance, Non-symmetrical correspondence analysis, log-ratio analysis, Taguchi inertia

Statnote 7: chi-square contingency tables

When the data are counts or the frequencies of particular events and can be expressed as a contingency table, then they can be analysed using the chi-square distribution. When applied to a 2 x 2 table, the test is approximate and care needs to be taken in analysing tables when the expected frequencies are small either by applying Yate’s correction or by using Fisher’s exact test. Larger contingency tables can also be analysed using this method. Note that it is a serious statistical error to use any of these tests on measurement data!

Thermostability of heterophile antibodies from human sera infected with Schistosoma mansoni and geo-helminths: an immuno-metric statistical analysis

Antibody in human sera that induces lysis of sheep erythrocytes in hemolytic assay was investigated. The present study showed that the presence in serum of the thermostable cytolytic anti-sheep red blood cells antibodies is dependent on the Schistosoma mansoni infection, and this is more frequent in adults than in children. The thermostable characteristic of hemolysins in normal sera was not dependent on the presence of Ascaris lumbricoides, Trichuris trichiura or hookworm geo-helminths. Further, thermostable complement-activating heterophile antibodies were noticed in children in association with massive number of S. mansoni eggs. The results were obtained by using the z- and the chi-square tests. The z-test allows us to formulate a one-sided alternative, i.e., a tendency of one of the attributes. On the other hand, the chi-square test analyzes the independence between attributes by using a contingency table. Besides the obtained results being interesting in the field of schistosomiasis mansoni, they can provide a new insight into the use of statistics in medical science.

One Anastomosis Gastric Bypass versus Roux en Y Gastric Bypass as Salvage Technique after Failed Gastric Band: a Retrospective Analysis of 80 Cases

The classical “Roux en Y Gastric Bypass” (RYGBP) is still the standard technique between all the ones being used nowadays. The “One anastomosis gastric bypass” (OAGBP), is an evolution of the “Minigastic bypass” described by Robert Rutledge in 2001, is a well known and progressively frequent but still controversial technique. In our group, after an experience of 10 years using the RYGBP as a salvage surgery after failed gastric banding, in 132 cases, we decided to adopt the OAGBP as our preferential bariatric technique also in this situation. The theoretical main reasons for that shift are related to the increased safety, maximized weight loss, long term weight loss maintenance and reversibility of the operation. Method: Retrospectively we evaluated data of the surgical management of revisional cases for conversion, after failed or complicated gastric bands to gastric bypass. We selected the last 40 cases of each technique since May 2010. Results: All cases were performed by laparoscopy without any conversion. In both groups the conversion has been performed in one single step (17 cases, 42,5%). Data showed lower morbidity with OAGBP (2,5% against 7,5%) and better weight loss in theOAGBP cohort after a median follow up of 16months (67%against 55%) in patients revised after gastric band failure or complications. None had statistic significance (p>0,1) by the chi-square contingency table analysis.Conclusion: It seems to there is a difference in favour of OAGBP for conversion of complicated gastric bands. In this study we didn’t found statistic significance probably because of the short numbers. Prospective and more powerful studies are necessary to evaluate the benefit of the studied procedure.

Bancroftian filariasis in an endemic area of Brazil: differences between genders during puberty

Gender differences in susceptibility to infectious diseases have been observed in various studies. A survey was performed in a bancroftian filariasis endemic area in the city of Olinda, Brazil. All residents aged 5 years or older were examined by thick blood film. People aged 9 to 16 years were interviewed and also tested for filarial antigenaemia. Data were analyzed by contingency table methods and regression models. The risk of microfilaraemia for males was significantly higher. Among those aged 9 to 16 years, the analysis of gender and filariasis by age showed that boys from 15 to 16 years had a higher risk of infection than girls. No association was found between menarche and filariasis in girls. The data suggest that variations between gender in filariasis could result, at least in part, from an increase in susceptibility of men. This epidemiologic feature needs to be considered while formulating elimination plans.

Comparing correspondence analysis and the log-ratio alternative for representing categorical data

We compare correspondance análisis to the logratio approach based on compositional data. We also compare correspondance análisis and an alternative approach using Hellinger distance, for representing categorical data in a contingency table. We propose a coefficient which globally measures the similarity between these approaches. This coefficient can be decomposed into several components, one component for each principal dimension, indicating the contribution of the dimensions to the difference between the two representations. These three methods of representation can produce quite similar results. One illustrative example is given

The prevalence of schistosomiasis in school-aged children as an appropriate indicator of its prevalence in the community

School-aged children (6-15 years) from the endemic area of Pernambuco were evaluated both as a target group for and an indicator of schistosomiasis control in the community. Parasitological data were drawn from baseline stool surveys of whole populations that were obtained to diagnose Schistosoma mansoni infection. Nineteen representative localities were selected for assessing the prevalence of schistosomiasis among individuals in the following age groups: 0-5, 6-15, 16-25, 26-40 and 41-80 years. For each locality, the prevalence in each age group was compared to that of the overall population using contingency table analysis. To select a reference group, the operational difficulties of conducting residential surveys were considered. School-aged children may be considered to be the group of choice as the reference group for the overall population for the following reasons: (i) the prevalence of schistosomiasis in this age group had the highest correlation with the prevalence in the overall population (r = 0.967), (ii) this age group is particularly vulnerable to infection and plays an important role in parasite transmission and (iii) school-aged children are the main target of the World Health Organization in terms of helminth control. The Schistosomiasis Control Program should consider school-aged children both as a reference group for assessing the need for intervention at the community level and as a target group for integrated health care actions of the Unified Health System that are focused on high-risk groups.

Compositional analysis of bivariate discrete probabilities

A joint distribution of two discrete random variables with finite support can be displayed as a two way table of probabilities adding to one. Assume that this table hasn rows and m columns and all probabilities are non-null. This kind of table can beseen as an element in the simplex of n · m parts. In this context, the marginals areidentified as compositional amalgams, conditionals (rows or columns) as subcompositions. Also, simplicial perturbation appears as Bayes theorem. However, the Euclideanelements of the Aitchison geometry of the simplex can also be translated into the tableof probabilities: subspaces, orthogonal projections, distances.Two important questions are addressed: a) given a table of probabilities, which isthe nearest independent table to the initial one? b) which is the largest orthogonalprojection of a row onto a column? or, equivalently, which is the information in arow explained by a column, thus explaining the interaction? To answer these questionsthree orthogonal decompositions are presented: (1) by columns and a row-wise geometric marginal, (2) by rows and a columnwise geometric marginal, (3) by independenttwo-way tables and fully dependent tables representing row-column interaction. Animportant result is that the nearest independent table is the product of the two (rowand column)-wise geometric marginal tables. A corollary is that, in an independenttable, the geometric marginals conform with the traditional (arithmetic) marginals.These decompositions can be compared with standard log-linear models.Key words: balance, compositional data, simplex, Aitchison geometry, composition,orthonormal basis, arithmetic and geometric marginals, amalgam, dependence measure,contingency table

Ratio maps and correspondence analysis

We compare two methods for visualising contingency tables and developa method called the ratio map which combines the good properties of both.The first is a biplot based on the logratio approach to compositional dataanalysis. This approach is founded on the principle of subcompositionalcoherence, which assures that results are invariant to considering subsetsof the composition. The second approach, correspondence analysis, isbased on the chi-square approach to contingency table analysis. Acornerstone of correspondence analysis is the principle of distributionalequivalence, which assures invariance in the results when rows or columnswith identical conditional proportions are merged. Both methods may bedescribed as singular value decompositions of appropriately transformedmatrices. Correspondence analysis includes a weighting of the rows andcolumns proportional to the margins of the table. If this idea of row andcolumn weights is introduced into the logratio biplot, we obtain a methodwhich obeys both principles of subcompositional coherence and distributionalequivalence.

Análise de Dados de Desemprego em Cabo Verde

Este trabalho baseia-se na análise de dados do desemprego em Cabo Verde nos anos de 2006 e 2008, usando informação da base de dados do INE e IEFP. Partindo da análise dos dados em estudo vai-se procurar descrever e perspectivar metodologias que contemplam as variáveis qualitativas e quantitativas com significado social positivo para a sociedade deste país. Após a introdução no capítulo 1, fez-se, no capítulo 2, a análise exploratória dos dados do desemprego em Cabo Verde referente aos anos 2006 e 2008. No capítulo 3 estudam-se associações entre variáveis, usando a metodologia de tabelas contingência, através da realização de testes de independência e testes de homogeneidade, e análise de medidas de associação. As variáveis usadas, vão ser essencialmente, o escalão etário, o género e o ano. O capítulo 4 é dedicado ao estudo de modelos Log - lineares em tabela de contingência, finalizando-se o trabalho com a apresentação das principais conclusões.

Tying up the loose ends in simple correspondence analysis

Although correspondence analysis is now widely available in statistical software packages and applied in a variety of contexts, notably the social and environmental sciences, there are still some misconceptions about this method as well as unresolved issues which remain controversial to this day. In this paper we hope to settle these matters, namely (i) the way CA measures variance in a two-way table and how to compare variances between tables of different sizes, (ii) the influence, or rather lack of influence, of outliers in the usual CA maps, (iii) the scaling issue and the biplot interpretation of maps,(iv) whether or not to rotate a solution, and (v) statistical significance of results.

A limited area model intercomparison on the "Montserrat-2000" flash-flood event using statistical and deterministic methods

In the scope of the European project Hydroptimet, INTERREG IIIB-MEDOCC programme, limited area model (LAM) intercomparison of intense events that produced many damages to people and territory is performed. As the comparison is limited to single case studies, the work is not meant to provide a measure of the different models' skill, but to identify the key model factors useful to give a good forecast on such a kind of meteorological phenomena. This work focuses on the Spanish flash-flood event, also known as "Montserrat-2000" event. The study is performed using forecast data from seven operational LAMs, placed at partners' disposal via the Hydroptimet ftp site, and observed data from Catalonia rain gauge network. To improve the event analysis, satellite rainfall estimates have been also considered. For statistical evaluation of quantitative precipitation forecasts (QPFs), several non-parametric skill scores based on contingency tables have been used. Furthermore, for each model run it has been possible to identify Catalonia regions affected by misses and false alarms using contingency table elements. Moreover, the standard "eyeball" analysis of forecast and observed precipitation fields has been supported by the use of a state-of-the-art diagnostic method, the contiguous rain area (CRA) analysis. This method allows to quantify the spatial shift forecast error and to identify the error sources that affected each model forecasts. High-resolution modelling and domain size seem to have a key role for providing a skillful forecast. Further work is needed to support this statement, including verification using a wider observational data set.

Genes involved in the metabolism of fatty acids and risk for Crohn's disease in children: a candidate gene study

Contexte - La prévalence de la maladie de Crohn (MC), une maladie inflammatoire chronique du tube digestif, chez les enfants canadiens se situe parmi les plus élevées au monde. Les interactions entre les réponses immunes innées et acquises aux microbes de l'hôte pourraient être à la base de la transition de l’inflammation physiologique à une inflammation pathologique. Le leucotriène B4 (LTB4) est un modulateur clé de l'inflammation et a été associé à la MC. Nous avons postulé que les principaux gènes impliqués dans la voie métabolique du LTB4 pourrait conférer une susceptibilité accrue à l'apparition précoce de la MC. Dans cette étude, nous avons exploré les associations potentielles entre les variantes de l'ADN des gènes ALOX5 et CYP4F2 et la survenue précoce de la MC. Nous avons également examiné si les gènes sélectionnés montraient des effets parent-d'origine, influençaient les phénotypes cliniques de la MC et s'il existait des interactions gène-gène qui modifieraient la susceptibilité à développer la MC chez l’enfant. Méthodes – Dans le cadre d’une étude de cas-parents et de cas-témoins, des cas confirmés, leurs parents et des contrôles ont été recrutés à partir de trois cliniques de gastro-entérologie à travers le Canada. Les associations entre les polymorphismes de remplacement d'un nucléotide simple (SNP) dans les gènes CYP4F2 et ALOX5 ont été examinées. Les associations allélique et génotypiques ont été examinées à partir d’une analyse du génotype conditionnel à la parenté (CPG) pour le résultats cas-parents et à l’aide de table de contingence et de régression logistique pour les données de cas-contrôles. Les interactions gène-gène ont été explorées à l'aide de méthodes de réduction multi-factorielles de dimensionnalité (MDR). Résultats – L’étude de cas-parents a été menée sur 160 trios. L’analyse CPG pour 14 tag-SNP (10 dans la CYP4F2 et 4 dans le gène ALOX5) a révélé la présence d’associations alléliques ou génotypique significatives entre 3 tag-SNP dans le gène CYP4F2 (rs1272, p = 0,04, rs3093158, p = 0.00003, et rs3093145, p = 0,02). Aucune association avec les SNPs de ALOX5 n’a pu être démontrée. L’analyse de l’haplotype de CYP4F2 a montré d'importantes associations avec la MC (test omnibus p = 0,035). Deux haplotypes (GAGTTCGTAA, p = 0,05; GGCCTCGTCG, p = 0,001) montraient des signes d'association avec la MC. Aucun effet parent-d'origine n’a été observé. Les tentatives de réplication pour trois SNPs du gene CYP4F2 dans l'étude cas-témoins comportant 225 cas de MC et 330 contrôles suggèrent l’association dans un de ceux-ci (rs3093158, valeur non-corrigée de p du test unilatéral = 0,03 ; valeur corrigée de p = 0.09). La combinaison des ces deux études a révélé des interactions significatives entre les gènes CYP4F2, ALOX et NOD2. Nous n’avons pu mettre en évidence aucune interaction gène-sexe, de même qu’aucun gène associé aux phénotypes cliniques de la MC n’a pu être identifié. Conclusions - Notre étude suggère que la CYP4F2, un membre clé de la voie métabolique LTB4 est un gène candidat potentiel pour MC. Nous avons également pu mettre en évidence que les interactions entre les gènes de l'immunité adaptative (CYP4F2 et ALOX5) et les gènes de l'immunité innée (NOD2) modifient les risques de MC chez les enfants. D'autres études sur des cohortes plus importantes sont nécessaires pour confirmer ces conclusions.