Adaptively weighted maximum likelihood estimation of discrete distributions

SummaryDiscrete data arise in various research fields, typically when the observations are count data.I propose a robust and efficient parametric procedure for estimation of discrete distributions. The estimation is done in two phases. First, a very robust, but possibly inefficient, estimate of the model parameters is computed and used to indentify outliers. Then the outliers are either removed from the sample or given low weights, and a weighted maximum likelihood estimate (WML) is computed.The weights are determined via an adaptive process such that if the data follow the model, then asymptotically no observation is downweighted.I prove that the final estimator inherits the breakdown point of the initial one, and that its influence function at the model is the same as the influence function of the maximum likelihood estimator, which strongly suggests that it is asymptotically fully efficient.The initial estimator is a minimum disparity estimator (MDE). MDEs can be shown to have full asymptotic efficiency, and some MDEs have very high breakdown points and very low bias under contamination. Several initial estimators are considered, and the performances of the WMLs based on each of them are studied.It results that in a great variety of situations the WML substantially improves the initial estimator, both in terms of finite sample mean square error and in terms of bias under contamination. Besides, the performances of the WML are rather stable under a change of the MDE even if the MDEs have very different behaviors.Two examples of application of the WML to real data are considered. In both of them, the necessity for a robust estimator is clear: the maximum likelihood estimator is badly corrupted by the presence of a few outliers.This procedure is particularly natural in the discrete distribution setting, but could be extended to the continuous case, for which a possible procedure is sketched.RésuméLes données discrètes sont présentes dans différents domaines de recherche, en particulier lorsque les observations sont des comptages.Je propose une méthode paramétrique robuste et efficace pour l'estimation de distributions discrètes. L'estimation est faite en deux phases. Tout d'abord, un estimateur très robuste des paramètres du modèle est calculé, et utilisé pour la détection des données aberrantes (outliers). Cet estimateur n'est pas nécessairement efficace. Ensuite, soit les outliers sont retirés de l'échantillon, soit des faibles poids leur sont attribués, et un estimateur du maximum de vraisemblance pondéré (WML) est calculé.Les poids sont déterminés via un processus adaptif, tel qu'asymptotiquement, si les données suivent le modèle, aucune observation n'est dépondérée.Je prouve que le point de rupture de l'estimateur final est au moins aussi élevé que celui de l'estimateur initial, et que sa fonction d'influence au modèle est la même que celle du maximum de vraisemblance, ce qui suggère que cet estimateur est pleinement efficace asymptotiquement.L'estimateur initial est un estimateur de disparité minimale (MDE). Les MDE sont asymptotiquement pleinement efficaces, et certains d'entre eux ont un point de rupture très élevé et un très faible biais sous contamination. J'étudie les performances du WML basé sur différents MDEs.Le résultat est que dans une grande variété de situations le WML améliore largement les performances de l'estimateur initial, autant en terme du carré moyen de l'erreur que du biais sous contamination. De plus, les performances du WML restent assez stables lorsqu'on change l'estimateur initial, même si les différents MDEs ont des comportements très différents.Je considère deux exemples d'application du WML à des données réelles, où la nécessité d'un estimateur robuste est manifeste : l'estimateur du maximum de vraisemblance est fortement corrompu par la présence de quelques outliers.La méthode proposée est particulièrement naturelle dans le cadre des distributions discrètes, mais pourrait être étendue au cas continu.

Ab initio absorption spectrum of NiO combining molecular dynamics with the embedded cluster approach in a discrete reaction field

We developed a procedure that combines three complementary computational methodologies to improve the theoretical description of the electronic structure of nickel oxide. The starting point is a Car-Parrinello molecular dynamics simulation to incorporate vibrorotational degrees of freedom into the material model. By means ofcomplete active space self-consistent field second-order perturbation theory (CASPT2) calculations on embedded clusters extracted from the resulting trajectory, we describe localized spectroscopic phenomena on NiO with an efficient treatment of electron correlation. The inclusion of thermal motion into the theoretical description allowsus to study electronic transitions that, otherwise, would be dipole forbidden in the ideal structure and results in a natural reproduction of the band broadening. Moreover, we improved the embedded cluster model by incorporating self-consistently at the complete active space self-consistent field (CASSCF) level a discrete (or direct) reaction field (DRF) in the cluster surroundings. The DRF approach offers an efficient treatment ofelectric response effects of the crystalline embedding to the electronic transitions localized in the cluster. We offer accurate theoretical estimates of the absorption spectrum and the density of states around the Fermi level of NiO, and a comprehensive explanation of the source of the broadening and the relaxation of the charge transferstates due to the adaptation of the environment

Implementació d'una plataforma client/servidor i un sistema gestor d'events pel projecte M-CISNE

El projecte M-CISNE es presenta con una eina de planificació per entorns multiclúster no dedicats i heterogenis, amb una aquitectura jeràrquica de dos nivells. En aquest treball es vol dotar a aquest sistema d'una plataforma client/servidor que permeti la comunicació remota entre els dos subsitemes a través de la xarxa i d'una eina que gestioni els events generats per cadascun d'ells.

Rates and pattern of evolution among Cretaceous radiolarians: relations with global paleoceanographic events

We present a new approach for analyzing the turnover rates of Cretaceous radiolarians recorded in pelagic sequences of western Tethys, The analysis of major extinction-radiation events and the fluctuation of diversity are compared with major paleoceanographic events and variation of diversity in dinoflagellates, calcareous nannoplankton and ammonites. There is an extraordinary correlation between biotic changes and sea level changes, temperatures, O, C and Sr isotopes, phosphorus accumulation rates and anoxic episodes. This reveals a predominantly abiotic control on the evolution of radiolarians. The rate of turnover and the diversity through time of two major orders of radiolarians (nassellarians and spumellarians) exhibits (1) the quasi-parallelism of their diversity curves, excluding a direct competition between them, (2) greater resistance of spumellarians to extinction during the early stage of extinction intervals and (3) a stronger post-extinction recovery of nassellarians. Evolutionary rates of radiolarians can be a good means of monitoring global environmental changes and allowing us to understand more clearly the relationship between plankton evolution, climate and pale oceanographic processes.

Adverse life events among community-dwelling persons aged 65-70 years: gender differences in occurrence and perceived psychological consequences.

OBJECTIVES: To describe the occurrence of selected adverse life events in young-old men and women, as well as their perceived psychological consequences. METHODS: In 2005, 1,422 participants in the Lausanne Cohort 65+ study, born in 1934-1938, self-reported whether they experienced any of 26 life events during the preceding year. Most participants (N = 1,309, 92%) completed the geriatric adverse life events scale during a face-to-face interview, by rating the level of stress associated with each event, as well as its impact on their psychological well-being. RESULTS: Overall, 72% of the participants experienced at least one of the 26 events in the preceding year (range 1-9). Disease affecting the respondent (N = 525) or a close relative (N = 276) was most frequent, as well as the death of a friend or non-close relative (N = 274). Women indicated a higher frequency of events (mean 2.1 vs. 1.7 events, P < 0.001), as well as a higher level of stress and a stronger negative impact on well-being than men. In multivariate analyses adjusting for self-rated health, depressive symptoms and comorbidity, female gender remained significantly associated with the level of stress and negative impact on psychological well-being. CONCLUSION: This exploratory study shows that several types of adverse life events frequently occur at age 65-70, with gender differences both in the frequency of reporting and consequences of these events. However, information on this topic is limited and studies based on different populations and designs are needed to better understand the impact of such events.

Subclinical thyroid dysfunction and the risk of heart failure events: an individual participant data analysis from 6 prospective cohorts.

BACKGROUND: American College of Cardiology/American Heart Association guidelines for the diagnosis and management of heart failure recommend investigating exacerbating conditions such as thyroid dysfunction, but without specifying the impact of different thyroid-stimulation hormone (TSH) levels. Limited prospective data exist on the association between subclinical thyroid dysfunction and heart failure events. METHODS AND RESULTS: We performed a pooled analysis of individual participant data using all available prospective cohorts with thyroid function tests and subsequent follow-up of heart failure events. Individual data on 25 390 participants with 216 248 person-years of follow-up were supplied from 6 prospective cohorts in the United States and Europe. Euthyroidism was defined as TSH of 0.45 to 4.49 mIU/L, subclinical hypothyroidism as TSH of 4.5 to 19.9 mIU/L, and subclinical hyperthyroidism as TSH <0.45 mIU/L, the last two with normal free thyroxine levels. Among 25 390 participants, 2068 (8.1%) had subclinical hypothyroidism and 648 (2.6%) had subclinical hyperthyroidism. In age- and sex-adjusted analyses, risks of heart failure events were increased with both higher and lower TSH levels (P for quadratic pattern <0.01); the hazard ratio was 1.01 (95% confidence interval, 0.81-1.26) for TSH of 4.5 to 6.9 mIU/L, 1.65 (95% confidence interval, 0.84-3.23) for TSH of 7.0 to 9.9 mIU/L, 1.86 (95% confidence interval, 1.27-2.72) for TSH of 10.0 to 19.9 mIU/L (P for trend <0.01) and 1.31 (95% confidence interval, 0.88-1.95) for TSH of 0.10 to 0.44 mIU/L and 1.94 (95% confidence interval, 1.01-3.72) for TSH <0.10 mIU/L (P for trend=0.047). Risks remained similar after adjustment for cardiovascular risk factors. CONCLUSION: Risks of heart failure events were increased with both higher and lower TSH levels, particularly for TSH ≥10 and <0.10 mIU/L.

Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

Simulation methods with extended stability for stiff biochemical kinetics

Background: With increasing computer power, simulating the dynamics of complex systems in chemistry and biology is becoming increasingly routine. The modelling of individual reactions in (bio)chemical systems involves a large number of random events that can be simulated by the stochastic simulation algorithm (SSA). The key quantity is the step size, or waiting time, τ, whose value inversely depends on the size of the propensities of the different channel reactions and which needs to be re-evaluated after every firing event. Such a discrete event simulation may be extremely expensive, in particular for stiff systems where τ can be very short due to the fast kinetics of some of the channel reactions. Several alternative methods have been put forward to increase the integration step size. The so-called τ-leap approach takes a larger step size by allowing all the reactions to fire, from a Poisson or Binomial distribution, within that step. Although the expected value for the different species in the reactive system is maintained with respect to more precise methods, the variance at steady state can suffer from large errors as τ grows. Results: In this paper we extend Poisson τ-leap methods to a general class of Runge-Kutta (RK) τ-leap methods. We show that with the proper selection of the coefficients, the variance of the extended τ-leap can be well-behaved, leading to significantly larger step sizes.Conclusions: The benefit of adapting the extended method to the use of RK frameworks is clear in terms of speed of calculation, as the number of evaluations of the Poisson distribution is still one set per time step, as in the original τ-leap method. The approach paves the way to explore new multiscale methods to simulate (bio)chemical systems.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Distribution of events of positive selection and population differentiation in a metabolic pathway: the case of asparagine N-glycosylation

Background: Asparagine N-Glycosylation is one of the most important forms of protein post-translational modification in eukaryotes. This metabolic pathway can be subdivided into two parts: an upstream sub-pathway required for achieving proper folding for most of the proteins synthesized in the secretory pathway, and a downstream sub-pathway required to give variability to trans-membrane proteins, and involved in adaptation to the environment andinnate immunity. Here we analyze the nucleotide variability of the genes of this pathway in human populations, identifying which genes show greater population differentiation and which genes show signatures of recent positive selection. We also compare how these signals are distributed between the upstream and the downstream parts of the pathway, with the aim of exploring how forces of population differentiation and positive selection vary among genes involved in the same metabolic pathway but subject to different functional constraints. Results:Our results show that genes in the downstream part of the pathway are more likely to show a signature of population differentiation, while events of positive selection are equally distributed among the two parts of the pathway. Moreover, events of positive selection arefrequent on genes that are known to be at bifurcation points, and that are identified as beingin key position by a network-level analysis such as MGAT3 and GCS1.Conclusions: These findings indicate that the upstream part of the Asparagine N-Glycosylation pathway has lower diversity among populations, while the downstream part is freer to tolerate diversity among populations. Moreover, the distribution of signatures of population differentiation and positive selection can change between parts of a pathway, especially between parts that are exposed to different functional constraints. Our results support the hypothesis that genes involved in constitutive processes can be expected to show lower population differentiation, while genes involved in traits related to the environment should show higher variability. Taken together, this work broadens our knowledge on how events of population differentiation and of positive selection are distributed among different parts of a metabolic pathway.

Involvement of the extracellular signal-regulated kinase cascade for cocaine-rewarding properties

A central feature of drugs of abuse is to induce gene expression in discrete brain structures that are critically involved in behavioral responses related to addictive processes. Although extracellular signal-regulated kinase (ERK) has been implicated in several neurobiological processes, including neuronal plasticity, its role in drug addiction remains poorly understood. This study was designed to analyze the activation of ERK by cocaine, its involvement in cocaine-induced early and long-term behavioral effects, as well as in gene expression. We show, by immunocytochemistry, that acute cocaine administration activates ERK throughout the striatum, rapidly but transiently. This activation was blocked when SCH 23390 [a specific dopamine (DA)-D1 antagonist] but not raclopride (a DA-D2 antagonist) was injected before cocaine. Glutamate receptors of NMDA subtypes also participated in ERK activation, as shown after injection of the NMDA receptor antagonist MK 801. The systemic injection of SL327, a selective inhibitor of the ERK kinase MEK, before cocaine, abolished the cocaine-induced ERK activation and decreased cocaine-induced hyperlocomotion, indicating a role of this pathway in events underlying early behavioral responses. Moreover, the rewarding effects of cocaine were abolished by SL327 in the place-conditioning paradigm. Because SL327 antagonized cocaine-induced c-fos expression and Elk-1 hyperphosphorylation, we suggest that the ERK intracellular signaling cascade is also involved in the prime burst of gene expression underlying long-term behavioral changes induced by cocaine. Altogether, these results reveal a new mechanism to explain behavioral responses of cocaine related to its addictive properties.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

Bridging the Faraoni and Selli oceanic anoxic events: late Hauterivian to early Aptian dysaerobic to anaerobic phases in the Tethys

A detailed geochemical analysis was performed on the upper part of the Maiolica Formation in the Breggia (southern Switzerland) and Capriolo sections (northern Italy). The analysed sediments consist of well-bedded, partly siliceous, pelagic carbonate, which lodges numerous thin, dark and organic-rich layers. Stable-isotope, phosphorus, organic-carbon and a suite of redox-sensitive trace-element contents (RSTE: Mo, U, Co, V and As) were measured. The RSTE pattern and C-org:P-tot ratios indicate that most organic-rich layers were deposited under dysaerobic rather than anaerobic conditions and that latter conditions were likely restricted to short intervals in the latest Hauterivian, the early Barremian and the pre-Selli early Aptian. Correlations are both possible with organic-rich intervals in central Italy (the Gorgo a Cerbara section) and the Boreal Lower Saxony Basin, as well as with the facies and drowning pattern in the Helvetic segment of the northern Tethyan carbonate platform. Our data and correlations suggest that the latest Hauterivian witnessed the progressive installation of dysaerobic conditions in the Tethys, which went along with the onset in sediment condensation, phosphogenesis and platform drowning on the northern Tethyan margin, and which culminated in the Faraoni anoxic episode. This episode is followed by further episodes of dysaerobic conditions in the Tethys and the Lower Saxony Basin, which became more frequent and progressively stronger in the late early Barremian. Platform drowning persisted and did not halt before the latest early Barremian. The late Barremian witnessed diminishing frequencies and intensities in dysaerobic conditions, which went along with the progressive installation of the Urgonian carbonate platform. Near the Barremian-Aptian boundary, the increasing density in dysaerobic episodes in the Tethyan and Lower Saxony Basins is paralleled by a change towards heterozoan carbonate production on the northern Tethyan shelf. The following return to more oxygenated conditions is correlated with the second phase of Urgonian platform growth and the period immediately preceding and corresponding to the Selli anoxic episode is characterised by renewed platform drowning and the change to heterozoan carbonate production. Changes towards more humid climate conditions were the likely cause for the repetitive installation of dys- to anaerobic conditions in the Tethyan and Boreal basins and the accompanying changes in the evolution of the carbonate platform towards heterozoan carbonate-producing ecosystems and platform drowning.

The Toarcian in the Subbetic basin (southern Spain): Bio-events (ammonite and calcareous nannofossils) and carbon-isotope stratigraphy

A detailed carbon-isotope stratigraphic study for the uppermost Pliensbachian lowermost Aalenian interval in the Median Subbetic palaeogeographic domain (External zones of the Betic Cordillera, southern Spain) has been carried out. During the Early Jurassic, the Median Subbetic, which represents a typical basin of the Hispanic Corridor connecting the Tethys and the Eastern Pacific, was located in the westernmost Tethys. The analyzed sections encompass the entire Toarcian stage as represented in the southern Iberian palaeomargin. Rocks are mainly rhythmic sequences of grey marls and marly limestones containing a rich ammonite fauna, nannofossils, and benthic foraminifers-all these provide an accurate biostratigraphic control. The lower and upper Toarcian boundaries are well represented in some of these sections and therefore represent optimal sites to link the carbon-isotope curves to ammonite zones, and to nannofossil events. delta C-13 values of bulk carbonates from the different localities of the Subbetic basin have similar variations from the uppermost Pliensbachian to the lowermost Aalenian, suggesting changes in the original DIC carbon isotope composition along the Hispanic corridor. The transition from Pliensbachian to Toarcian is marked by increasing delta C-13 values from similar to 12 to 2.0 parts per thousand, interrupted in the Serpentinum Zone by a negative shift concomitant with the Toarcian oceanic anoxic event (T-OAE), with the major ammonite extinction event of the Toarcian, and an important turnover of calcareous nannoplankton. The negative shift observed in the Serpentinum Zone confirms the global perturbation of the carbon cycling documented along the Tethys and the palaeo-Pacific in organic material and in marine carbonates. However, the amplitude of the negative excursion (similar to - 1.5 parts per thousand) is not compatible with an isotopic homogeneous seawater DIC and/or CO2 atmospheric reservoirs. The interval from the middle to the top of the Toarcian delta C-13 shows relatively constant values, minor ammonite turnovers, and is associated with increasing diversity of calcareous nannoplankton. (c) 2012 Elsevier B.V. All rights reserved.

The risk of subsequent cancer and arterial cardiovascular events in patients with superficial vein thrombosis in the legs.

Although it has been clearly demonstrated that venous thromboembolism is associated with an increased risk of subsequent overt cancer and arterial cardiovascular events in comparison with control populations, whether this association also applies to patients with isolated (ie, without concomitant involvement of the deep vein system) superficial vein thrombosis (SVT) in the legs is unknown. In 737 consecutive patients with isolated SVT not involving the sapheno-femoral junction, we conducted a retrospective investigation to assess the rate of cancer and that of arterial cardiovascular events occurring during follow-up. The event rates were compared with those occurring in 1438 controls having comparable characteristics. Both cases and controls were followed-up for an average period of 26 ± 8 months (range, 3-45). Malignancy was diagnosed in 26 cases (3.5%) and 56 controls (3.9%), leading to a hazard ratio of 0.86 (95% confidence interval, 0.55%-1.35%). Arterial cardiovascular events occurred in 32 cases (4.3%) and 63 controls (4.4%), leading to a hazard ratio of 0.97 (95% confidence interval, 0.63%-1.50%). We conclude that the occurrence of isolated SVT in the legs does not place patients at an increased risk of malignancies or arterial cardiovascular events. Whether this conclusion also applies to patients whose thrombosis involves the sapheno-femoral junction remains to be demonstrated.