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

CODAPACK3D. A new version of Compositional Data Package

The statistical analysis of compositional data should be treated using logratios of parts, which are difficult to use correctly in standard statistical packages. For this reason a freeware package, named CoDaPack was created. This software implements most of the basic statistical methods suitable for compositional data. In this paper we describe the new version of the package that now is called CoDaPack3D. It is developed in Visual Basic for applications (associated with Excel©), Visual Basic and Open GL, and it is oriented towards users with a minimum knowledge of computers with the aim at being simple and easy to use. This new version includes new graphical output in 2D and 3D. These outputs could be zoomed and, in 3D, rotated. Also a customization menu is included and outputs could be saved in jpeg format. Also this new version includes an interactive help and all dialog windows have been improved in order to facilitate its use. To use CoDaPack one has to access Excel© and introduce the data in a standard spreadsheet. These should be organized as a matrix where Excel© rows correspond to the observations and columns to the parts. The user executes macros that return numerical or graphical results. There are two kinds of numerical results: new variables and descriptive statistics, and both appear on the same sheet. Graphical output appears in independent windows. In the present version there are 8 menus, with a total of 38 submenus which, after some dialogue, directly call the corresponding macro. The dialogues ask the user to input variables and further parameters needed, as well as where to put these results. The web site http://ima.udg.es/CoDaPack contains this freeware package and only Microsoft Excel© under Microsoft Windows© is required to run the software. Kew words: Compositional data Analysis, Software

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 has n rows and m columns and all probabilities are non-null. This kind of table can be seen as an element in the simplex of n · m parts. In this context, the marginals are identified as compositional amalgams, conditionals (rows or columns) as subcompositions. Also, simplicial perturbation appears as Bayes theorem. However, the Euclidean elements of the Aitchison geometry of the simplex can also be translated into the table of probabilities: subspaces, orthogonal projections, distances. Two important questions are addressed: a) given a table of probabilities, which is the nearest independent table to the initial one? b) which is the largest orthogonal projection of a row onto a column? or, equivalently, which is the information in a row explained by a column, thus explaining the interaction? To answer these questions three orthogonal decompositions are presented: (1) by columns and a row-wise geometric marginal, (2) by rows and a columnwise geometric marginal, (3) by independent two-way tables and fully dependent tables representing row-column interaction. An important result is that the nearest independent table is the product of the two (row and column)-wise geometric marginal tables. A corollary is that, in an independent table, 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

Robust Factor Analysis for Compositional Data

Factor analysis as frequent technique for multivariate data inspection is widely used also for compositional data analysis. The usual way is to use a centered logratio (clr) transformation to obtain the random vector y of dimension D. The factor model is then y = Λf + e (1) with the factors f of dimension k < D, the error term e, and the loadings matrix Λ. Using the usual model assumptions (see, e.g., Basilevsky, 1994), the factor analysis model (1) can be written as Cov(y) = ΛΛT + ψ (2) where ψ = Cov(e) has a diagonal form. The diagonal elements of ψ as well as the loadings matrix Λ are estimated from an estimation of Cov(y). Given observed clr transformed data Y as realizations of the random vector y. Outliers or deviations from the idealized model assumptions of factor analysis can severely effect the parameter estimation. As a way out, robust estimation of the covariance matrix of Y will lead to robust estimates of Λ and ψ in (2), see Pison et al. (2003). Well known robust covariance estimators with good statistical properties, like the MCD or the S-estimators (see, e.g. Maronna et al., 2006), rely on a full-rank data matrix Y which is not the case for clr transformed data (see, e.g., Aitchison, 1986). The isometric logratio (ilr) transformation (Egozcue et al., 2003) solves this singularity problem. The data matrix Y is transformed to a matrix Z by using an orthonormal basis of lower dimension. Using the ilr transformed data, a robust covariance matrix C(Z) can be estimated. The result can be back-transformed to the clr space by C(Y ) = V C(Z)V T where the matrix V with orthonormal columns comes from the relation between the clr and the ilr transformation. Now the parameters in the model (2) can be estimated (Basilevsky, 1994) and the results have a direct interpretation since the links to the original variables are still preserved. The above procedure will be applied to data from geochemistry. Our special interest is on comparing the results with those of Reimann et al. (2002) for the Kola project data

Modelització numèrica d’un sistema de descàrrega d’una columna de la catedral d Girona mitjançant el mètode dels elements finits (anàlisi no lineal)

Aquest projecte consisteix en aplicar el càlcul no lineal en la modelització volumètrica numèrica de l’estructura del sistema de descàrrega d’una columna del claustre de la catedral de Girona mitjançant el mètode dels elements finits. A la Universitat de Girona s’ha fet diferents estudis del claustre de la catedral de Girona però sempre simulant un comportament lineal de les característiques dels materials. El programa utilitzat és la versió docent del programa ANSYS disponible al Dept. d’EMCI i l’element emprat ha sigut el SOLID65. Aquest element permet introduir característiques de no linealitat en els models i és adequat per a anàlisi no lineal d’elements com la pedra de Girona

Select-divide-and-conquer method for large-scale configuration interaction

A select-divide-and-conquer variational method to approximate configuration interaction (CI) is presented. Given an orthonormal set made up of occupied orbitals (Hartree-Fock or similar) and suitable correlation orbitals (natural or localized orbitals), a large N-electron target space S is split into subspaces S0,S1,S2,...,SR. S0, of dimension d0, contains all configurations K with attributes (energy contributions, etc.) above thresholds T0={T0egy, T0etc.}; the CI coefficients in S0 remain always free to vary. S1 accommodates KS with attributes above T1≤T0. An eigenproblem of dimension d0+d1 for S0+S 1 is solved first, after which the last d1 rows and columns are contracted into a single row and column, thus freezing the last d1 CI coefficients hereinafter. The process is repeated with successive Sj(j≥2) chosen so that corresponding CI matrices fit random access memory (RAM). Davidson's eigensolver is used R times. The final energy eigenvalue (lowest or excited one) is always above the corresponding exact eigenvalue in S. Threshold values {Tj;j=0, 1, 2,...,R} regulate accuracy; for large-dimensional S, high accuracy requires S 0+S1 to be solved outside RAM. From there on, however, usually a few Davidson iterations in RAM are needed for each step, so that Hamiltonian matrix-element evaluation becomes rate determining. One μhartree accuracy is achieved for an eigenproblem of order 24 × 106, involving 1.2 × 1012 nonzero matrix elements, and 8.4×109 Slater determinants

Proves de càrrega al claustre de la catedral de Girona

El projecte tracta de la realització de proves de càrrega al Claustre de la Catedral de Girona per tal de conèixer la càrrega que suporten les columnes de pedra del claustre. Un cop coneguda aquesta càrrega, es podrà decidir si son aptes o no per al canvi d’ús previst a la planta superior i es podrà aportar una resposta definitiva a una de les dues hipòtesis plantejades: si es forma o no arc de descàrrega a les columnes del Claustre. D’altre banda, es realitzarà una proposta de tractament de les columnes del Claustre que consistirà a escollir les que hauran de ser reparades i les que hauran de ser substituïdes, en base a una sèrie de càlculs numèrics

Aprofitament de residus vegetals per a la concentració i separació d'ions metàl·lics de solucions aquoses

La adsorció sobre carbó actiu és un tractament habitual para a l'eliminació de metalls contaminants d'afluents aquosos. Amb la finalitat de reduir el cost dels tractament actuals s'està estudiant la utilització d'adsorbents de baix cost com a substituts dels sorbents convencionals que s'utilitzen usualment en els tractaments d'aigües. En aquesta tesi s'ha avaluat la possibilitat d'aprofitar diversos residus vegetals procedents de processos agrícoles o residus industrials, com a sorbents per a l'eliminació de metalls de solucions aquoses. En concret, el treball s'ha centrat en assolir dos grans objectius, primer de tot, conèixer la possibilitat d'utilitzar el pinyol d'oliva com a adsorbent de metalls divalents de solucions aquoses, com són el Cd(II), Cu(II), Ni(II) i Pb(II), i en segon lloc, avaluar la utilització de quatre residus vegetals (escorça de yohimbe, rapa de raïm, suro i pinyol d'oliva) per a l'extracció de Cr(VI) de solucions aquoses i utilitzar el residu vegetal més apropiat com a base per tal d'obtenir un nou sorbent, amb més bona morfologia per ser utilitzat en sistemes d'eliminació de Cr(VI) en continu. Les dades cinètiques d'adsorció es van tractar amb el model cinètic de pseudo-segon ordre, amb el que es van determinar la velocitat inicial d'adsorció dels metalls. Els models d'isotermes d'equilibri de Langmuir i Freundlich es van utilitzar per descriure els equilibris d'adsorció i calcular els paràmetres de les isotermes. En els assajos amb mescles binàries de metalls, els resultats es van tractar amb el model ampliat de Langmuir. En aquest treball s'ha determinat, primerament, que el pinyol d'oliva és un bon sorbent dels metalls divalents estudiats. S'ha determinat que l'adsorció dels metalls varia amb el pH de la solució i que la presència de sals en la solució provoca una davallada en l'adsorció. També es va determinar la competència en l'adsorció entre els metalls al realitzar assajos en mescles binàries. Amb els estudis d'extracció de Cr(VI) amb els quatre materials s'ha determinat que aquests vegetals poden adsorbir el Cr(VI) en solució, tot i que la rapa i la yohimbe són els materials que presenten unes capacitats d'adsorció més elevades. Les anàlisis amb espectroscòpia electrònica de Raig-X realitzats en la rapa de raïm i la yohimbe en contacte amb Cr(VI) varen demostrar la presència del metall en els dos estats d'oxidació, Cr(VI) i Cr(III), evidenciant la participació del mecanisme de reducció de Cr(VI) en l'adsorció del metall sobre els materials. Amb l'encapsulació de la rapa de raïm es va incrementar significativament la capacitat d'adsorció del metall respecta a la seva utilització en la forma original. Les dades experimentals es van tractar amb un model de difusió i es va determinar el coeficient de difusió en l'adsorció en el sorbent amb diferent percentatge de rapa en l'interior. Finalment s'han realitzat estudis preliminars d'adsorció de crom hexavalent en continu utilitzant columnes de petites dimensions farcides amb el sorbent format per boles de rapa encapsulada en alginat càlcic. En aquests assajos s'ha determinat l'efecte en l'adsorció de metall de diversos paràmetres d'operació, com la concentració inicial de metall, el cabal i l'alçada de llit. El model de transferència de matèria i difusió en columna emprat pel tractament de les dades experimentals va mostrar l'efecte de la concentració inicial i del cabal en els coeficients de difusió i transferència de matèria. El model BDST no va predir correctament el temps de servei en els diferents paràmetres d'operació experimentats a causa de la complexitat del procés d'adsorció. Amb aquests assajos preliminars s'han establert les bases de l'adsorció de Cr(VI) amb el nou sorbent en sistemes en continu amb la finalitat de realitzar estudis posteriors d'eliminació del metall en columnes de majors dimensions.

Millora de la qualitat espermàtica de dosis seminals de mascles reproductors Piétrain per filtració en diverses reïnes

En aquest treball s'ha dissenyat un mètode ràpid i fiable de tinció amb fluorocroms per a l'anàlisi de la integritat i viabilitat espermàtiques a partir del marcatge de la beina mitocondrial amb MitoTracker®Green FM, de l'acrosoma amb la lectina Trypsin inhibitor from Soybean (SBTI) conjugada amb el fluorocrom Alexa Fluor®488 específic per la proacrosina i del nucli amb els fluorocroms bis-benzimida (específic per a cèl·lules viables) i iodur de propidi (específic per a cèl·lules no viables). També s'ha determinat l'efecte de la filtració de dosis seminals de mascles astentoteratonecrospèrmics en columnes de Sephadex neutre i de dosis de mascles amb baixa qualitat espermàtica per filtració en columnes de Sephadex iònic, llana de vidre i glass beads sobre la qualitat espermàtica dels diferents grups de mascles analitzats. Els resultats obtinguts han mostrat que diversos paràmetres de qualitat espermàtica milloren després de la filtració en les diferents reïnes segons la patologia que presentin.

Desenvolupament de mètodes de preconcentració emprant membranes líquides suportades i extracció en fase sòlida per a la determinació de l'herbicida glifosat i el seu metabòlit AMPA en aigües naturals

El glifosat, N-(fosfonometil) glicina, és un dels herbicides més utilitzats arreu del món a causa de la seva baixa toxicitat i al seu ampli espectre d'aplicació. A conseqüència del gran ús que se'n fa, és necessari monitoritzar aquest compost i el seu principal metabòlit, l'àcid aminometilfosfònic (AMPA), en el medi ambient. S'han descrit diversos mètodes instrumentals basats en cromatografia de gasos (GC) i de líquids (HPLC), sent aquesta darrera l'opció més favorable a causa del caràcter polar dels anàlits. Per assolir nivells de concentració baixos cal, però, la preconcentració dels anàlits. En aquest treball s'estudien diferents alternatives amb aquest objectiu. S'ha avaluat la tècnica de membrana líquida suportada (SLM) on la membrana consisteix en una dissolució orgànica, que conté un transportador (en el nostre cas, un bescanviador d'anions comercial, Aliquat 336), que impregna un suport polimèric microporós que se situa entre dues solucions aquoses: la de càrrega, que conté els anàlits inicialment, i la receptora, on es retenen els anàlits després del seu transport a través de la membrana. Les condicions d'extracció més adequades s'obtenen treballant en medi bàsic amb NaOH on els anàlits estan en forma aniònica i les majors recuperacions s'obtenen amb HCl 0,1 M o NaCl 0,5 M, la qual cosa indica que l'ió clorur és la força impulsora del transport. Un cop dissenyat el sistema, es duen a terme experiments de preconcentració amb dues geometries diferents: un sistema de membrana laminar (LSLM) on recircula la fase receptora i un sistema de fibra buida (HFSLM). Els millors resultats s'obtenen amb el mòdul de fibra buida, amb factors de concentració de 25 i 3 per a glifosat i AMPA, respectivament, fent recircular durant 24 hores 100 ml de solució de càrrega i 4 ml de solució receptora. També s'aplica una tècnica més selectiva, la cromatografia d'afinitat amb ió metàl·lic immobilitzat (IMAC), basada en la interacció entre els anàlits i un metall immobilitzat en una resina a través d'un grup funcional d'aquesta. En aquest estudi s'immobilitza pal·ladi al grup funcional 8-hidroxiquinoleïna de la resina amb matriu acrílica Spheron Oxine 1000 i s'avalua per a l'extracció i preconcentració de glifosat i AMPA. Per a ambdós anàlits l'adsorció és del 100 % i les recuperacions són superiors al 80 % i al 60 % per a glifosat i AMPA, respectivament, utilitzant HCl 0,1 M + NaCl 1 M com a eluent. Aquests resultats es comparen amb els obtinguts amb dues resines més, també carregades amb pal·ladi: Iontosorb Oxin 100, que té el mateix grup funcional però matriu de cel·lulosa, i Spheron Thiol 1000, on el grup funcional és un tiol i la matriu també és acrílica. Per al glifosat els resultats són similars amb totes les resines, però per a l'AMPA la resina Spheron Thiol és la única que proporciona recuperacions superiors al 93 %. Finalment, una altra opció estudiada és l'acoblament de dues columnes de cromatografia líquida (LC-LC). En l'estudi l'objectiu és millorar el mètode existent per a glifosat i AMPA en aigües naturals on el LOD era de 0,25 ug/l. El mètode consisteix en la derivatització precolumna amb el reactiu fluorescent FMOC i l'anàlisi amb l'acoblament LC-LC-fluorescència. Variant lleugerament les condicions de derivatització s'aconsegueix quantificar 0,1 ug/l de glifosat i AMPA. Es fortifiquen aigües naturals amb 0,1, 1 i 10 ug/l dels anàlits per validar el mètode. S'obtenen recuperacions d'entre el 85 % i el 100 %, amb desviacions estàndard relatives inferiors al 8 %. Aplicant una tècnica de preconcentració prèvia a la derivatització i anàlisi utilitzant una resina de bescanvi aniònic, Amberlite IRA-900, es millora la sensibilitat del mètode i s'assoleix un LOD per al glifosat de 0,02 ug/l.