Clustering compositional data trajectories

Our essay aims at studying suitable statistical methods for the clustering ofcompositional data in situations where observations are constituted by trajectories ofcompositional data, that is, by sequences of composition measurements along a domain.Observed trajectories are known as “functional data” and several methods have beenproposed for their analysis.In particular, methods for clustering functional data, known as Functional ClusterAnalysis (FCA), have been applied by practitioners and scientists in many fields. To ourknowledge, FCA techniques have not been extended to cope with the problem ofclustering compositional data trajectories. In order to extend FCA techniques to theanalysis of compositional data, FCA clustering techniques have to be adapted by using asuitable compositional algebra.The present work centres on the following question: given a sample of compositionaldata trajectories, how can we formulate a segmentation procedure giving homogeneousclasses? To address this problem we follow the steps described below.First of all we adapt the well-known spline smoothing techniques in order to cope withthe smoothing of compositional data trajectories. In fact, an observed curve can bethought of as the sum of a smooth part plus some noise due to measurement errors.Spline smoothing techniques are used to isolate the smooth part of the trajectory:clustering algorithms are then applied to these smooth curves.The second step consists in building suitable metrics for measuring the dissimilaritybetween trajectories: we propose a metric that accounts for difference in both shape andlevel, and a metric accounting for differences in shape only.A simulation study is performed in order to evaluate the proposed methodologies, usingboth hierarchical and partitional clustering algorithm. The quality of the obtained resultsis assessed by means of several indices

Possible solution of some essential zero problems in compositional data analysis

One of the tantalising remaining problems in compositional data analysis lies in how to deal with data sets in which there are components which are essential zeros. By anessential zero we mean a component which is truly zero, not something recorded as zero simply because the experimental design or the measuring instrument has not been sufficiently sensitive to detect a trace of the part. Such essential zeros occur inmany compositional situations, such as household budget patterns, time budgets,palaeontological zonation studies, ecological abundance studies. Devices such as nonzero replacement and amalgamation are almost invariably ad hoc and unsuccessful insuch situations. From consideration of such examples it seems sensible to build up amodel in two stages, the first determining where the zeros will occur and the secondhow the unit available is distributed among the non-zero parts. In this paper we suggest two such models, an independent binomial conditional logistic normal model and a hierarchical dependent binomial conditional logistic normal model. The compositional data in such modelling consist of an incidence matrix and a conditional compositional matrix. Interesting statistical problems arise, such as the question of estimability of parameters, the nature of the computational process for the estimation of both the incidence and compositional parameters caused by the complexity of the subcompositional structure, the formation of meaningful hypotheses, and the devising of suitable testing methodology within a lattice of such essential zero-compositional hypotheses. The methodology is illustrated by application to both simulated and real compositional data

Monitoring procedures in environmental geochemistry and compositional data analysis theory

First discussion on compositional data analysis is attributable to Karl Pearson, in 1897. However, notwithstanding the recent developments on algebraic structure of the simplex, more than twenty years after Aitchison’s idea of log-transformations of closed data, scientific literature is again full of statistical treatments of this type of data by using traditional methodologies. This is particularly true in environmental geochemistry where besides the problem of the closure, the spatial structure (dependence) of the data have to be considered. In this work we propose the use of log-contrast values, obtained by asimplicial principal component analysis, as LQGLFDWRUV of given environmental conditions. The investigation of the log-constrast frequency distributions allows pointing out the statistical laws able togenerate the values and to govern their variability. The changes, if compared, for example, with the mean values of the random variables assumed as models, or other reference parameters, allow definingmonitors to be used to assess the extent of possible environmental contamination. Case study on running and ground waters from Chiavenna Valley (Northern Italy) by using Na+, K+, Ca2+, Mg2+, HCO3-, SO4 2- and Cl- concentrations will be illustrated

The neuronal basis of attention: rate versus synchronization modulation

Extensive theoretical and experimental work on the neuronal correlates of visual attention raises two hypotheses about the underlying mechanisms. The first hypothesis, named biased competition, originates from experimental single-cell recordings that have shown that attention upmodulates the firing rates of the neurons encoding the attended features and downregulates the firing rates of the neurons encoding the unattended features. Furthermore, attentional modulation of firing rates increases along the visual pathway. The other, newer hypothesis assigns synchronization a crucial role in the attentional process. It stems from experiments that have shown that attention modulates gamma-frequency synchronization. In this paper, we study the coexistence of the two phenomena using a theoretical framework. We find that the two effects can vary independently of each other and across layers. Therefore, the two phenomena are not concomitant. However, we show that there is an advantage in the processing of information if rate modulation is accompanied by gamma modulation, namely that reaction times are shorter, implying behavioral relevance for gamma synchronization.

Silver space-time trellis-coded modulation

Silver Code (SilC) was originally discovered in [1–4] for 2×2 multiple-input multiple-output (MIMO) transmission. It has non-vanishing minimum determinant 1/7, slightly lower than Golden code, but is fast-decodable, i.e., it allows reduced-complexity maximum likelihood decoding [5–7]. In this paper, we present a multidimensional trellis-coded modulation scheme for MIMO systems [11] based on set partitioning of the Silver Code, named Silver Space-Time Trellis Coded Modulation (SST-TCM). This lattice set partitioning is designed specifically to increase the minimum determinant. The branches of the outer trellis code are labeled with these partitions. Viterbi algorithm is applied for trellis decoding, while the branch metrics are computed by using a sphere-decoding algorithm. It is shown that the proposed SST-TCM performs very closely to the Golden Space-Time Trellis Coded Modulation (GST-TCM) scheme, yetwith a much reduced decoding complexity thanks to its fast-decoding property.

Some last thoughts on compositional data analysis

One of the disadvantages of old age is that there is more past than future: this,however, may be turned into an advantage if the wealth of experience and, hopefully,wisdom gained in the past can be reflected upon and throw some light on possiblefuture trends. To an extent, then, this talk is necessarily personal, certainly nostalgic,but also self critical and inquisitive about our understanding of the discipline ofstatistics. A number of almost philosophical themes will run through the talk: searchfor appropriate modelling in relation to the real problem envisaged, emphasis onsensible balances between simplicity and complexity, the relative roles of theory andpractice, the nature of communication of inferential ideas to the statistical layman, theinter-related roles of teaching, consultation and research. A list of keywords might be:identification of sample space and its mathematical structure, choices betweentransform and stay, the role of parametric modelling, the role of a sample spacemetric, the underused hypothesis lattice, the nature of compositional change,particularly in relation to the modelling of processes. While the main theme will berelevance to compositional data analysis we shall point to substantial implications forgeneral multivariate analysis arising from experience of the development ofcompositional data analysis…

Compositional Data in Biomedical Research

Modern methods of compositional data analysis are not well known in biomedical research.Moreover, there appear to be few mathematical and statistical researchersworking on compositional biomedical problems. Like the earth and environmental sciences,biomedicine has many problems in which the relevant scienti c information isencoded in the relative abundance of key species or categories. I introduce three problemsin cancer research in which analysis of compositions plays an important role. Theproblems involve 1) the classi cation of serum proteomic pro les for early detection oflung cancer, 2) inference of the relative amounts of di erent tissue types in a diagnostictumor biopsy, and 3) the subcellular localization of the BRCA1 protein, and it'srole in breast cancer patient prognosis. For each of these problems I outline a partialsolution. However, none of these problems is \solved". I attempt to identify areas inwhich additional statistical development is needed with the hope of encouraging morecompositional data analysts to become involved in biomedical research

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.

Dynamic graphics of parametrically linked multivariate methods used in compositional data analysis

Many multivariate methods that are apparently distinct can be linked by introducing oneor more parameters in their definition. Methods that can be linked in this way arecorrespondence analysis, unweighted or weighted logratio analysis (the latter alsoknown as "spectral mapping"), nonsymmetric correspondence analysis, principalcomponent analysis (with and without logarithmic transformation of the data) andmultidimensional scaling. In this presentation I will show how several of thesemethods, which are frequently used in compositional data analysis, may be linkedthrough parametrizations such as power transformations, linear transformations andconvex linear combinations. Since the methods of interest here all lead to visual mapsof data, a "movie" can be made where where the linking parameter is allowed to vary insmall steps: the results are recalculated "frame by frame" and one can see the smoothchange from one method to another. Several of these "movies" will be shown, giving adeeper insight into the similarities and differences between these methods.

Biplots of compositional data

The singular value decomposition and its interpretation as alinear biplot has proved to be a powerful tool for analysing many formsof multivariate data. Here we adapt biplot methodology to the specifficcase of compositional data consisting of positive vectors each of whichis constrained to have unit sum. These relative variation biplots haveproperties relating to special features of compositional data: the studyof ratios, subcompositions and models of compositional relationships. Themethodology is demonstrated on a data set consisting of six-part colourcompositions in 22 abstract paintings, showing how the singular valuedecomposition can achieve an accurate biplot of the colour ratios and howpossible models interrelating the colours can be diagnosed.

Neuroprotective role of PrPC against kainate-induced epileptic seizures and cell death depends on the modulation of JNK3 activation by GluR6/7-PSD-95 binding

Cellular prion protein (PrPC) is a glycosyl-phosphatidylinositol¿anchored glycoprotein. When mutated or misfolded, the pathogenic form (PrPSC) induces transmissible spongiform encephalopathies. In contrast, PrPC has a number of physiological functions in several neural processes. Several lines of evidence implicate PrPC in synaptic transmission and neuroprotection since its absence results in an increase in neuronal excitability and enhanced excitotoxicity in vitro and in vivo. Furthermore, PrPC has been implicated in the inhibition of N-methyl-D-aspartic acid (NMDA)¿mediated neurotransmission, and prion protein gene (Prnp) knockout mice show enhanced neuronal death in response to NMDA and kainate (KA). In this study, we demonstrate that neurotoxicity induced by KA in Prnp knockout mice depends on the c-Jun N-terminal kinase 3 (JNK3) pathway since Prnpo/oJnk3o/o mice were not affected by KA. Pharmacological blockage of JNK3 activity impaired PrPC-dependent neurotoxicity. Furthermore, our results indicate that JNK3 activation depends on the interaction of PrPC with postsynaptic density 95 protein (PSD-95) and glutamate receptor 6/7 (GluR6/7). Indeed, GluR6¿PSD-95 interaction after KA injections was favored by the absence of PrPC. Finally, neurotoxicity in Prnp knockout mice was reversed by an AMPA/KA inhibitor (6,7-dinitroquinoxaline-2,3-dione) and the GluR6 antagonist NS-102. We conclude that the protection afforded by PrPC against KA is due to its ability to modulate GluR6/7-mediated neurotransmission and hence JNK3 activation.

Dynamic and Regulated Association of Caveolin with Lipid Bodies: Modulation of Lipid Body Motility and Function by a Dominant Negative Mutant

Caveolins are a crucial component of caveolae but have also been localized to the Golgi complex, and, under some experimental conditions, to lipid bodies (LBs). The physiological relevance and dynamics of LB association remain unclear. We now show that endogenous caveolin-1 and caveolin-2 redistribute to LBs in lipid loaded A431 and FRT cells. Association with LBs is regulated and reversible; removal of fatty acids causes caveolin to rapidly leave the lipid body. We also show by subcellular fractionation, light and electron microscopy that during the first hours of liver regeneration, caveolins show a dramatic redistribution from the cell surface to the newly formed LBs. At later stages of the regeneration process (when LBs are still abundant), the levels of caveolins in LBs decrease dramatically. As a model system to study association of caveolins with LBs we have used brefeldin A (BFA). BFA causes rapid redistribution of endogenous caveolins to LBs and this association was reversed upon BFA washout. Finally, we have used a dominant negative LB-associated caveolin mutant (cavDGV) to study LB formation and to examine its effect on LB function. We now show that the cavDGV mutant inhibits microtubule-dependent LB motility and blocks the reversal of lipid accumulation in LBs.

Neuroprotective role of PrPC against kainate-induced epileptic seizures and cell death depends on the modulation of JNK3 activation by GluR6/7-PSD-95 binding

Cellular prion protein (PrPC) is a glycosyl-phosphatidylinositol¿anchored glycoprotein. When mutated or misfolded, the pathogenic form (PrPSC) induces transmissible spongiform encephalopathies. In contrast, PrPC has a number of physiological functions in several neural processes. Several lines of evidence implicate PrPC in synaptic transmission and neuroprotection since its absence results in an increase in neuronal excitability and enhanced excitotoxicity in vitro and in vivo. Furthermore, PrPC has been implicated in the inhibition of N-methyl-D-aspartic acid (NMDA)¿mediated neurotransmission, and prion protein gene (Prnp) knockout mice show enhanced neuronal death in response to NMDA and kainate (KA). In this study, we demonstrate that neurotoxicity induced by KA in Prnp knockout mice depends on the c-Jun N-terminal kinase 3 (JNK3) pathway since Prnpo/oJnk3o/o mice were not affected by KA. Pharmacological blockage of JNK3 activity impaired PrPC-dependent neurotoxicity. Furthermore, our results indicate that JNK3 activation depends on the interaction of PrPC with postsynaptic density 95 protein (PSD-95) and glutamate receptor 6/7 (GluR6/7). Indeed, GluR6¿PSD-95 interaction after KA injections was favored by the absence of PrPC. Finally, neurotoxicity in Prnp knockout mice was reversed by an AMPA/KA inhibitor (6,7-dinitroquinoxaline-2,3-dione) and the GluR6 antagonist NS-102. We conclude that the protection afforded by PrPC against KA is due to its ability to modulate GluR6/7-mediated neurotransmission and hence JNK3 activation.

Direct modulation of electroluminescence from silicon nanocrystals beyond radiative recombination rates

We propose a light emitting transistor based on silicon nanocrystals provided with 200 Mbits/ s built-in modulation. Suppression of electroluminescence from silicon nanocrystals embedded into the gate oxide of a field effect transistor is achieved by fast Auger quenching. In this process, a modulating drain signal causes heating of carriers in the channel and facilitates the charge injection into the nanocrystals. This excess of charge enables fast nonradiative processes that are used to obtain 100% modulation depths at modulating voltages of 1 V.