Distributional equivalence and subcompositional coherence in the analysis of contingency tables, ratio-scale measurements and compositional data

We consider two fundamental properties in the analysis of two-way tables of positive data: the principle of distributional equivalence, one of the cornerstones of correspondence analysis of contingency tables, and the principle of subcompositional coherence, which forms the basis of compositional data analysis. For an analysis to be subcompositionally coherent, it suffices to analyse the ratios of the data values. The usual approach to dimension reduction in compositional data analysis is to perform principal component analysis on the logarithms of ratios, but this method does not obey the principle of distributional equivalence. We show that by introducing weights for the rows and columns, the method achieves this desirable property. This weighted log-ratio analysis is theoretically equivalent to spectral mapping , a multivariate method developed almost 30 years ago for displaying ratio-scale data from biological activity spectra. The close relationship between spectral mapping and correspondence analysis is also explained, as well as their connection with association modelling. The weighted log-ratio methodology is applied here to frequency data in linguistics and to chemical compositional data in archaeology.

Weighted metric multidimensional scaling

This paper establishes a general framework for metric scaling of any distance measure between individuals based on a rectangular individuals-by-variables data matrix. The method allows visualization of both individuals and variables as well as preserving all the good properties of principal axis methods such as principal components and correspondence analysis, based on the singular-value decomposition, including the decomposition of variance into components along principal axes which provide the numerical diagnostics known as contributions. The idea is inspired from the chi-square distance in correspondence analysis which weights each coordinate by an amount calculated from the margins of the data table. In weighted metric multidimensional scaling (WMDS) we allow these weights to be unknown parameters which are estimated from the data to maximize the fit to the original distances. Once this extra weight-estimation step is accomplished, the procedure follows the classical path in decomposing a matrix and displaying its rows and columns in biplots.

Analysis of matched matrices

We consider the joint visualization of two matrices which have common rowsand columns, for example multivariate data observed at two time pointsor split accord-ing to a dichotomous variable. Methods of interest includeprincipal components analysis for interval-scaled data, or correspondenceanalysis for frequency data or ratio-scaled variables on commensuratescales. A simple result in matrix algebra shows that by setting up thematrices in a particular block format, matrix sum and difference componentscan be visualized. The case when we have more than two matrices is alsodiscussed and the methodology is applied to data from the InternationalSocial Survey Program.

An enhanced concave program relaxation for choice network revenue management

The network choice revenue management problem models customers as choosing from an offer-set, andthe firm decides the best subset to offer at any given moment to maximize expected revenue. The resultingdynamic program for the firm is intractable and approximated by a deterministic linear programcalled the CDLP which has an exponential number of columns. However, under the choice-set paradigmwhen the segment consideration sets overlap, the CDLP is difficult to solve. Column generation has beenproposed but finding an entering column has been shown to be NP-hard. In this paper, starting with aconcave program formulation based on segment-level consideration sets called SDCP, we add a class ofconstraints called product constraints, that project onto subsets of intersections. In addition we proposea natural direct tightening of the SDCP called ?SDCP, and compare the performance of both methodson the benchmark data sets in the literature. Both the product constraints and the ?SDCP method arevery simple and easy to implement and are applicable to the case of overlapping segment considerationsets. In our computational testing on the benchmark data sets in the literature, SDCP with productconstraints achieves the CDLP value at a fraction of the CPU time taken by column generation and webelieve is a very promising approach for quickly approximating CDLP when segment consideration setsoverlap and the consideration sets themselves are relatively small.

Two-stage, extreme albitization of A-type granites from Rajasthan, NW India

Albitization is a common process during which hydrothermal fluids convert plagioclase and/or K-feldspar into nearly pure albite; however, its specific mechanism in granitoids is not well understood. The c. 1700 Ma A-type metaluminous ferroan granites in the Khetri complex of Rajasthan, NW India, have been albitized to a large extent by two metasomatic fronts, an initial transformation of oligoclase to nearly pure albite and a subsequent replacement of microcline by albite, with sharp contacts between the microcline-bearing and microcline-free zones. Albitization has bleached the original pinkish grey granite and turned it white. The mineralogical changes include transformation of oligoclase (similar to An(12)) and microcline (similar to Or(95)) to almost pure albite (similar to An(0 center dot 5-2)), amphibole from potassian ferropargasite (X-Fe 0 center dot 84-0 center dot 86) to potassic hastingsite (X-Fe 0 center dot 88-0 center dot 97) and actinolite (X-Fe 0 center dot 32-0 center dot 67), and biotite from annite (X-Fe 0 center dot 71-0 center dot 74) to annite (X-Fe 0 center dot 90-0 center dot 91). Whole-rock isocon diagrams show that, during albitization, the granites experienced major hydration, slight gain in Si and major gain in Na, whereas K, Mg, Fe and Ca were lost along with Rb, Ba, Sr, Zn, light rare earth elements and U. Whole-rock Sm-Nd isotope data plot on an apparent isochron of 1419 +/- 98 Ma and reveal significant disturbance and at least partial resetting of the intrusion age. Severe scatter in the whole-rock Rb-Sr isochron plot reflects the extreme Rb loss in the completely albitized samples, effectively freezing Sr-87/Sr-86 ratios in the albite granites at very high values (0 center dot 725-0 center dot 735). This indicates either infiltration of highly radiogenic Sr from the country rock or, more likely, radiogenic ingrowth during a considerable time lag (estimated to be at least 300 Myr) between original intrusion and albitization. The albitization took place at similar to 350-400 degrees C. It was caused by the infiltration of an ascending hydrothermal fluid that had acquired high Na/K and Na/Ca ratios during migration through metamorphic rocks at even lower temperatures in the periphery of the plutons. Oxygen isotope ratios increase from delta O-18 = 7 parts per thousand in the original granite to values of 9-10 parts per thousand in completely albitized samples, suggesting that the fluid had equilibrated with surrounding metamorphosed crust. A metasomatic model, using chromatographic theory of fluid infiltration, explains the process for generating the observed zonation in terms of a leading metasomatic front where oligoclase of the original granite is converted to albite, and a second, trailing front where microcline is also converted to albite. The temperature gradients driving the fluid infiltration may have been produced by the high heat production of the granites themselves. The confinement of the albitized granites along the NE-SW-trending Khetri lineament and the pervasive nature of the albitization suggest that the albitizing fluids possibly originated during reactivation of the lineament. More generally, steady-state temperature gradients induced by the high internal heat production of A-type granites may provide the driving force for similar metasomatic and ore-forming processes in other highly enriched granitoid bodies.

Synthesis of a radioiodinated photoreactive MAGE-1 peptide derivative and photoaffinity labeling of cell-associated human leukocyte antigen-A1 molecules.

The synthesis of a photoreactive derivative of the human leukocyte antigen-A1 (HLA-A1)-restricted MAGE-1 peptide 161-169 (EADPTGHSY) is described. Using conventional automated solid-phase peptide synthesis, a photoreactive derivative of this peptide was synthesized by replacing histidine-167 with photo-reactive N-beta-4-azidosalicyloyl-L-2,3-diaminopropionic acid. The C-terminal tyrosine was incorporated as phosphotyrosine. This peptide derivative was radioiodinated in the presence of chloramine T. This iodination took place selectively at the photoreactive group, because the phosphate ester prevented tyrosine iodination. Following dephosphorylation with alkaline phosphatase and chromatographic purification, the radiolabeled peptide derivative was incubated with cells expressing HLA-A1 or other HLA molecules. Photoactivation resulted in efficient photoaffinity labeling of HLA-A1. Other HLA molecules or other cellular components were not detectably labeled. This labeling was inhibited by HLA-A1 but not by HLA-A2-binding peptides. This synthesis is generally applicable and can also be adapted to the synthesis of well-defined radiolabeled nonphotoreactive peptide derivatives.

Correspondence analysis of two transition tables

The case of two transition tables is considered, that is two squareasymmetric matrices of frequencies where the rows and columns of thematrices are the same objects observed at three different timepoints. Different ways of visualizing the tables, either separatelyor jointly, are examined. We generalize an existing idea where asquare matrix is descomposed into symmetric and skew-symmetric partsto two matrices, leading to a decomposition into four components: (1)average symmetric, (2) average skew-symmetric, (3) symmetricdifference from average, and (4) skew-symmetric difference fromaverage. The method is illustrated with an artificial example and anexample using real data from a study of changing values over threegenerations.

Weighted Euclidean biplots

We construct a weighted Euclidean distance that approximates any distance or dissimilarity measure between individuals that is based on a rectangular cases-by-variables data matrix. In contrast to regular multidimensional scaling methods for dissimilarity data, the method leads to biplots of individuals and variables while preserving all the good properties of dimension-reduction methods that are based on the singular-value decomposition. The main benefits are the decomposition of variance into components along principal axes, which provide the numerical diagnostics known as contributions, and the estimation of nonnegative weights for each variable. The idea is inspired by the distance functions used in correspondence analysis and in principal component analysis of standardized data, where the normalizations inherent in the distances can be considered as differential weighting of the variables. In weighted Euclidean biplots we allow these weights to be unknown parameters, which are estimated from the data to maximize the fit to the chosen distances or dissimilarities. These weights are estimated using a majorization algorithm. Once this extra weight-estimation step is accomplished, the procedure follows the classical path in decomposing the matrix and displaying its rows and columns in biplots.

Multicentre validation study of nucleic acids extraction from FFPE tissues.

In most pathology laboratories worldwide, formalin-fixed paraffin embedded (FFPE) samples are the only tissue specimens available for routine diagnostics. Although commercial kits for diagnostic molecular pathology testing are becoming available, most of the current diagnostic tests are laboratory-based assays. Thus, there is a need for standardized procedures in molecular pathology, starting from the extraction of nucleic acids. To evaluate the current methods for extracting nucleic acids from FFPE tissues, 13 European laboratories, participating to the European FP6 program IMPACTS (www.impactsnetwork.eu), isolated nucleic acids from four diagnostic FFPE tissues using their routine methods, followed by quality assessment. The DNA-extraction protocols ranged from homemade protocols to commercial kits. Except for one homemade protocol, the majority gave comparable results in terms of the quality of the extracted DNA measured by the ability to amplify differently sized control gene fragments by PCR. For array-applications or tests that require an accurately determined DNA-input, we recommend using silica based adsorption columns for DNA recovery. For RNA extractions, the best results were obtained using chromatography column based commercial kits, which resulted in the highest quantity and best assayable RNA. Quality testing using RT-PCR gave successful amplification of 200 bp-250 bp PCR products from most tested tissues. Modifications of the proteinase-K digestion time led to better results, even when commercial kits were applied. The results of the study emphasize the need for quality control of the nucleic acid extracts with standardised methods to prevent false negative results and to allow data comparison among different diagnostic laboratories.

2002 Iowa Youth Survey State of Iowa Results, April 2003

The 2002 Iowa Youth Survey (IYS) State of Iowa report was designed to help state-level planners identify youth development-related needs, develop relevant programs, and assess the outcomes of those programs. These data can help us better understand our youth and their needs. They can help us assess the strengths and weaknesses of our schools, families and communities from the young person’s perspective. In addition, the data in this report help the state obtain funds for a wide variety of programs. At every step in the process – from needs identification, to program development and implementation, to program assessment – the 2002 IYS data will provide a valuable resource. The state report can also help Iowa’s schools, area education agencies and counties assess their relative strengths and weaknesses. The grades 6, 8, and 11, as well as male and female percentages reported in each of these reports can be compared with the respective state report percentages. The higher the proportion of students in each of these columns that completed a usable IYS questionnaire, the more likely the comparisons with the state report percentages will be unbiased. Such comparisons should be considered exploratory, but for the most part are likely to prove useful.

2005 Iowa Youth Survey State of Iowa Results, Revised December 2006

The 2005 Iowa Youth Survey (IYS) State of Iowa report was designed to help state-level planners identify youth development-related needs, develop relevant programs, and assess the outcomes of those programs. These data can help you better understand our youth and their needs. They can also help you assess the strengths and weaknesses of our schools, families, and communities from the young person’s perspective. In addition, the data in this report help the state obtain funds for a wide variety of programs. At every step of the process–from needs identification, to program development and implementation, to program assessment–the 2005 IYS data will prove to be a valuable resource. The state report can also help Iowa’s schools, area education agencies, and counties assess their relative strengths and weaknesses. The grades 6, 8, and 11, as well as male and female percentages reported in district-level, AEA-level, county-level, and other 2005 IYS reports can be compared with the respective state report percentages. The higher the proportion of students in each of these columns that completed usable IYS questionnaires, the more likely the comparisons with the state report percentages will be unbiased. Such comparisons should be considered exploratory, but for the most part are likely to be useful.

Isolation and characterization of carcinoembryonic antigen (CEA) extracted from normal human colon mucosa.

Carcinoembryonic antigen (CEA) was identified in perchloric acid (PCA)_extract from normal colon mucosa by 2 immunological criteria: a line of identity in double diffusion and a parallel inhibition curve in radioimmunoassay (RIA), both with reference colon carcinoma-CEA (CEA-Tu). The average concentration of CEA in normal colon mucosa (CEA-No) was 35 times lower than in primary large bowel carcinomas and 230 times lower than in metastatic colon or rectum carcinomas. CEA-No was purified from PCA extracts of normal colon mucosa by Sephadex G-200 filtration and immunoadsorbent columns. Purified CEA-No had quatitatively the same inhibition activity in RIA as the British Standard CEA coded 73/601. Purified CEA-No was labelled with 125I. The percentage of binding of labelled CEA-No to a specific goat anti-CEA-Tu antiserum was similar to that of CEA-Tu. Labelled CEA-No could be used as radioactive tracer in RIA as well as labelled CEA-Tu. The physico-chemical properties of purified CEA-Tu as demonstrated by Sepharose 6 B filtration, SDS Polyacrylamide gel analysis and cesium chloride density gradient, were found to be almost identical to those of reference CEA-Tu. Preliminary results showed that CEA-No and CEA-Tu contained the same types of carbohydrates in similar proportions. A rabbit antiserum against CEA-No was obtained which demonstrated the same specificity as conventional anti-CEA-Tu antisera.

Mobility of inorganic and organic phosphorus forms under different levels of phosphate and poultry litter fertilization in soils

The eutrophication of aquifers is strongly linked to the mobility of P in soils. Although P mobility was considered irrelevant in a more distant past, more recent studies have shown that P, both in organic (Po) and inorganic forms (Pi), can be lost by leaching and eluviation through the soil profile, particularly in less weathered and/or sandier soils with low P adsorption capacity. The purpose of this study was to determine losses of P forms by leaching and eluviation from soil columns. Each column consisted of five PVC rings (diameter 5 cm, height 10 cm), filled with two soil types: a clayey Red-Yellow Latosol and a sandy loam Red-Yellow Latosol, which were exposed to water percolation. The soils were previously treated with four P rates (as KH2PO4 ) to reach 0, 12.5, 25.0 and 50 % of the maximum P adsorption capacity (MPAC). The P source was homogenized with the whole soil volume and incubated for 60 days. After this period the soils were placed in the columns; the soil of the top ring was mixed with five poultry litter rates of 0, 20, 40, 80, and 160 t ha-1 (dry weight basis). Treatments consisted of a 4 x 5 x 2 factorial scheme corresponding to four MPAC levels, five poultry litter rates, two soils, with three replications, arranged in a completely randomized block design. Deionized water was percolated through the columns 10 times in 35 days to simulate about 1,200 mm rainfall. In the leachate of each column the inorganic P (reactive P, Pi) and organic P forms (unreactive P, Po) were determined. At the end of the experiment, the columns were disassembled and P was extracted with the extractants Mehlich-1 (HCl 0.05 mol L-1 and H2SO4 0.0125 mol L-1) and Olsen (NaHCO3 0.5 mol L-1; pH 8.5) from the soil of each ring. The Pi and Po fractions were measured by the Olsen extractant. It was found that under higher poultry litter rates the losses of unreactive P (Po) were 6.4 times higher than of reactive P (Pi). Both the previous P fertilization and increasing poultry litter rates caused a vertical movement of P down the soil columns, as verified by P concentrations extracted by Mehlich-1 and NaHCO3 (Olsen). The environmental critical level (ECL), i.e., the P soil concentration above which P leaching increases exponentially, was 100 and 150 mg dm-3 by Mehlich-1 and 40 and 60 mg dm-3 by Olsen, for the sandy loam and clay soils, respectively. In highly weathered soils, where residual P is accumulated by successive crops, P leaching through the profile can be significant, particularly when poultry litter is applied as fertilizer.

Spatial relationships between soil attributes and corn yield in no-tillage system

Soil properties play an important role in spatial variability of crop yield. However, a low spatial correlation has generally been observed between maps of crop yield and of soil properties. The objectives of the present investigation were to assess the spatial pattern variability of soil properties and of corn yield at the same sampling intensity, and evaluate its cause-and-effect relationships. The experimental site was structured in a grid of 100 referenced points, spaced at 10 m intervals along four parallel 250 m long rows spaced 4.5 m apart. Thus, points formed a rectangle containing four columns and 25 rows. Therefore, each sampling cell encompassed an area of 45 m² and consisted of five 10 m long crop rows, in which the referenced points represented the center. Samples were taken from the layers 0-0.1 m and 0.1-0.2 m. Soil physical and chemical properties were evaluated. Statistical analyses consisted of data description and geostatistics. The spatial dependence of corn yield and soil properties was confirmed. The hypothesis of this study was confirmed, i.e., when sampling the soil to determine the values of soil characteristics at similar to sampling intensity as for crop yield assessments, correlations between the spatial distribution of soil characteristics and crop yield were observed. The spatial distribution pattern of soil properties explained 65 % of the spatial distribution pattern of corn yield. The spatial distribution pattern of clay content and percentage of soil base saturation explained most of the spatial distribution pattern of corn yield.

ANIBAL, stable isotope-based quantitative proteomics by aniline and benzoic acid labeling of amino and carboxylic groups

Identification and relative quantification of hundreds to thousands of proteins within complex biological samples have become realistic with the emergence of stable isotope labeling in combination with high throughput mass spectrometry. However, all current chemical approaches target a single amino acid functionality (most often lysine or cysteine) despite the fact that addressing two or more amino acid side chains would drastically increase quantifiable information as shown by in silico analysis in this study. Although the combination of existing approaches, e.g. ICAT with isotope-coded protein labeling, is analytically feasible, it implies high costs, and the combined application of two different chemistries (kits) may not be straightforward. Therefore, we describe here the development and validation of a new stable isotope-based quantitative proteomics approach, termed aniline benzoic acid labeling (ANIBAL), using a twin chemistry approach targeting two frequent amino acid functionalities, the carboxylic and amino groups. Two simple and inexpensive reagents, aniline and benzoic acid, in their (12)C and (13)C form with convenient mass peak spacing (6 Da) and without chromatographic discrimination or modification in fragmentation behavior, are used to modify carboxylic and amino groups at the protein level, resulting in an identical peptide bond-linked benzoyl modification for both reactions. The ANIBAL chemistry is simple and straightforward and is the first method that uses a (13)C-reagent for a general stable isotope labeling approach of carboxylic groups. In silico as well as in vitro analyses clearly revealed the increase in available quantifiable information using such a twin approach. ANIBAL was validated by means of model peptides and proteins with regard to the quality of the chemistry as well as the ionization behavior of the derivatized peptides. A milk fraction was used for dynamic range assessment of protein quantification, and a bacterial lysate was used for the evaluation of relative protein quantification in a complex sample in two different biological states