Past seminars

2020

2019

  • Lundi 2 dĂ©cembre 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Aline Marguet
    (INRIA, IBIS)
    TBA
  • Lundi 18 novembre 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Andrea Rau
    (INRA, GABI)
    TBA
  • Lundi 21 octobre 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Julie Aubert
    (AgroParisTech-INRA, UMR 518 MIA)
    Unraveling biotic interactions determining soil microbial community assembly and functioning

    Joint work with S. Romdhane, A. Spor, S. Ouadah, L. Philippot

    Microbial communities play important roles in all ecosystems and yet, a sound understanding of the ecological processes governing the assembly of these communities in the environment is missing. To address the role of biotic interactions in assembly and functioning of the soil microbiota, we used a top down manipulation approach based on the removal of various populations in a natural microbial community. Suspensions of the soil microbiota were subjected to various biocidal and filtration treatments before being inoculated into the same sterilized soil. We hypothesized that if biotic interactions is an important force shaping the assembly of microbiota, removal of microbial groups should largely affect the fitness of the remaining ones during soil recolonization. We show that nearly 50 % of the dominant bacterial taxa were subjected to competitive interactions, underlining the importance of biotic interactions in the assembly of microbial communities in soil. Moreover, evidence for competitive exclusion between members of Bacillales and Proteobacteriales suggests that potential general rules of microbial community assembly can be identified. Reassembly after removal resulted in greater changes in activities related to N- than to C-cycling indicating functional differences. Our approach can provide a new avenue to study microbial interactions in complex ecosystems.

  • Lundi 7 octobre 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Blaise Hanczar
    (Université d'Evry, IBISC)
    Apprentissage profond pour la prĂ©diction de phĂ©notypes Ă  partir de donnĂ©es d’expression

    L’apprentissage profond (deep learning) est une mĂ©thode d’intelligence artificielle qui a rĂ©cemment permis de faire d’impressionnants progrĂšs dans de nombreux domaines scientifiques. Parmi ses applications les plus prometteuses se trouvent la mĂ©decine personnalisĂ©e et l’analyse des donnĂ©es «omiques».  Aujourd’hui, des efforts croissants sont dĂ©ployĂ©s pour utiliser l’apprentissage profond sur des donnĂ©es gĂ©nomiques, transcriptomiques, protĂ©omiques et mĂ©ta-gĂ©nomiques afin de mieux caractĂ©riser les patients. AprĂšs avoir exposĂ© briĂšvement le principe des rĂ©seaux de neurones, je prĂ©senterai nos travaux sur la prĂ©diction de phĂ©notypes Ă  partir de donnĂ©es d’expression de gĂšnes. Dans cette tĂąche, les deux principaux verrous scientifiques sont la petite taille de la base d’apprentissage et l’interprĂ©tation du rĂ©seau. Les  mĂ©thodes basĂ©es sur l’apprentissage par transfert et l’apprentissage semi-supervisĂ© seront utilisĂ©es afin de pallier au premier problĂšme. Pour l’interprĂ©tation, nous proposons de dĂ©composer le calcul de chaque prĂ©diction dans le rĂ©seau de neurones afin d'associer Ă  chaque neurone de la connaissance biologique.

  • Lundi 23 septembre 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Timothée Tabouy
    (AgroParisTech-INRA, UMR 518 MIA)
    Impact de l’échantillonnage sur l’infĂ©rence de structures dans les rĂ©seaux. Application aux rĂ©seaux d’échanges de graines et Ă  l’écologie

    Dans cette thĂšse nous nous intĂ©ressons Ă  l’étude du modĂšle Ă  bloc stochastique (SBM) en prĂ©sence de donnĂ©es manquantes. Nous proposons une classification des donnĂ©es manquantes en deux catĂ©gories Missing At Random et Not Missing At Random pour les modĂšles Ă  variables latentes suivant le modĂšle dĂ©crit par D. Rubin. De plus, nous nous sommes attachĂ©s Ă  dĂ©crire plusieurs stratĂ©gies d’échantillonnages de rĂ©seau et leurs lois. L’infĂ©rence des modĂšles de SBM avec donnĂ©es manquantes est faite par l’intermĂ©diaire d’une adaptation de l’algorithme EM : l’EM avec approximation variationnelle. L’identifiabilitĂ© de plusieurs des SBM avec donnĂ©es manquantes a pu ĂȘtre dĂ©montrĂ©e ainsi que la consistance et la normalitĂ© asymptotique des estimateurs du maximum de vraisemblance et des estimateurs avec approximation variationnelle dans le cas oĂč chaque dyade (paire de nƓuds) est Ă©chantillonnĂ©e indĂ©pendamment et avec mĂȘme probabilitĂ©. Nous nous sommes aussi intĂ©ressĂ©s aux modĂšles de SBM avec covariables, Ă  leurs infĂ©rence en prĂ©sence de donnĂ©es manquantes et comment procĂ©der quand les covariables ne sont pas disponibles pour conduire l’infĂ©rence. Finalement, toutes nos mĂ©thodes ont Ă©tĂ© implĂ©mentĂ©es dans un package R disponible sur le CRAN. Une documentation complĂšte sur l’utilisation de ce package a Ă©tĂ© Ă©crite en complĂ©ment.

  • Lundi 9 septembre 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Aurélien Latouche
    (CNAM, MSDMA)
    Semiparametric approach for covariate-specific time dependent ROC curves for correlated survival data / Validation d'un marqueur pronostique à partir de données de survie en grappe

    Joint work with Alessandra Meddis and Paul Blanche

    1 Introduction
    Considerable research has focused on the development of new biomarkers. The first step in developing a clinically useful biomarker is to identify its ability in discriminating patients at high risk of dying within the next t-years (e.g. 5-years) from those who will not. The standard methodology to quantify the discrimination performance of a biomarker, with right censored data, is to estimate time dependent ROC curves, ROC(t). In presence of clustered failure times, the common strategy is to ignore heterogeneity in the phase of evaluation of the performance of a candidate biomarker, but to confirm its discriminatory capacity, it is important to account for heterogeneity while adjusting for clinical covariates. The usefulness of our approach is illustrated on our motivating example, which consists in the first meta-analysis on individual data of more than 2000 patients from 15 centers with non metastatic breast cancer. Its objective was to quantify the clinical usefulness of circulating tumor cells (CTCs) count as a prognostic marker of survival.
    2 Methodology
    ROC(t) allows to study the capacity of a biomarker Y to discriminate between patients who experience event prior time t (cumulative cases) from those who do not up to time t (dynamic controls). The current methodol- ogy does not account for heterogeneity while estimating ROC(t). In this work, we fill this gap by proposing an extension to clustered data of the Song & Zhou method (Statistica Sinica, 2008). To estimate the covariate- specific time dependent ROC curve we consider a joint model: (i) shared frailty model which links the covariates and the biomarker to the time-to-event, (ii) location scale model to link the covariates to the biomarker. We evaluate the performance of the proposed method in a simulation study. We demonstrate an application of the estimator to data derived from a meta-analysis on individual patient data with non metastatic breast cancer where the goal is to understand the clinical usefulness of CTCs count for this scenario. In particular, we estimate the covariate-specific ROC curves that quantify the discrimination performance of CTCs count within subgroups of patients having the same tumor stage at time of diagnosis, since subjects with inflammatory tumor show a higher number of CTCs and a poorer prognosis. A bootstrap method is proposed for calculating confidence intervals.
    3 Results
    The estimator is computationally simple and the simulation results highlighted the robustness of the method at varying of censoring with negligible bias (≈ 10−3). Moreover, we provide the results for the motivating example with the time dependent ROC curves and respective AUCs for different tumor stage. The wide confidence intervals highlighted that having inflammatory tumor does not influence the discrimination of the CTCs count.
    4 Conclusions
    In presence of clustered failure times it is important to take into account heterogeneity. In fact, the introduction of a random effect (frailty) is needed to estimate the performance of the biomarker in the general population. In this scenario, the covariate-specific time dependent ROC curve can be easily estimated with the proposed approach.

  • Mardi 25 juin 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Mathieu Roche
    Analyse spatiale et thématique par des méthodes de fouille de textes

    Les dynamiques spatiales et thĂ©matiques font l'objet d'un intĂ©rĂȘt croissant notamment dans les milieux Ă  enjeux et de tension (Ă©pidĂ©miologie, sĂ©curitĂ© alimentaire, etc.). L’objectif de nos travaux est de proposer un cadre mĂ©thodologique permettant l’apprĂ©hension de ces phĂ©nomĂšnes Ă  partir de donnĂ©es textuelles hĂ©tĂ©rogĂšnes. Les analyses produites s'appuient, en gĂ©nĂ©ral, sur trois types d'informations (thĂ©matiques, spatiales et temporelles) qui sont extraites et exploitĂ©es par des mĂ©thodes de fouille de donnĂ©es et de traitement automatique du langage naturel (TALN). Les contributions mĂ©thodologiques seront prĂ©sentĂ©es Ă  travers diffĂ©rents projets pluridisciplinaires et le dĂ©ploiement d'outils dĂ©diĂ©s, en particulier le systĂšme PADI-Web propre Ă  la veille automatique en Ă©pidĂ©miologie animale.

  • Lundi 17 juin 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Arnaud Gloaguen
    (CentraleSupelec, L2S)
    Joint Matrix/Tensor Factorization with MGCCA

    Regularized Generalized Canonical Correlation Analysis (RGCCA) is a general multiblock data analysis framework that encompasses several important multivariate analysis methods such as principal component analysis, partial least squares regression and several versions of generalized canonical correlation analysis. In this paper, we extend RGCCA to the case where at least one block has a tensor structure. This method is called Multiway Generalized Canonical Correlation Analysis (MGCCA). Convergence properties of the MGCCA algorithm are studied and computation of higher-level components are discussed. The usefulness of MGCCA is shown on simulation and on the analysis of a cognitive study in human infants using high-density electro-encephalography (EEG).

  • Lundi 6 mai 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Pierre Alquier
    (ENSAE)
    Generalization bounds for variational inference.

    Bayesian inference provides an attractive learning framework to analyze data, and to sequentially update knowledge on streaming data. Unfortunately, exact Bayesian inference is rarely feasible in practice and approximation methods are usually employed, but do such methods preserve the generalization properties of Bayesian inference? In this talk, I will show that it is indeed the case for some variational inference (VI) algorithms. First, I will show generalization bounds for estimation in the batch setting. These results can be seen as extensions of the "concentration of the posterior" theorems to variational approximations of the posterior. I will then focus on the sequential case (streaming data). I will propose various online VI algorithms and derive generalization bounds. In this case, our theoretical result relies on the convexity of the variational objective, but we argue that our result should hold more generally and present empirical evidence in support of this.

    Joint works with James Ridgway (https://arxiv.org/abs/1706.09293), Badr-Eddine Chérief-Abdellatif (https://projecteuclid.org/euclid.ejs/1537344604) and Emti Khan (will appear on arXiv in the next few days)

  • Lundi 15 avril 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Frédérique Clément
    (INRIA Saclay)
    ModĂ©lisation multi-Ă©chelles en biologie du dĂ©veloppement et de la reproduction : une approche “middle-out” basĂ©e sur la dynamique cellulaire

    Dans cet exposé, nous présenterons des approches de modélisation multi-échelles, basées sur la
    description des dynamiques cellulaires à l’Ɠuvre dans un processus relevant à la fois de la biologie du
    développement et de la biologie de la reproduction, la folliculogenÚse ovarienne.
    Le formalisme sous-jacent est celui des dynamiques de populations structurĂ©es, considĂ©rĂ©es soit d’un
    point de vue dĂ©terministe (Ă©quations de transport non conservatives), soit d’un point de vue stochastique
    (chaĂźnes de Markov Ă  temps continu). Ce formalisme permet d’établir des liens mĂ©canistes
    entre le contrÎle du devenir cellulaire (prolifération, différenciation apoptose) et celui du destin fol-
    liculaire (ovulation ou dĂ©gĂ©nĂ©rescence), ou d’étudier comment les interactions entre l’ovocyte et les
    cellules somatiques environnantes façonnent la morphogenÚse du follicule.
    Nous illustrerons ces principes, d’une part, sur le dĂ©veloppement terminal qui est intĂ©grĂ© dans la
    dynamique du cycle ovarien rĂ©sultant des interactions hormonales au sein de l’axe hypothalamo-
    hypophyso-gonadique, et, d’autre part, sur une Ă©tape du dĂ©veloppement prĂ©coce pendant laquelle
    le follicule ovarien se stratifie en plusieurs couches cellulaires concentriques à l’ovocyte.
    Dans le premier cas, nous introduirons un systĂšme d’équations de transport faiblement couplĂ©es,
    permettant de rendre compte du processus de sélection des follicules ovulatoires. Dans le deuxiÚme cas,
    nous décrirons en parallÚle un processus de branchement multi-type et son équivalent déterministe.

  • Lundi 1er avril 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Pierre Gloaguen
    (AgroParisTech, MIAP)
    La diffusion de Langevin comme modÚle pour le déplacement animal et la sélection d'habitat

    The utilisation distribution describes the relative probability of use of a spatial unit by an
    animal. It is natural to think of it as the long-term consequence of the animal's short-term
    movement decisions: it is the accumulation of small displacements which, over time, gives
    rise to global patterns of space use. However, most utilisation distribution models either
    ignore the underlying movement, assuming the independence of observed locations, or are
    based on simplistic Brownian motion movement rules.
    We introduce a new continuous-time model of animal movement, based on the Langevin
    diffusion. This stochastic process has an explicit stationary distribution, conceptually anal-
    ogous to the idea of the utilisation distribution, and thus provides an intuitive framework to
    integrate movement and space use.
    We model the stationary (utilisation) distribution with a resource selection function to link
    the movement to spatial covariates, and allow inference into habitat selection.
    Standard approximation techniques can be used to derive the pseudo-likelihood of the Langevin
    diffusion movement model, and thus to estimate habitat preference and movement parame-
    ters from tracking data. We investigate the performance of the method on simulated data,
    and discuss its sensitivity to the time scale of the sampling. We present an example of its
    application to tracking data of Stellar sea lions (fEumetopias jubatusg).
    Due to its continuous-time formulation, this method can be applied to irregular teleme-
    try data. It provides a rigorous framework to estimate long-term habitat selection from
    correlated movement data.

    Joint work with Marie-Pierre Etienne (AgroCampus Ouest) and Theo Michelot (University of Shefield)

  • Lundi 18 mars 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Amandine Cornille
    (INRA - UniversitĂ© Paris-Sud - CNRS - AgroParisTech, GĂ©nĂ©tique Quantitative et Évolution)
    Ecological drivers and genomics bases of the aphid-apple-bacteria interaction in the context of domestication

    Amandine’s research aims at unravelling the evolutionary forces underlying the adaptation of populations and the emergence of new species, using population genetics and genomics approaches combined with field experiments. She has studied both short and long evolutionary timescales (domestication and wild species divergence, respectively) to investigate the ecological factors and genomic processes underlying adaptation. Her research is now also moving toward the study of plant biotic interactions, which are fundamental to our understanding of biodiversity and adaptation; research over the past decades has indeed clearly demonstrated that interactions between species are major drivers of adaptive evolution. Amandine will present the main results of her past research on the apple domestication and the link with her starting project on the impact of apple domestication on the adaptation of a major apple aphid pest Dysaphis plantaginea (the rosy apple aphid) to the cultivated apple (Malus domestica).

  • Lundi 4 mars 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Ralf Steuer
    (Humboldt, UniversitÀt zu Berlin, Institut of Theoretical Biology)
    Cellular resource allocation and the architecture of phototrophic growth

    The potential offered by cyanobacteria will undoubtedly play a major role in mastering the challenges of the 21st century – from securing global food supply to the synthesis of renewable raw materials. As yet, however, fundamental questions regarding the metabolic principles of cyanobacterial phototrophic growth are not resolved: How are metabolic, photosynthetic, and ribosomal proteins optimally partitioned during phototrophic growth? What is the highest growth rate a cyanobacterium can attain?
    The presentation will describe our recent efforts to answer these questions using experimental and constraint-based computational models: we describe phototrophic growth as a cellular resource allocation problem and estimate the costs and benefits of all metabolic constituents of a cyanobacterial cell. Of particular interest are the cellular organization that enables fast phototrophic growth and the corresponding intracellular limits on growth rates. The model-derived resource allocation is in good agreement with experimental findings. I will provide an outlook how such analyses has implications for ecology and models of global biogeochemical cycles.

    References:
    [1] Zavƙel T(*), Faizi M, Loureiro C, Poschmann G, StĂŒhler K, Sinetova M, Zorina A, Steuer R(*), ČervenĂœ J (2019) Quantitative insights into the cyanobacterial cell economy. eLife  [* corresponding authors]
    [2] Faizi M, Zavrel T, Loureiro C, Cerveny J, Steuer R (2018) A model of optimal protein allocation during phototrophic growth. Biosystems 166, 26-36.
    [3] Reimers AM, Knoop H, Bockmayr A, Steuer R (2017) Cellular trade-offs and optimal resource allocation during cyanobacterial diurnal growth. Proc Natl Acad Sci U S A. Pii: 201617508.
    [4] Westermark S and Steuer R (2016) Toward multiscale models of cyanobacterial growth: a modular approach. Front. Bioeng. Biotechnol. 4:95.
    [5] Knoop H, Gruendel M, Zilliges Y, Lehmann R, Hoffmann S, Lockau W, Steuer R (2013) Flux balance analysis of cyanobacterial metabolism: The metabolic network of Synechocystis sp. PCC 6803. PLoS Comput Biol 9(6): e1003081.

  • Lundi 18 fĂ©vrier 2019 - - Salle de rĂ©union 142, bĂątiment 210
    SĂ©bastien Picault
    (INRA-Oniris BIOEPAR Nantes et Univ. Lille-CNRS CRIStAL)
    L'Intelligence Artificielle pour l'épidémiologie prédictive : faciliter la co-construction et la révision des modÚles

    La modĂ©lisation mĂ©caniste permet de mieux comprendre la propagation de pathogĂšnes dans et entre populations, d’en prĂ©dire la dynamique dans des situations contrastĂ©es, et d'identifier les mesures de gestion les plus pertinentes. Les enjeux sanitaires actuels demandent toutefois une diversitĂ© d'expertises qui accroĂźt la complexitĂ© des modĂšles et rend coĂ»teuse leur conception, au risque de rĂ©duire leur fiabilitĂ©, maintenabilitĂ© et rĂ©utilisabilitĂ©.

    Mes travaux à BIOEPAR (Nantes) visent à élaborer un nouveau cadre de conception de modÚles mécanistes en épidémiologie via des recherches en Intelligence Artificielle. L'utilisation de systÚmes multi-agents multi-niveaux permet d'encapsuler les paradigmes classiques de modélisation (modÚles à compartiments vs. centrés individus) à différentes échelles (individu, population, métapopulation) dans un formalisme homogÚne et modulaire, facilitant le changement de
    reprĂ©sentation pour un mĂȘme systĂšme. SimultanĂ©ment, l'Ă©laboration d'un langage domaine-spĂ©cifique (DSL) pour l’épidĂ©miologie permet d’expliciter les hypothĂšses, paramĂštres, processus mobilisĂ©s, et donc de renforcer l'implication des experts non modĂ©lisateurs dĂšs la conception des modĂšles et tout au long du processus de dĂ©veloppement.

    La combinaison d’une approche multi-agents et d'un DSL a ainsi conduit au framework EMULSION, dĂ©veloppĂ© depuis fin 2016 Ă  BIOEPAR. Je montrerai
    une application d'EMULSION Ă  l'Ă©tude d'une zoonose endĂ©mique du bĂ©tail (la fiĂšvre Q), pour comparer diverses hypothĂšses (Ă©tats de santĂ© pertinents, forme de la fonction de dispersion par le vent, etc.), eu Ă©gard aux donnĂ©es disponibles, et ainsi mieux comprendre la contribution des mouvements commerciaux vs. des conditions environnementales dans la propagation du pathogĂšne Ă  l'Ă©chelle d’un territoire.

    Enfin, je prĂ©senterai comment ces travaux ouvrent la voie Ă  une production automatisĂ©e d’outils d'aide Ă  la dĂ©cision dĂ©diĂ©s aux gestionnaires de la santĂ© (techniciens conseils, vĂ©tĂ©rinaires, dĂ©cideurs publics).

  • Lundi 4 fĂ©vrier 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Vincent Miele
    (UCB Lyon1, UMR CNRS 5558 LBBE)
    Exploring multiplex ecological networks

    Within an ecosystem, species interact with each other in many different ways, including predation, competition, and facilitation, and this can be modelled as a network of multiple interaction types (called "multiplex" network). The variety of interaction types that link species to each other has long been recognized but has rarely been synthesized for entire multi-species ecosystems. In a first part, I will present our recent study [1] on a unique marine ecological network that integrates thousands of trophic and non-trophic interactions. Using network statistical modeling, we showed that the diverse interaction types have a suprinsingly simple non-random organization. In a second part, I will present our recent study [2] where we investigated the role of the diversity of interaction types per se for species diversity and the functioning of ecological communities, using a bio-energetic consumer-resource dynamical model in which we incorporated the non-trophic interactions.

    En collaboration avec Sonia KĂ©fi (CNRS, Institut des sciences de l'Ă©volution, Montpellier).

    [1] KĂ©fi et al, PLoS Biol (2016), https://doi.org/10.1371/journal.pbio.1002527
    [2] Miele et al, biorxiv (2018), https://doi.org/10.1101/411249

  • Lundi 21 janvier 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Paul Bastide
    (KU Leuven)
    A flexible Bayesian framework to study viral trait evolution

    During the course of an outbreak or epidemic, many viral pathogens are known to evolve rapidly, leaving imprint of the pattern of spread in their genomes. Uncovering the molecular footprint of this transmission process is a key goal of phylodynamic inference. Relatively less focus has been put on the evolution of quantitative phenotypic traits of viruses. Traits such as geographical location or virulence can be studied using phylogenetic Comparative Methods (PCMs) that account for a shared evolutionary history among the set of non-independent samples. Conditioning on such an history, the observed traits can be seen as the result of a stochastic process running on the branches of a phylogenetic tree. The Ornstein-Uhlenbeck (OU) process is often used to model stabilizing selection toward an optimal trait value. For a multivariate trait, the dynamics of the trajectory is controlled by a selection strength matrix, that is only constrained to have positive eigenvalues. Depending on the form of this matrix, the OU can have a variety of behaviors, and is hence suited to model various biological processes.

    We propose a Bayesian inference framework for the study of this flexible model. Using a Markov Chain Monte Carlo (MCMC) based method, one critical aspect is to be able to sample uniformly in the space of constrained matrices, both for the selection strength and the variance matrix, in a context where traditional Gibbs sampling cannot be used. This can be done using a smooth transformation that maps the parameters to an unconstrained space. We investigated the use of two such maps, along with adequate prior distributions. MCMC methods also rely on multiple likelihood evaluations, at each step of the chain. Exploiting the tree structure, we studied a fast and flexible algorithm to compute both the likelihood and its gradient for a wide class of processes, that contains but is not limited to the OU. This makes it possible to use efficient sampling methods, such as the Hamiltonian Monte Carlo (HMC).

    We implemented the new framework in BEAST, a widely used and flexible phylodynamics software. This allows us to leverage on the many other tools of the BEAST ecosystem, such as the phylogenetic factor model, that can be used to model extra-environmental variation, or the marginal likelihood estimation for model selection. It also offers us the possibility to integrate the results over the space of all probable trees, in an integrated analysis that directly starts from the genomic sequences, instead of relying on a fixed tree. We illustrate the use of this framework for the study of the heritability of virulence of the human immunodeficiency virus (HIV), a question that has attracted a lot of attention recently, and for which model choice is a recognized critical aspect.

  • Lundi 7 janvier 2019 - - Salle de rĂ©union 142, bĂątiment 210
    Nathanael Hozé
    (Institut Pasteur, Unité de Modélisation mathématique des maladies infectieuses)
    Assessing the level of virus circulation in the context of high cross-reactivity: the case of chikungunya and mayaro fever in French Guiana

    Dengue, chikungunya and Zika have recently caused large outbreaks in the Americas. Besides those, other arboviruses may have been overlooked. For example, Mayaro virus (MAYV), which was first identified in 1954 in Trinidad, is very rarely detected by Sentinel surveillance in French Guiana (FG). However, we cannot rule out that some MAYV infections may have been misdiagnosed because Mayaro fever shares symptoms with other arboviral infections and there is cross-reactivity between MAYV and chikungunya virus (CHIKV). We propose a statistical approach to assess the circulation of both MAYV and CHIKV in FG using a multiplexed cross-sectional serological survey that includes 2,697 individuals sampled over the FG territory. We account explicitly for the cross-reactivity between MAYV and CHIKV. We find that 41% (95% CI: 35%, 48%) of MAYV cases might be misdiagnosed as CHIKV positive with serology assays. Our framework also allows testing and exploring different scenarios of virus circulation for each region of the territory. We find that MAYV has been circulating for years in the southern part of FG, affecting mostly adult males, and we estimate seroprevalence as high as 27% (95% CI: 19%, 35%) in some regions. On the opposite, a CHIKV epidemic hit FG in 2014, was more present in the northern, coastal regions and infected women more than men. 

2018

  • Lundi 17 dĂ©cembre 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Sarah Lemler
    (CentraleSupelec, MICS)
    Estimation on a jump diffusion process with jumps driven by a Hawkes process

    In neurosciences, two kinds of data are available: a continuous signal coming from the action potential of one single neuron along time, and several discrete signals which are the spike train of several neurons (the occurring times when the action potential overcomes a threshold). Usually these two kinds of data are investigated separately using either diffusion processes to describe the dynamic of the action potential or counting processes to describe the graph of connectivity between neurons (see [3]). We propose a new model which allows to deal with both type of signals in a single procedure to take into accounts the whole information. This new approach can be applied to neurosciences issues but also to other fields as in finance to model financial assets. Our model is a jump-diffusion model driven by a multidimensional Hawkes process (see [1]).  Assuming that we have high frequencies observations of the process together with observations of the multidimensionnal Hawkes process, we focus first on the nonparametric estimation of the drift function, when the other coefficients are assumed to be known (see [2]) and then we consider a simultaneous estimation of both parameters the jump parameter and the drift, assuming that $\sigma$ is known. The proposed method is based on a mean squared regression approach. The novelty here lies in the contrast function which takes into account the multidimensional Hawkes process. We will also discuss about the existence of a stationary measure for our process and ergodicity results that we need to obtain oracle inequalities for our estimators. Finally, a numerical study illustrates our purpose.

    [1] P. Brémaud and L. Massoulié. Stability of nonlinear Hawkes processes. The Annals of Probability, 24, 1996, 1563-1588
    [2] F. Comte and V. Genon-Catalot and Y. Rozenholc. Penalized nonparametric mean square estimation of the coefficients
    of diffusion processes. Bernoulli, 13, 2007, 514-543
    [3] P. Reynaud-Bouret and V. Rivoirard and C. Tuleau-Malot. Inference of functional connectivity in neurosciences via Hawkes processes. Global Conference on Signal and Information Processing (GlobalSIP), 2013 IEEE,  317-320
  • Lundi 3 dĂ©cembre 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Radu Stoica
    (Université de Lorraine, Institut Elie Cartan de Lorraine, Nancy)
    ABC Shadow algorithm: a tool for statistical analysis of spatial patterns

    This talk presents an original ABC algorithm, {\it ABC Shadow}, that can be applied to sample posterior densities that are continuously differentiable. The proposed algorithm solves the main condition to be fulfilled by any ABC algorithm, in order to be useful in practice. This condition requires enough samples in the parameter space region, induced by the observed statistics. The algorithm is tuned on the posterior of a Gaussian model which is entirely known, and then it is applied for the statistical analysis of several spatial patterns. These patterns are issued or assumed to be outcomes of point processes. The considered models are: Strauss, Candy and area-interaction. Results on real data, representing the galaxies distribution in our Universe, are also shown.

  • Lundi 19 novembre 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Amandine VĂ©ber
    (Ecole Polytechnique, CMAP)
    Le processus Lambda-Fleming-Viot spatial : un modĂšle d'Ă©volution en espace continu

    Ce processus a Ă©tĂ© introduit en 2008 par N. Barton (IST Austria) et A. Etheridge (U. Oxford) pour modĂ©liser l'Ă©volution de la diversitĂ© gĂ©nĂ©tique au sein d'une population ayant une structure spatiale continue (un champ, une forĂȘt, les Alpes, ...). Dans ce modĂšle, la rĂ©gulation locale de la taille de population est assurĂ©e en encodant les reproductions non par des "horloges" individuelles, mais par une suite alĂ©atoire d'Ă©vĂšnements couvrant une certaine zone gĂ©ographique et ayant un certain impact. AprĂšs avoir prĂ©sentĂ© le modĂšle, nous discuterons de tout le travail qui reste Ă  faire pour mettre au point des mĂ©thodes statistiques permettant de reconstruire certains de ses paramĂštres-clĂ©s Ă  partir de donnĂ©es.

  • Lundi 5 novembre 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Zhanwu Dai
    (INRA, UMR EGFV, Centre Inra de Nouvelle-Aquitaine-Bordeaux)
    SÉMINAIRE ANNULÉ
  • Lundi 22 octobre 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Christian Vestergaard
    (CNRS-Institut Pasteur, Decision and Bayesian Computation Group)
    Randomized reference models for studying contagion processes in temporal networks

    Many dynamical systems can successfully be analyzed using the temporal network formalism. This is notably the case for the human interaction networks that support the spread of contagion and information processes in the population. Empirical temporal networks and dynamic processes that take place in these situations show heterogeneous, non-Markovian, and intrinsically correlated dynamics, making their analysis particularly challenging. Randomized reference models (RRMs) constitute a versatile toolbox for studying such systems. Defined as ensembles of random networks with given features constrained to match those of an input (empirical) network, they may be used to identify statistically significant features in empirical temporal networks (i.e.\ different from the null random networks) and to infer the effects of such features on dynamical processes unfolding in the network. However, the effects of most randomization procedures on temporal network features remain poorly understood, rendering their use non-trivial and susceptible to misinterpretation.

    Here we propose a unified framework for classifying and understanding microcanonical RRMs (MRRMs), which constrain chosen features to take exactly the same value as in the empirical network but are otherwise random. The framework lets us order MRRMs and deduce their effects on important temporal network features, and we use it to show how we may generate new MRRMs from existing ones by sequential composition of independent MRRMs. We show how to apply the framework to unravel how different features of an empirical network of mobile-phone calls influence the spread of information.

  • Lundi 8 octobre 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Elie Desmond-Le Quéméner
    (INRA LBE)
    Microbial thermodynamics: a tool for the understanding of mixed microbial consortia

    The research conducted at LBE (Laboratoire de Biotechnologie de l'Environnement) in Narbonne aims to develop the concept of environmental biorefinery, which consists in treating the by-products of human activities (waste, agricultural residues, effluents) and valorizing them as resources of industrial interest (bioenergies, biomolecules, organic amendment), meanwhile minimizing their environmental and health impact. In this context, my research focus on the understanding of the rules governing microbial ecosystems at the heart of the functioning of these bioprocesses, in order to develop explanatory and predictive models for their control and optimization. Here I will present how microbial thermodynamics can be used:
    -together with microbial transition state theory to predict microbial dynamics and functional community patterns
    -to reveal extracellular electron transfer mediated parasitism
    -to predict individual Operational Taxonomic Units (OTU) metabolism in mixed culture processes
    I will also introduce few perspectives for the modelling of biotic interactions in microbial ecosystems (time-series analyses, individual based modelling...)

  • Lundi 10 septembre 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Liubov Tupikina
    (Ecole Polytechnique, Laboratoire de Physique de la MatiÚre Condensée)
    Heterogeneous network models
    I will present some of our results on the general topic of random walks on networks. In particular, I will discuss the results about how heterogeneities influence the dynamics of, so-called, continuous time random walk model. Then I will talk about some further applications of the heterogeneous random walk model to the epidemiological models on static and temporal networks. Moreover, I will present some results of the discrete-state models on small graphs and discuss the challenges of the presented methods applications to the real data.
     
  • Lundi 11 juin 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Tatiana Giraud
    (CNRS-UPS-AgroParisTech, UMR UMR 8079 Ecologie, Systématique et Evolution)
    Domestication of cheese-making fungi

    Domestication is an excellent model for studies of adaptation because it involves recent and strong selection on a few, identified traits. Few studies have focused on the domestication of fungi, despite their importance to bioindustry and to a general understanding of adaptation in eukaryotes. Penicillium fungi are ubiquitous molds among which two distantly related species have been independently selected for cheese-making, P. roqueforti for blue cheeses like Roquefort, and P. camemberti for soft cheeses like Camembert. The selected traits include morphology, aromatic profile, lipolytic and proteolytic activities, and ability to grow at low temperatures, in a matrix containing bacterial and fungal competitors. By comparing the genomes of ten Penicillium species, we show that adaptation to cheese was associated with multiple recent horizontal transfers of large genomic regions carrying crucial metabolic genes. We identified seven horizontally-transferred regions (HTRs) spanning more than 10 kb each, flanked by specific transposable elements, and displaying nearly 100% identity between distant Penicillium species. Two HTRs carried genes with functions involved in the utilization of cheese nutrients or competition and were found nearly identical in multiple strains and species of cheese-associated Penicillium fungi, indicating recent selective sweeps; they were experimentally associated with faster growth and greater competitiveness on cheese and contained genes highly expressed in the early stage of cheese maturation.

     We also used population genomics to reconstruct the evolutionary history of Penicillium roqueforti. Four populations were identified, including two containing only cheese strains (one corresponding to the emblematic Roquefort “protected designation of origin” strains), and two non-cheese populations including silage and food-spoiling strains. Approximate Bayesian computation analyses indicated that the two cheese populations were derived from independent domestication events. The non-Roquefort population had experienced a stronger genetic bottleneck and displayed greater fitness for traits related to industrial cheese maturation, such as greater lipolysis, cheese cavity colonization and salt tolerance. It probably originated from the industrial selection of a single clonal lineage and is used worldwide for the production of all types of blue cheese other than Roquefort. The Roquefort population resulted from a softer domestication event, probably due to the ancient use of different strains across multiple farms, with possible selection for slower growth before the invention of refrigeration and for greater spore production on the traditional multiplication medium (bread). We detected genomic regions affected by recent positive selection and genomic islands specific to one of the cheese populations, some of which corresponded to putative horizontal gene transfer events. This study sheds light on the processes of adaptation to rapid environmental changes, has industrial implications and raises questions about the conservation of genetic resources.

  • Lundi 28 mai 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Zhanwu Dai
    (INRA, UMR EGFV, Centre Inra de Nouvelle-Aquitaine-Bordeaux )
    Reporté à une date ultérieure
  • Lundi 14 mai 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Elie Desmond-Le Quéméner
    (INRA, UR 050 LBE)
    Séminaire repoussé à une date ultérieure
  • Lundi 23 avril 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Victor Picheny
    (INRA, MIA Toulouse)
    Séminaire reporté à une date ultérieure
  • Lundi 9 avril 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Andreas DrÀger
    (University of Tubingen, Applied Bioinformatics Group)
    TBA
  • Lundi 26 mars 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Julie Josse
    (École Polytechnique, CMAP)
    Distributed Multi-Level Matrix Completion for Medical Databases

    Gathering the information contained in the databases of several hospitals is a step toward personalized medical care as it increases the chances of finding similar patient profiles and therefore provinding them better treatment. However, there are technical (computations and storage issues) and social barriers (privacy concerns) to the aggregation of medical data. Both obstacles can be overcome by turning to distributed computations so that hospitals only share some intermediate results instead of the raw data. As it is often the case, the medical databases are incomplete. One aim of the project is to impute the data of one hospital using the data of the other hospitals. This could also be an incentive to encourage the hospitals to participate in the project. In this talk, we will describe a single imputation method for multi-level (hierarchical) data that can be used to impute both quantitative, categorical and mixed data. This method is based on multi-level simultaneous component analysis (MLSCA) which basically decomposes the variability in both a between and within (hospitals) variability and performs a SVD on each part. The imputation method can be seen as an extension of matrix completion methods. The methods and their distributed versions are implemented in an R package.

  • Lundi 12 mars 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Marie-Laure Martin-Magniette
    (AgroParisTech-INRA, UMR 518 MIA / IPS2)
    ModÚles de mélange gaussiens et de graphes pour mieux comprendre la réponse aux stress de la plante modÚle Arabidopsis
    Les donnĂ©es transcriptomiques sont la source la plus importante de donnĂ©es disponibles actuellement sur les organismes. Pour la plante modĂšle Arabidopsis, son transcriptome a pu ĂȘtre mesurĂ© dĂšs le dĂ©but des annĂ©es 2000 par des puces Ă  ADN. Depuis 2010, nous avons entrepris de faire une Ă©tude globale de la rĂ©ponse aux stress de la plante. AprĂšs une prĂ©sentation du contexte biologique, je prĂ©senterai une analyse de coexpression de plusieurs catĂ©gories de stress rĂ©alisĂ©e Ă  l'aide de mĂ©langes gaussiens. Dans la seconde partie de l'exposĂ©, je vous prĂ©senterai la construction d'un rĂ©seau de corĂ©gulation et comment l'Ă©tude de ce rĂ©seau avec un modĂšle de mĂ©lange de graphes (Stochastic Block Models) a permis de prĂ©dire la fonction de certains gĂšnes mal annotĂ©s et d'identifier une voie de rĂ©ponse aux stress.
  • Lundi 12 mars 2018 - - Salle de rĂ©union 142, bĂątiment 210
    El Houcine Bergou
    (INRA MaIAGE et KAUST-VCC)
    Random Direct Search Method for Unconstrained Smooth Minimization

    In this work we consider the problem of unconstrained minimization of a smooth function in $R^n$ in a setting where only function evaluations are possible. We design a novel random direct search (RDS) method and analyze its complexity. At each iteration, RDS generates a random search direction according to a certain fixed probability law. Our assumptions on this law are very mild. For instance, we allow for the uniform distribution on the sphere and also distributions that concentrate all measure on a positive spanning set. Given a current iterate $x$, the objective function is compared at three points: $x$, $x+\alpha s$ and $x-\alpha s$, where $\alpha>0$ is a stepsize parameter and $s$ is the random search direction. The best of these three points is the next iterate. The complexity of RDS depends on the probability law via a simple characteristic closely related to the cosine measure which is used in the analysis of deterministic direct search (DDS) methods. Unlike in DDS, where $O(n)$ function evaluations must be performed in each iteration in the worst case, our random search method only requires two new function evaluations per iteration. Consequently, while DDS depends quadratically on $n$, our method depends linearly on $n$.

  • Lundi 12 fĂ©vrier 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Séminaire annulé
    TBA
  • Lundi 29 janvier 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Elad Noor
    (ETH Zurich, Institute of Molecular Systems Biology)
    The Hidden costs of enzymatic catalysis
    The existence of a trade-off between the biomass yield and growth rate of cells has been used to explain aerobic fermentation in cancer cells (Warburg effect), yeast cells (Crabtree effect) and in bacteria such as E. coli. This trade-off relies on the assumption that even though fermentation pathways produce 5-10 times less ATP per glucose, respiration requires so much more resources and is therefore inefficient when carbon is not limiting. Is this trade-off a universal constraint imposed by thermodynamics, or a coincidental feature of the specific enzyme kinetic parameters that evolved in these organisms? To answer this question we developed a new method called Enzyme-Flux Cost Minimization (EFCM) to model the costs of both respiration and fermentation (along with ~1000 other flux combinations called elementary flux modes). We find that the trade-off in E. coli is not universal and depends strongly on the availability of oxygen. This framework successfully predicts in vivo enzyme concentrations, and has applications in metabolic engineering where similar candidate pathways can be compared not just by their yields, but also by their costs.
     
  • Lundi 15 janvier 2018 - - Salle de rĂ©union 142, bĂątiment 210
    Nicolas Brunel
    (ENSIIE, UEVE-CNRS UMR 8071 LaMME)
    Optimal control and additive perturbations helps in estimating ill-posed and uncertain dynamical systems

    Ordinary Differential Equations (ODE) are routinely calibrated on real data for estimating unknown parameters or for reverse-engineering of biological systems. Nevertheless, standard statistical technics can give disappointing results because of the complex relationship between parameters and states, that makes the corresponding estimation problem ill-posed. Moreover, ODE are mechanistic models that are prone to modelling errors, whose influences on inference are often neglected during statistical analysis. We propose a regularised estimation framework that consists in adding a perturbation to the original ODE. This perturbation facilitates data fitting and represents also possible model misspecifications, so that parameter estimation is done by solving a trade-off between data fidelity and model fidelity. We show that the underlying optimisation problem is an optimal control problem, that can be solved by the Pontryagin Maximum Principle for general nonlinear and partially observed ODE. The same methodology can be used for the joint estimation of finite and time-varying parameters. We show, in the case of a well-specified parametric model, that our estimator is consistent and reaches the root-$n$ rate. Numerical experiments considering various sources of model misspecifications shows that Tracking still furnish accurate estimates.

    In a second part of the talk, we present a novel algorithm that deals directly with the log-likelihood of the observations and avoid the use of a nonparametric proxy. The inference still uses a perturbed model that is estimated based on the discretisation in time and of the perturbation function (piecewise constant). We focus on linear ODEs and show that we can compute efficiently the parameter estimator by dynamic programming.

    The computational speed enables to address the estimation of (relatively) high-dimensional systems, and to implement standard computationally intensive procedure such as cross-validation and bootstrap. Finally, we show that we can use our approach to estimate nonlinear ODEs used for modelling bacteria interactions in microbiome.

  • Lundi 18 dĂ©cembre 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Jessica Tressou
    (AgroParisTech-INRA, UMR 518 MIA)
    Principes du bootstrap en iid et pour les chaĂźnes de Markov : application aux valeurs extrĂȘmes
    AprĂšs avoir prĂ©sentĂ© briĂšvement le principe du bootstrap dans le cadre indĂ©pendant et identiquement distribuĂ© et montrĂ© sur un exemple ses limites dans le cas des valeurs extrĂȘmes, je montrerai comment on peut dĂ©finir une procĂ©dure de bootstrap rĂ©gĂ©nĂ©ratif sur des chaines de Markov. Cette mĂ©thode, dĂ©veloppĂ©e par Patrice Bertail et StĂ©phane ClĂ©mençon dans un article publiĂ© dans Bernoulli en 2006, a ensuite Ă©tĂ© appliquĂ©e Ă  plusieurs types d'estimateurs (U-statistiques, R-statistiques, L-statistiques).
    Je prĂ©senterai plus en dĂ©tails les applications en lien avec les valeurs extrĂȘmes.
     
    Bertail P, Clémençon S, Tressou J, 2009. Extreme value statistics for Markov chains via the (pseudo-)regenerative method. Extremes, 12(4), 327-360.
    Bertail P, Clémençon S, Tressou J, 2013. Regenerative Block-Bootstrap Confidence Intervals for Tails and Extremal Indexes. Electronic Journal of Statistics, 7, 1224-1248.

2017

  • Lundi 4 dĂ©cembre 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Corinne Vacher
    (Université de Bordeaux - INRA, UMR Biogeco)
    La phyllosphĂšre: une jungle microbienne Ă  l'interface plante-climat

    Les plantes interagissent avec des microorganismes tout au long de leur cycle de vie. Ces communautés microbiennes modulent la résistance des plantes contre de multiples stress, influencent leur capacité à se reproduire et contribuent au fonctionnement des écosystÚmes. Les communautés microbiennes de la phyllosphÚre, qui sont à l'interface entre la plante et l'atmosphÚre, pourraient jouer un rÎle non-négligeable dans la réponse des plantes au changement climatique. Une meilleure prise en compte de ces communautés, voire leur pilotage, pourrait permettre de renforcer la résilience des écosystÚmes. Afin de prédire la dynamique et l'évolution des communautés microbiennes de la phyllosphÚre, il est nécessaire de disposer d'un cadre conceptuel qui intÚgre la plante, l'atmosphÚre et les microorganismes qui se développent à leur interface. En s'appuyant sur la théorie de Mark Vellend (2010), nous avons proposé d'analyser la structure de ces communautés à l'aide de quatre processus: la sélection, la diversification évolutive, la dispersion et la dérive. La sélection, qui est exercée par les traits et le microclimat foliaires et par les interactions entre microorganismes, filtre parmi tous les microorganismes présents sur les feuilles ceux qui parviendront à y survivre. La diversification évolutive génÚre des microorganismes possédant de nouvelles fonctions, tandis que la dispersion permet à certains microorganismes d'atteindre l'habitat foliaire. La dérive écologique, quant à elle, crée des variations aléatoires dans les abondances microbiennes. Si de nombreuses études démontrent l'influence de ces quatre processus, il reste encore de nombreux défis à relever. L'un d'entre eux, que j'illustrerai plus en détail, est de parvenir à déchiffrer les réseaux d'interactions entre les microorganismes à partir de données de séquençage haut-débit et d'évaluer l'influence de ces réseaux sur la santé de la plante.

    Travail joint avec Charlie PAUVERT, Arndt HAMPE, Annabel PORTE, Ursula SAUER, Stéphane COMPANT et Cindy MORRIS et basé sur les publications

    Vacher, C., Hampe, A., Porté, A.J., Sauer, U., Compant, S., Morris, C.E. 2016. The phyllosphere: microbial jungle at the plant-climate interface. Annual Review of Ecology, Evolution, and Systematics, 47: 1-24

    Stéphane A.P. Derocles, David A Bohan, Alex J. Dumbrell, James J.N. Kitson, François Massol, Charlie Pauvert, Manuel Plantegenest, Corinne Vacher, Darren M. Evans. Biomonitoring for the 21st century: integrating Next Generation Sequencing with Ecological Network Analysis. Submitted to Advances in Ecological Research.

  • Lundi 20 novembre 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Marco Bellinzoni
    (Institut Pasteur, Unité de Microbiologie Structurale )
    Unusual features in actinobacterial α-ketoacid dehydrogenase complexes: challenging old beliefs
    Pyruvate dehydrogenase (PDH) and α-ketoglutarate dehydrogenase (KDH) are conserved, tripartite enzyme complexes that carry out essential reactions in central metabolism. Despite their overall conservation, an increasing amount of evidence points to the presence of unique features in the structure and regulation of the actinobacterial versions of these enzymes. I will resume our ongoing studies on the components of a mixed PDH/KDH supercomplex from the species Mycobacterium smegmatis and Corynebacterium glutamicum, and show how our findings challenge old paradigms about the architecture and the regulation of these fascinating enzymatic machines.
  • Lundi 6 novembre 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Sylvain Billiard
    (Université de Lille - CNRS, Unité Evo-Eco-Paléo)
    Rejuvenating functional responses with the renewal theory

    How fast do interactions occur between individuals is central in ecology. Functional responses are classically used to describe the number of predation, mating, competition, etc. in a given timeframe. Hundreds of different forms of functional responses have been proposed in the ecological and mathematical literature. It is well known that this form can dramatically affect the stability and dynamics of populations. Yet, the forms given to functional responses are generally poorly justified from the individual point of view, most ecologists generally adopting a phenomenological approach, in a purely deterministic framework. Here, we propose a novel and original stochastic approach based on the renewal theory. We show how it is possible to derive classical and novel functional responses from the behaviors of the individuals by modelling the time taken by all activities that must be fulfilled for an interaction to be successful. We show how a stochastic approximation of the functional responses can be obtained thanks to the renewal theory. We give applications of our theoretical framework and discuss the importance of interactions as a source of stochasticity in ecological models.

  • Lundi 16 octobre 2017 - - Salle de rĂ©union 142, bĂątiment 210
    séance annulée
    -
  • Lundi 2 octobre 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Pierre Latouche
    (Université Paris 1, Laboratoire SAMM)
    Multiple change points detection and clustering in dynamic networks

    The increasing amount of data stored in form of dynamic interactions between actors necessitates the use of methodologies to automatically extract relevant information. The interactions can be represented by dynamic networks in which most existing methods look for clusters of vertices to summarize the data. In this work, a new framework is proposed in order to cluster the vertices while detecting change points in the intensities of the interactions. These change points are key in the understanding of the temporal interactions. The model used involves non-homogeneous Poisson point processes with cluster dependent piecewise constant intensity functions and common discontinuity points.  A variational expectation maximization algorithm is derived for inference. We show that the pruned exact linear time method, originally developed for change points detection in univariate time series, can be considered for the maximization step. This allows the detection of both the number of change points  and their location.  Experiments on artificial and real datasets are carried out and the proposed approach is compared with related methods.

  • Lundi 18 septembre 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Wolfram Liebermeister
    (INRA, MaIAGE)
    Enzyme and flux cost functions for metabolic modelling

    Metabolic fluxes in cells are sometimes assumed to reflect an economical use of enzymes. To apply this principle of "minimal enzyme cost" and to study its consequences, I consider two types of computational metabolic models: models in which fluxes are given and metabolite levels are chosen to minimise enzyme cost; and models in which the fluxes themselves are optimised. In the first case, optimal enzyme and metabolite levels can be computed by solving a convex optimization problem, and connections between enzyme cost and reaction thermodynamics become very apparent. In the second case, we obtain a nonlinear version of flux balance analysis, a method commonly used in flux prediction. Optimal metabolic fluxes can be found be screening the elementary flux modes, a well-defined set of maximally sparse flux distributions. Altogether, we obtain a tractable method for predicting optimal fluxes, metabolite levels, and enzyme levels in kinetic models of central metabolism.

  • Lundi 26 juin 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Annulé (AG du département MIA)
    TBA
  • Lundi 12 juin 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Christophe Biernacki
    (Université de Lille 1, Modal, INRIA Lille Nord-Europe)
    About two disinherited sides of statistics: data units and computational saving
    Statistics often focuses on designing models, theoretical estimates, related algorithms and model selection. However, some sides of this whole process are somewhat not really tackled by statisticians, leaving the practitioner with some empirically choices, thus poor theoretical warranties. In this context, we identify two situations of interest which are firstly the data unit definition, in case where the practitioner hesitates between few, and secondly the way of saving computational time, for instance by early stopping rules of some estimating algorithms.
    In the first case (data units), we highlight that it is possible to embed data unit selection into a classical model selection principle. We introduce the problem in a regression context before to focus on the model-based clustering and co-clustering context, for data of different kinds (continuous, categorical). It is a joint work with Alexandre Lourme (University of Bordeaux).
    In the second case (computational saving), we recall that an increasingly recurrent statistical question is to design a trade-off between estimate accuracy and computation time. Most estimates practically arise from algorithmic processes aiming at optimizing some standard, but usually only asymptotically relevant, criteria. Thus, the quality of the resulting estimate is a function of both the iteration number and also the involved sample size. We focus on estimating an early stopping time of a gradient descent estimation process aiming at maximizing the likelihood in the simplified context of linear regression (with some discussion in other contexts). It is a joint work with Alain Célisse and Maxime Brunin (University of Lille and Inria, both).
     
  • Lundi 29 mai 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Jean-Christophe Palauqui
    (Institut Jean-Pierre Bourgin, UMR1318 INRA-AgroParisTech)
    Etude descriptive et modélisation du développement embryonnaire chez Arabidopsis thaliana
    L’acquisition des grands plans d’organisation d’une plante s’opĂšre au cours de l’embryogenĂšse prĂ©coce. Ce processus, trĂšs stĂ©rĂ©otypĂ© chez Arabidopsis, conduit Ă  la mise en place d’un organisme pluricellulaire, permettant d’assurer les fonctions essentielles aprĂšs germination de la graine. Au travers d’une Ă©tude descriptive 3D de l’organisation des plans de divisions, qui retrace l’historique des Ă©vĂ©nements de divisions cellulaires, nous dĂ©veloppons une approche de modĂ©lisation de la division cellulaire basĂ©e sur la gĂ©omĂ©trie de la cellule.
     
  • Lundi 24 avril 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Vincent Briane
    (INRIA, IRISA Rennes)
    An adaptive statistical test to detect non Brownian diffusion from particle trajectories

    (This is a joint work with the supervisors of my PhD Myriam Vimond and Charles Kervrann.)

    Assessing the dynamics of particles inside live cell is of paramount interest to understand cell mechanisms. In this presentation, we assume that the motions of particles follow a certain class of random process: the diffusion processes. Our contribution is to propose a statistical method able to classify the motion of the observed trajectories into three groups: subdiffusion (the particle is trapped in a confined domain or moves in a crowded area), superdiffusion (the particle moves in a specific direction thanks to a molecular motor) and free diffusion (namely Brownian motion). This method is an alternative to Mean Square Displacement (MSD) analysis. We assess our procedure on both simulations and real cases.

  • Lundi 15 mai 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Melina Gallopin
    (UPSay, I2BC)
    Nonlinear network-based quantitative trait prediction from transcriptomic data

    Quantitatively predicting phenotype variables by the expression changes in a set of candidate genes is of great interest in molecular biology but it is also a challenging task for several reasons. First, the collected biological observations might be heterogeneous and correspond to different biological mechanisms. Secondly, the gene expression variables used to predict the phenotype are potentially highly correlated since genes interact through unknown regulatory networks. In this talk, we present a novel approach designed to predict quantitative traits from transcriptomic data, taking into account the heterogeneity in biological observations and the hidden gene regulatory networks. The proposed model performs well on prediction but it is also fully model-based, which facilitates the downstream biological interpretation. The model provides clusters of individuals based on the relation between gene expression data and the phenotype, and also leads to infer a gene regulatory network specific for each cluster of individuals.
    We perform numerical simulations to demonstrate that our model is competitive with other prediction models, and we demonstrate the predictive performance and the interpretability of our model to predict olfactory behavior from transcriptomic data on real data from Drosophila Melanogaster Genetic Reference Panel (DGRP).

  • Lundi 20 mars 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Simon Cauchemez
    (Institut Pasteur, Unité Modélisation mathématique des maladies infectieuses)
    Statistical analysis and modelling of epidemics

    Looking at a number of examples from recent outbreaks such as Zika in the Americas, Ebola in West Africa or MERS-CoV in the Middle, I will discuss the challenges associated with the analysis of epidemic data as well as the role modelling can play in the management of these epidemics.

  • Lundi 6 mars 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Séminaire annulé
    -
  • Lundi 20 fĂ©vrier 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Coralie Fritsch
    (Institut Elie Cartan Lorraine, Inria, Ă©quipe Tosca)
    Dynamique adaptative de populations bactériennes dans un bioréacteur

    Je présenterai une approche numérique pour déterminer la possibilité d'invasion de populations bactériennes mutantes dans un bioréacteur. Pour cela, j'introduirai tout d'abord un modÚle stochastique individu-centré et un modÚle déterministe EDP de croissance-fragmentation. Je présenterai le lien entre les fitness d'invasion de ces deux modÚles, définies respectivement par la probabilité de survie de la population et la valeur propre principale d'un certain opérateur. Je donnerai ensuite des résultats sur les variations de ces fitness d'invasion en fonction d'un paramÚtre environnemental du modÚle (agissant sur la croissance et la division des bactéries). Enfin, j'appliquerai ces résultats à un modÚle de bioréacteur afin de décrire une méthode numérique pour l'étude de la dynamique adaptative de la population bactérienne.

  • Lundi 13 fĂ©vrier 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Vyacheslav Kungurtsev
    (Czech Technical University in Prague, Department of Computer Science)
    Nonlinear Optimization, Algorithms for Problems Satisfying Weak Geometric Assumptions, and Problems Arising in Big Data

    In this talk I will present an introduction to the foundations of mathematical programming and optimization theory and review my work on algorithms provably convergent under weak problem assumptions and parallel algorithms suitable for large scale big data machine learning applications. Nonlinear continuous optimization is a mature and active field shown to be effective in solving problems arising from a myriad of applications, including engineering and data science. How wide a class of problems a particular algorithm is capable of solving depends on the algorithm formulation being able to take advantage of problem structure and geometric properties. Algorithms that converge reliably and quickly for a broad range of constrained medium-scale problems have been developed recently, advancing the state of the art. The age of big data necessitates the use of parallel architectures in the computation of an algorithm''s procedural steps. I will present a framework for using problem structure to quickly and reliably solve a large scale nonconvex optimization problem as would arise in machine learning.

  • Lundi 30 janvier 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Kevin Cohen
    (University of Colorado, School of Medicine, Biomedical Text Mining Group, Computational Bioscience Program, USA)
    Synthetic lethal screen reranking with natural language processing

    High-throughput assays are an experimental paradigm that has revolutionized biology in recent decades.  They are very powerful, but produce large numbers of false positives that impede the efficient analysis of experimental data.  This talk will discuss the use of natural language processing to deal with the false positive problem, focusing on a type of high-throughput assay called a synthetic lethal screen.

  • Lundi 16 Janvier 2017 - - Salle de rĂ©union 142, bĂątiment 210
    Etienne Birmelé
    (Paris Descartes, MAP5)
    Étude des perturbations des rĂ©gulations gĂ©niques dans le cancer de la vessie

    Certains types de cancer de la vessie se comportent de façon trÚs similaires à des cellules saines en train de proliférer ou de se différencier, laissant à penser que la tumeur utilise principalement le réseau de régulation préexistant. Le projet LIONS, en collaboration avec les universités d'Evry, Montpellier et York, ainsi qu'avec l'institut Curie, cherche à déterminer quels sont les facteurs de transcription (TF) clés de la tumorigénÚse, c'est-à-dire dont une altération dans le mécanisme de régulation rend le comportement cellulaire pathologique.

    D'un point de vue statistique, cette notion, différente de celle d'expression différentielle, nécessite de 1) disposer d'un réseau de référence correspondant aux cellules saines 2) établir un score de dérégulation indiquant quels gÚnes ont un comportement anormal dans la tumeur 3) déterminer les jeux minimaux de TF expliquant ces comportement anormaux. L'exposé développera les pistes étudiées actuellement pour chacun de ces points.

  • Lundi 21 novembre 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Olivier Martin
    (AgroParisTech-INRA, UMR 791 MOSAR)
    Modéliser la dynamique des phénotypes à l'échelle de la vie de l'animal : du systÚme biologique au systÚme d'élevage

    Le couplage de modÚles dynamiques est illustré à travers un exemple appliqué aux performances productives et reproductives des vaches.Ce travail s'inscrit dans le cadre du projet européen PROLIFIC et aborde la problématique générale de la fertilité des troupeaux bovins laitiers.

    Un modĂšle des performances Ă  l'Ă©chelle de la vie de la vache est d'abord briĂšvement prĂ©sentĂ© (GARUNS : Martin and Sauvant, 2010ab). Ce modĂšle est lui-mĂȘme le produit d'un couplage d'un modĂšle thĂ©orique de la dynamique des prioritĂ©s entre fonctions vitales et d'un modĂšle de partition de l'Ă©nergie. Dans le modĂšle GARUNS, il n'y a pas de rĂ©gulation de la reproduction et les dates d'insemination associĂ©es Ă  chaque cycle de reproduction sont contrĂŽlĂ©es par des paramĂštres fixes. Ce modĂšle du systĂšme animal n'intĂšgre donc pas de variabilitĂ© des performances de reproduction qui dĂ©coule en particulier des anomalies de cyclicitĂ© et des Ă©checs d'insĂ©mination.

    Un modÚle du systÚme reproducteur est ensuite présenté (RPM : Martin et al., 2012). Ce modÚle repose sur un cadre conceptuel générique de représentation du fonctionnement du systÚme reproducteur. Dans ce modÚle, des unités biologiques, comme l'ovaire, l'uterus ou l'embryon, intéragissent par des signaux hormonaux qui régulent leur dynamique de fonctionnement. Selon leur état de compétence, ces unités sont susceptibles de produire un signal ou de réagir à un signal en changeant d'état de compétence. Chaque unité biologique est ainsi décrite à travers des cycles d'états de compétence dont la dynamique est régulée par la signalisation hormonale. La dynamique d'ensemble du systÚme reproducteur émerge de la dynamique des intéractions entre les unités biologiques. Dans le cadre du projet PROLIFIC, deux modÚles ont été développés sur la base de ce cadre conceptuel: une version dite "heavy", conçue pour intégrer les connaissances sur les mécanismes physiologiques sous-jacents et une version dite "lite", conçue pour le couplage avec le modÚle GARUNS, permettant de simuler la variabilité des cycles de reproduction et opérationnelle pour des simulations rapides.

    Un modÚle du systÚme d'élevage est briÚvement décrit (BSM : Friggens et al., unpub.). A l'échelle du troupeau, ce modÚle contrÎle les individus vaches représentées explicitement par des versions du couple GARUNS-RPMlite. Ce modÚle simule les pratiques d'élevage en matiÚre d'alimentation, de détection d'oestrus, d'insémination, de réforme et de sélection génétique. En pratique, pour le couplage entre GARUNS et RPMlite à l'échelle d'un individu, BSM fournit le temps d'insémination au modÚle GARUNS, décidé sur la base de la dynamique d'estrus produite par RPMlite et de la stratégie de reproduction fixée.

    Le couplage des trois modÚles est ensuite abordé. Cette exemple correspond au couplage de modÚles basés sur des concepts similaires (priorié vs états de compétence), avec des pas de temps différents (d vs h) et dans un contexte du couplage d'un modÚle avec un modÚle pré-existant. Le modÚle GARUNS produit une dynamique de performances productives, en particulier le bilan énergétique et le niveau de turnover du métabolisme énergétique, qui régule RPMlite. RPMlite produit une signalisation hormonale (potentiellement lue et interprétée par BSM) et fournit les temps de conception au modÚle GARUNS. Le modÚle GARUNS enclenche alors un cycle de reproduction enchaßnant une gestation et une lactation (ou éventuellement un avortement qui interromp le cycle), ce qui modifie la dynamique du métabolisme énergétique. Le couplage repose ainsi sur une boucle entre les modÚles GARUNS et RPMlite.

    Des résultats de simulations sont enfin présentés pour illustrer la variabilité des performances de reproduction produites par le couplage GARUNS+RPMlite.

2016

  • Lundi 5 dĂ©cembre 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Michael Blum
    (Univ. J. Fourier & CNRS, Laboratoire TIMC-IMAG, Grenoble)
    Détection d'outliers en grande dimension: application à la génomique des populations

    Notre objectif est de dĂ©tecter quelles sont les variables outliers dans des jeux de donnĂ©es de grande dimension. Les mĂ©thodes de dĂ©tection d'outliers sont utilisĂ©es en gĂ©nomique pour dĂ©tecter quels sont les gĂšnes qui permettent aux individus de s’adapter Ă  leur environnement. Nous proposons une approche rapide basĂ©e sur l’analyse en composantes principales. Le principe est de considĂ©rer comme gĂšnes candidats ceux qui sont excessivement corrĂ©lĂ©s avec les composantes principales. Pour ce faire, nous calculons pour chaque marqueur gĂ©nĂ©tique un vecteur qui mesure l’association entre un marqueur gĂ©nĂ©tique et les composantes principales. Nous utilisons ensuite la distance de Mahalanobis pour trouver quels sont les vecteurs atypiques. En utilisant un jeu de donnĂ©es humains comprenant un peu plus d’un millier d’individus et des centaines de milliers de marqueurs gĂ©nĂ©tiques, nous montrons que cette approche permet de dĂ©tecter des exemples d’adaptation biologique chez l’homme.

  • Lundi 7 novembre 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Kevin Cohen
    (University of Colorado, School of Medicine, Biomedical Text Mining Group)
    Annulé
  • Lundi 17 octobre 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Florence DĂ©barre
    (CollĂšge de France, CIRB CNRS UMR 7241)
    Evolution in spatially heterogeneous environments

    Understanding the links between the diversity of habitats and biodiversity is a core topic in Ecology, Evolution and Conservation Biology; mathematical and computational models can help provide qualitative answers to this question. During my talk, I will present some ecological (short-term) and evolutionary (long-term) consequences of spatial structuring and environmental heterogeneities. I will start with the example of resistance to treatments, and will show how spatial heterogeneities can help limit the spread of resistance. Then I will move on to a more ecological model and will present results on the evolution of specialist and generalist strategies in a spatially heterogeneous environment. 

  • Lundi 26 septembre 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Tabea Rebafka
    (UPMC, LPMA)
    Estimation et clustering dans un modĂšle de processus de Poisson semiparamĂ©trique Ă  blocs stochastiques pour des rĂ©seaux d’interaction longitudinaux

    In this work, we introduce a Poisson process stochastic block model for recurrent interaction events, where each individual belongs to a latent group and interactions between two individuals follow a conditional inhomogeneous Poisson process whose intensity is driven by the individuals’ latent groups. The model is semiparametric as the intensities per group pair are modeled in a nonparametric way. First an identifiability result on the weights of the latent groups and the nonparametric intensities is established. Then we propose an estimation procedure, relying on a semi parametric version of a variational expectation-maximization algorithm. Two different versions of the method are proposed, using either histogram-type (with an adaptive choice of the partition size) or kernel intensity estimators. We also propose an integrated classification likelihood criterion to select the number of latent groups. Asymptotic consistency results are then explored, both for the estimators of the cumulative intensities per group pair and for the kernel procedures that estimate the intensities per group pair. Finally, we carry out synthetic experiments and analyse several real datasets to illustrate the strengths and weaknesses of our approach.
    This is joint work with Catherine Matias et Fanny Villers.

  • Lundi 12 Septembre 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Ludovic Cottret
    (INRA-CNRS, LIPM)
    Analyse de la robustesse phénotypique d'une bactérie phytopathogÚne par intégration du réseau métabolique et du réseau de régulation

    Dans l'analyse du rĂ©seau mĂ©tabolique, la robustesse d'un phĂ©notype est communĂ©ment dĂ©finie comme la capacitĂ© du mĂ©tabolisme Ă  maintenir ce phĂ©notype malgrĂ© des perturbations gĂ©nĂ©tiques ou environnementales. Plusieurs Ă©lĂ©ments peuvent ĂȘtre Ă  l'origine de la robustesse phĂ©notypique : la versatilitĂ©, i.e. la capacitĂ© du systĂšme Ă  fonctionner Ă  partir de diffĂ©rents nutriments; la redondance fonctionnelle comprenant la redondance gĂ©nĂ©tique et les voies mĂ©taboliques alternatives; et enfin le contrĂŽle du systĂšme qui intervient pour capter et compenser efficacement les perturbations.

    Nous nous sommes intéressés à l'analyse de la robustesse phénotypique chez une bactérie phytopathogÚne, Ralstonia solanacearum. Plus particuliÚrement, nous avons tenté de prédire l'influence du réseau de régulation de la virulence sur la robustesse phénotypique. Pour cela, nous avons d'abord effectué une reconstruction de haute qualité du réseau métabolique et du réseau de régulation de la virulence grùce à une suite d'outils automatiques et semi automatiques. Ensuite, nous avons développé une librairie Java, appelée FlexFlux, destinée à l'analyse de balance des flux (FBA). L'originalité de FlexFlux est d'intégrer de façon native le réseau de régulation et le réseau de métabolique dans chacune de ces fonctions. Nous verrons enfin comment nous avons utilisé FlexFlux pour mesurer la robustesse de plusieurs phénotypes (liés ou non à la virulence) et l'influence du réseau de régulation de la virulence sur celle-ci.

  • Lundi 12 septembre 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Catherine Larédo
    (INRA, MaIAGE)
    Estimation paramétrique pour des équations différentielles stochastiques à effets mixtes à partir de données longitudinales discrétisées.

    Les Ă©quations diffĂ©rentielles stochastiques fournissent un cadre naturel pour modĂ©liser la variabilitĂ© intrinsĂšque inhĂ©rente Ă  de nombreux processus physiques Ă  temps continu.Quand ces processus sont observĂ©s sur plusieurs individus ou unitĂ©s expĂ©rimentales, les Ă©quations diffĂ©rentielles stochastiques Ă  effets mixtes (SDEME) permettent de quantifier simultanĂ©ment la variabilitĂ© intrinsĂšque (variabilitĂ© intra) et la variabilitĂ© entre individus (variabilitĂ© inter). Ces dynamiques modĂ©lisĂ©es par des processus Ă  temps continu sont gĂ©nĂ©ralement observĂ©es avec un certain pas de temps (donnĂ©es discrĂ©tisĂ©es). Du fait de la difficultĂ© Ă  Ă©tudier la vraisemblance, faire l’infĂ©rence Ă  partir d’observations discrĂ©tisĂ©es de SDEME est un problĂšme ouvert d’un point de vue thĂ©orique pour des modĂšles gĂ©nĂ©raux de SDEME. Nous Ă©tudions ici des cas pour lesquels on dispose d’approximations explicites de la vraisemblance.

    Travail en collaboration avec Maud Delattre (AgroParisTech, France) et Valentine Genon-Catalot
    (UMR CNRS 8145, Laboratoire MAP5, Université Paris Descartes, Sorbonne Paris Cité, France).

  • Lundi 20 juin 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Julien Chiquet
    (AgroParisTech-INRA, UMR 518 MIA)
    Fast tree inference with weighted fusion penalties

    Given a data set with many features observed in a large number of conditions, it is desirable to fuse and aggregate conditions which are similar to ease the interpretation and extract the main characteristics of the data. This paper presents a multidimensional fusion penalty framework to address this question when the number of conditions is large. If the fusion penalty is encoded by an ℓq-norm, we prove for uniform weights that the path of solutions is a tree which is suitable for interpretability. For the ℓ1 and ℓ∞-norms, the path is piecewise linear and we derive a homotopy algorithm to recover exactly the whole tree structure. For weighted ℓ1-fusion penalties, we demonstrate that distance-decreasing weights lead to balanced tree structures. For a subclass of these weights that we call “exponentially adaptive”, we derive an O(n log(n)) homotopy algorithm and we prove an asymptotic oracle property. This guarantees that we recover the underlying structure of the data efficiently both from a statistical and a computational point of view. We provide a fast implementation of the homotopy algorithm for the single feature case, as well as an efficient embedded cross-validation procedure that takes advantage of the tree structure of the path of solutions. Our proposal outperforms its competing procedures on simulations both in terms of timings and prediction accuracy. As an example we consider phenotypic data: given one or several traits, we reconstruct a balanced tree structure and assess its agreement with the known taxonomy.

  • Lundi 23 mai 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Annulé
    TBA
  • Lundi 18 avril 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Reporté
    TBA
  • Lundi 4 avril 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Argyris Kalogeratos
    (ENS Cachan, CMLA)
    Algorithmes efficaces pour contenir des processus épidémiques sur réseaux à l'aide de ressources d'efficacité limitée / Suppressing epidemics on arbitrary networks using treatment resources of limited efficiency

    Résumé : Dans de nombreuses situations réelles, il est essentiel de supprimer un processus de diffusion indésirable (virus, information, comportements, etc.) en temps réel. Cette exposé proposera des méthodes pour l'allocation dynamique de ressources pour des épidémies de type SIS (en temps continu) à l'aide d'un budget de ressources limité à disposition des autorités.
    Dans cet exposé, nous montrerons que la structure macroscopique et microscopique du réseau joue un rÎle clé dans l'explication de ces phénomÚnes de propagation et deux stratégies seront présentées : a) une approche simple et focalisée sur l'évolution court-terme du processus, et b) une approche plus sophistiquée qui utilise un ordre de priorité (précalculée avant l'épidémie) spécifiant le déroulement de la stratégie de guérison.

    Abstract: In many real-life situations, it is critical to dynamically suppress or remove an undesired diffusion process (viruses, information, behaviors, etc.). The talk will present a framework for Dynamic Resource Allocation (DRA) assuming a continuous-time SIS epidemic model, and that a budget of treatment resources of limited efficiency are at the disposal of authorities.
    Special emphasis will be given on the macroscopic and microscopic (or local) properties of the network structure for the problem and two strategies will be presented that fall in this framework: a) a simple yet effective greedy approach, and b) a more sophisticated one that uses a precomputed priority plan of how the healing strategy should proceed on a specific network.

  • Lundi 21 mars 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Avner Bar-Hen
    (Université Paris Descartes - CNRS, Laboratoire MAP5, UMR 8145)
    DĂ©tection de cluster spatiale en utilisant la distance au plus proche voisin

    Afin de caractĂ©riser l'impact des processus Ă©cologiques sur la distribution spatiale des espĂšces d'arbres, nous prĂ©sentons une mĂ©thode pour dĂ©tecter des clusters de points. Notre mĂ©thode est basĂ©e sur une transformation itĂ©rative de la distance entre les points. Notre approche a l'avantage d'ĂȘtre indĂ©pendante d'une forme arbitraire de cluster et permet un ajustement pour les covariables. La comparaison de la distance observĂ©e entre points avec un processus de rĂ©fĂ©rence conduit Ă  une classification hiĂ©rarchique des clusters. Le choix du nombre optimal de clusters est effectuĂ©e en utilisant la statistique de Gap. Notre procĂ©dure est illustrĂ©e sur une rĂ©partition spatiale des espĂšces de la Dicorynia en Guyane française.

  • Lundi 7 mars 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Marie Suez
    (UPMC, IBPS / INRA, MaIAGE)
    DiversitĂ© gĂ©nĂ©tique des populations de cerfs Ă©laphe (Cervus elaphus) en Île-de-France, en liaison avec l'anthropisation.

    Au cours des 60 derniĂšres annĂ©es le dĂ©veloppement des infrastructures de transports (Autoroutes, Lignes Grandes Vitesse, Nationales doubles voies) a fragmentĂ© l’habitat des cerfs Ă©laphe (Cervus elaphus). D’aprĂšs les observations naturalistes, cette anthropisation a causĂ© la fragmentation de deux populations gĂ©ographiques existantes en sept dans la partie Sud et d’une en trois dans la partie Nord.

    Afin d’évaluer l’impact de ces infrastructures sur la structuration gĂ©nĂ©tique de ces populations de cerfs, nous avons Ă©chantillonnĂ© chacune de ces populations grĂące Ă  la coopĂ©ration de trois fĂ©dĂ©rations de chasse. Le cours laps de temps Ă©coulĂ© depuis la construction de ces infrastructures nous a conduits Ă  choisir comme marqueurs molĂ©culaires les microsatellites, efficaces dans l’infĂ©rence d’évĂšnements rĂ©cents. Les nouvelles techniques de sĂ©quençages (NGS) permettent d’obtenir d’importants jeux de donnĂ©es rapidement, nous avons choisi d’utiliser ces mĂ©thodes de sĂ©quençage pour obtenir nos donnĂ©es. Aucun logiciel ne permettant de traiter les donnĂ©es de sĂ©quençage haut dĂ©bit des microsatellites pour des espĂšces dont le gĂ©nome n’est pas complĂštement sĂ©quencĂ©, nous avons alors rĂ©alisĂ© un programme, MicNeSs qui permet de gĂ©notyper rapidement et objectivement (sans intervention humaine) un grand nombre d’individus et de locus. Nous avons utilisĂ© MicNeSs pour gĂ©notyper 345 individus pour 17 locus microsatellites. A partir de ce jeu de donnĂ©es, nous avons montrĂ© l’existence d’une structuration gĂ©nĂ©tique des populations de cerfs Ă©laphe en Île-de-France en liaison avec les infrastructures routiĂšres et ferroviaires. Nous avons mis en Ă©vidence un effet fort des jumelages autoroutes/LGV et une efficacitĂ© diffĂ©rentielle des passages grande faune de 2Ăšme et 3Ăšme gĂ©nĂ©ration sur les populations de cerfs Ă©laphe en Île-de-France.

  • Lundi 22 fĂ©vrier 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Youssef Diouane
    (ISAE - SupAĂ©ro)
    Globally convergent evolution strategies with application to an Earth imaging problem in geophysics.

    In recent years, there has been significant and growing interest in Derivative-Free Optimization (DFO). This field can be divided into two categories: deterministic and stochastic. Despite addressing the same problem domain, only few interactions between the two DFO categories were established in the existing literature. In this thesis, we attempt to bridge this gap by showing how ideas from deterministic DFO can improve the efficiency and the rigorousness of one of the most successful class of stochastic algorithms, known as Evolution Strategies (ES’s). We propose to equip a class of ES’s with known techniques from deterministic DFO. The modified ES’s achieve rigorously a form of global convergence under reasonable assumptions. By global convergence, we mean convergence to first-order stationary points independently of the starting point. The modified ES’s are extended to handle general constrained optimization problems. Furthermore, we show how to significantly improve the numerical performance of ES’s by incorporating a search step at the beginning of each iteration. In this step, we build a quadratic model using the points where the objective function has been previously evaluated. Motivated by the recent growth of high performance computing resources and the parallel nature of ES’s, an application of our modified ES’s to Earth imaging geophysics problem is proposed. The obtained results provide a great improvement to known solutions of this problem.

  • Jeudi 18 fĂ©vrier 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Stefanie Widder
    (Division of Computational Systems Biology, Univ. of Vienna)
    (Self)-Organization of the human microbiota in health and disease

    Central to understanding the behavior of microbial communities (MC) are microbial interactions and their organization. In my talk I will give two examples where community-wide interactions lead to complex behavior of the MC. Pathogenic microbiota cause chronic infections in the airways of cystic fibrosis patients. We show that organization into two competing communities and shifts in their interactions, as well as in their metabolic core processes are associated to a shift in the severity of the disease. Using network analysis and the keystone concept, we propose functional and taxonomic keystones as targets for novel drug development. The human gut microbiome is characterized by seemingly opposing trends: stability of a (functional) core community and pronounced variability throughout lifetime and between subjects. We explore self-organization as potential mechanism underlying these evidences. We use a modeling approach to analyze the time behavior of gut microbiota and find that the community stratifies into three sub-groups linked to lifetime and abundance. In time the gut community exhibits pink noise and scale-invariance, hallmarks for self-organization.

  • Lundi 8 fĂ©vrier 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Nicolas BOUSQUET
    (EDF, R&D)
    Estimer des indicateurs de risque par simulation de modÚles complexes "boßte noire" en tirant parti de contraintes de forme (monotonie, convexité...), avec applications en aide à la décision

    Les modĂšles boĂźte noire sont de plus nombreux dans les Ă©tudes d'impact et de gestion des ressources (industrielles, environnementales, etc.). Souvent dĂ©terministes, ils sont utilisĂ©s pour mener des simulations Ă  partir d'entrĂ©es rendues stochastiques, car considĂ©rĂ©es incertaines par essence ou mal connues. Des indicateurs classiques d'aide Ă  la dĂ©cision en sortie des modĂšles sont des probabilitĂ©s de dĂ©passement, ou des quantiles. Ne pouvant entrer dans le modĂšle par des mĂ©thodes intrusives, de nombreuses techniques d'estimation par rĂ©duction de variance (par rapport Ă  des approches de Monte Carlo classiques) de ces indicateurs sont maintenant disponibles. Cependant, elles nĂ©cessitent souvent des hypothĂšses de rĂ©gularité qui ne peuvent ĂȘtre vĂ©rifiĂ©es (et sont fausses lorsque le phĂ©nomĂšne modĂ©lisĂ© subit des effets falaises de perte de continuitĂ©), et ont un coĂ»t de simulation potentiellement trĂšs Ă©levĂ©, qui interdit de les utiliser dans de nombreux cas concrets. C'est pourquoi des mĂ©thodes ont Ă©tĂ© trĂšs rĂ©cemment dĂ©veloppĂ©es pour pallier ces difficultĂ©s en tirant parti des contraintes (ou propriĂ©tĂ©s) de forme s'exerçant sur les sorties de tels modĂšles. Un cas important est celui de la monotonie. La convexitĂ© est aussi Ă©voquĂ©e. L'exposĂ© prĂ©sentera la construction de telles mĂ©thodes, dĂ©taillera les outils thĂ©oriques nĂ©cessaires pour Ă©tudier le bien-fondĂ© de l'estimation, et montrera des rĂ©sultats thĂ©oriques et appliquĂ©s sur des cas industriels.

     

  • Lundi 25 janvier 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Sophie DONNET
    (AgroParisTech-INRA, UMR 518 MIA)
    Bayesian estimation for multidimensional Hawkes processes

    Multidimensional Hawkes processes are used to modelise multivariate neuron spike data. The estimation of intensity functions allows to understand the neuronal interaction structure. In a non-parametric frequentist framework, LASSO estimators have been proposed in the literature. In this work, we propose a Bayesian non-parametric estimation. We sample the posterior distribution through a Sequential Monte Carlo algorithm, well adapted to point processes.

  • Lundi 11 janvier 2016 - - Salle de rĂ©union 142, bĂątiment 210
    Jean-Benoist LEGER
    (INRA, MaIAGE)
    ModĂšle de graphes Ă  espace latent continu de type SBM

    De nombreuses donnĂ©es entre des Ă©lĂ©ments peuvent ĂȘtre prĂ©sentĂ©s sous
    forme de rĂ©seaux. Ces donnĂ©es peuvent ĂȘtre binaires, comme une
    prĂ©sence/absence de relation, quantifiĂ©es, continues, oĂč ĂȘtre valuĂ©es
    sur d'autres espaces. Pour construire un modÚle probabiliste adapté à
    ces données des hypothÚses sont nécessaires. Des hypothÚses
    d'appartenance des nƓuds Ă  des classes latentes et une indĂ©pendance de
    la loi sur les lien conditionnellement à la loi sur les nƓuds conduit à
    des modĂšles de type SBM (Stochastic Block Models). Il est possible de
    relacher la contrainte sur l'espace latent Ă  valeurs discrĂštes pour se
    placer dans un espace latent Ă  valeurs continues, ce qui conduit Ă  des
    modĂšles de type MMSBM. (Mixed Membership SBM).

    Il est également possible de disposer d'information extérieures pouvant
    ĂȘtre introduites sous forme de covariables.

    Cette présentation introduira les modÚles de graphes à classes latentes
    de type SBM, avec ou sans covariables pour diverses lois de
    probabilités sur les liens. Elle présentera l'extension à la classe
    latente continue, et introduira une méthode d'estimation basée sur le
    Variational-EM.

2015

  • Lundi 15 dĂ©cembre 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Edward IONIDES
    (University of Michigan, Department of Statistics )
    Inference for dynamic and latent variable models via iterated, perturbed Bayes maps

    Iterated filtering algorithms are stochastic optimization procedures for latent variable models that recursively combine parameter perturbations with latent variable reconstruction. Previously, theoretical support for these algorithms has been based on the use of conditional moments of perturbed parameters to approximate derivatives of the log likelihood function. We introduce a new theoretical approach based on the convergence of an iterated Bayes map. A new algorithm supported by this theory displays substantial numerical improvement on the computational challenge of inferring parameters of a partially observed Markov process.

  • Lundi 30 novembre 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Pierre Larmande
    (IRD, South Green, Montpellier)
    Enabling knowledge management in the Agronomic Domain

    The drastic growth in data in the recent years, within the Agronomic sciences has brought the concept of knowledge management to the forefront. Some of the factors that contribute to this change include a) conducting high-throughput experiments have become affordable, the time spent in generating data through these experiments are minuscule when compared to its integration and analysis; b) publishing data over the web is fairly trivial and c) multiple databases exist for each type of data (i.e. ‘omics’ data) with a possible overlap or slight variation in its coverage [1, 2]. In most cases these sources remain autonomous and disconnected. Hence, efficiently managed data and the underlying knowledge in principle will make data analysis straightforward aiding in more efficient decision making. We are involved in developing methods to aid data integration and knowledge management within the domain of Agronomic sciences to improve information accessibility and interoperability. To this end, we address the challenge by pursuing several complementary research directions towards: distributed, heterogeneous data integration.

    References:
    Goble, C. and Stevens, R. (2008) State of the nation in data integration for bioinformatics. Journal of Biomedical Informatics, 41(5), 687-693.
    Antezana, E., et al. (2009) Biological knowledge management: the emerging role of the Semantic Web technologies.Brief. in Bioinformatics,10(4), 392-407.

  • Lundi 16 novembre 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Fabrice Rossi
    (équipe SAMM, Université Paris 1 Panthéon-Sorbonne)
    Analyse exploratoire de graphes dynamiques

    Nous étudions dans ce travail des données d'interaction, constituées de triplets source, destination, instant. Ce type de données est assez fréquent dans les relations intermediées informatiquement, comme par exemple les emails (expéditeur, récepteur, horodatage de la connexion au serveur STMP), les appels téléphoniques, les SMS, etc. On peut les voir comme un graphe dynamique : les sources et destinations forment les sommets du graphe, alors que les arcs sont les ntraces des interactions. Une fonction de présence indique si un arc est actif à un instant donné.

    Nous proposons une méthode d'analyse exploratoire de ce type de données par tri-classification : nous construisons des classes de sources, des classes de destinations et des intervalles de temps qui garantissent une forme de stationnarité locale des interactions à l'intersection de trois classes. La méthode proposée ne demande aucun paramÚtre utilisateur et donne des résultats trÚs satisfaisants sur des données réelles volumineuses.

  • Lundi 5 octobre 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Christopher Quince
    (Warwick Medical School, University of Warwick)
    Probabilistic Modelling of Microbial Community Structure

    I will give an overview of methods that use generative probabilistic models to describe microbial community structure as determined through next generation sequencing. I will discuss the concept that the human gut microbiota is derived from a finite number of discrete types or enterotypes. I will show that Dirichlet-multinomial mixtures allow a more nuanced description of enterotypes as diffuse peaks in community configurations. I will then extend the idea of Dirichlet priors for community configurations to hierarchical Dirichlet processes which allow ‘unseen’ species to be modelled. I will use these both as a means for fitting the ecological model, Hubbell’s Unified Neutral Theory of Biodiversity, and for avoiding rarefaction in microbial community diversity estimation.

  • Lundi 2 novembre 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Isabelle Bloch
    (UMR 5141 LTCI, Telecom ParisTech - CNRS)
    ModÚles symboliques pour la reconnaissance de structures dans les images et l'interprétation de scÚnes
  • Vendredi 2 octobre 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Philipp W. Messer
    (Department of Biological Statistics and Computational Biology, Cornell University)
    Understanding the rapid evolution of pesticide and drug resistance

    We typically think of evolution as a slow and gradual process, driven by the accumulation of small changes over millions of years. While it is well known that evolution can be much faster when humans impose artificial selection, for instance during animal breeding, such rapid evolutionary responses are generally thought to rely on the availability of standing genetic variation. However, recent studies of the evolution of pesticide and drug resistance revealed that adaptation can be rapid despite requiring complex alleles that are not initially present in a population. Furthermore, it appears that rapid adaptation does not always follow the classic selective sweep model, but often produces so-called soft selective sweeps, where multiple adaptive alleles of independent mutational origin sweep through the population at the same time. In my talk, I will show how the emerging field of population genomics can help us uncover the mechanisms that underlie these rapid evolutionary responses and explain the frequent occurrence of soft selective sweeps. I will also discuss the challenges this poses for computational approaches aimed at identifying adaptive loci, as well as for our theoretical understanding of adaptive dynamics, and present new strategies for tackling these problems.

  • Lundi 21 septembre 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Jean-Michel Marin
    (UMR CNRS 5149, Institut de Mathématiques et Modélisation, Université de Montpellier)
    MĂ©thodes d’infĂ©rence de l’histoire dĂ©mographique de populations structurĂ©es Ă  partir de donnĂ©es de polymorphisme gĂ©nĂ©tique

    Un des principaux dĂ©veloppements de la modĂ©lisation en gĂ©nĂ©tique des populations est l’utilisation des mĂ©thodes dites coalescentes ou gĂ©nĂ©alogiques. Le but est de reconstruire des Ă©lĂ©ments de l'histoire de populations. Pour examiner la structure des donnĂ©es gĂ©nĂ©tiques, ces mĂ©thodes utilisent l'arbre gĂ©nĂ©alogique des gĂšnes. La formulation d’un modĂšle est contrainte par un scĂ©nario Ă©volutif qui imite la rĂ©alitĂ© historique et dĂ©mographique de l'espĂšce. Un tel scĂ©nario rĂ©sume l’histoire Ă©volutive des populations par une suite d'Ă©vĂ©nements dĂ©mographiques depuis une population ancestrale. Ces Ă©vĂ©nements sont constituĂ©s de divergences, des migrations et des variations de tailles entre les populations. Les jeux de donnĂ©es que l’on considĂšre sont constituĂ©s d'informations gĂ©nĂ©tiques issues de plusieurs locus. Les modĂšles que nous Ă©tudions sont sous l'hypothĂšse de neutralitĂ© qui implique l'absence d'effet de sĂ©lection. Avec ces modĂšles, nous pouvons infĂ©rer de quelle sources ancestrale provient une population rĂ©cente, dĂ©crire des voies d’invasion de populations... Il faut alors utiliser une procĂ©dure de choix de modĂšle, chaque hypothĂšse correspond Ă  un scĂ©nario dĂ©mographique. La plupart du temps, on ne sait pas calculer la vraisemblance de donnĂ©es de polymorphisme. Dans cet exposĂ©, nous prĂ©senterons les dĂ©fis statistiques vĂ©hiculĂ©s par ces modĂšles sans vraisemblance explicite. Puis, nous montrerons comment certaines mĂ©thodes bayĂ©siennes approchĂ©es permettent d’y rĂ©pondre.

  • Lundi 7 septembre 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Wolfram Liebermeister
    (Institut fĂŒr Biochemie, CharitĂ© - UniversitĂ€tsmedizin Berlin)
    Enzyme economy in metabolic networks

    It often assumed - sometimes tacitly - that cells use their enzyme resources economically. This hypothesis can be studied by analysing enzyme requirements in kinetic models. I discuss models in which enzyme levels must realize a given flux distribution. The search for cost-optimal enzyme and metabolite levels can be formulated as a convex optimization problem. The enzyme profiles are shaped by opposing tendencies: a need for sufficient saturation with substrate, and an avoidance of small thermodynamic forces. The resulting prediction of metabolite and enzyme levels can complement constraint-based approaches for flux prediction, such as Resource Balance Analysis.

  • Lundi 29 juin 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Bertrand Cloez
    (INRA-SupAgro, UMR 729 MISTEA, Montpellier)
    Comportement en temps long de processus avec extinction

    Dans cet exposé, nous nous intéresserons à un processus de Markov possédant un état absorbant comme par exemple le nombre d'individus ou la proportion de gÚnes dans une population. L'équilibre d'un tel processus est l'état absorbant (0 pour le nombre d'individu, 0 ou 100% pour la proportion de gÚnes). Il arrive, sous certaines conditions, que ce processus atteigne une sorte d'équilibre avant l'extinction. On parle d'équilibre quasi-stationnaire. Nous décrirons quelques propriétés de celle-ci et donnerons deux algorithmes pour la déterminer. Le premier est une méthode particulaire proche des algorithmes génétiques utilisés en filtrage non-linéaire. Le second est une marche aléatoire renforcée dont l'étude est basée sur les algorithmes stochastiques et la méthode de l'EDO.

  • Lundi 15 juin 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Evelyne Lutton
    (INRA - GMPA)
    Modéliser, visualiser, optimiser.
    Pour mieux comprendre et maßtriser les procédés agro-alimentaires, les modÚles numériques deviennent un support de plus en plus utile. Cependant, la construction de ces modÚle reste encore une étape ardue pour de nombreuses raisons, principalement à cause de la diversité des mécanismes en jeu (physiques, chimiques, biologiques), des relations non linéaires entre variables, du couplages inter-échelles des phénomÚnes, de la diversité des connaissances expertes ou empiriques, de l'incertitude des mesures et des connaissances, ainsi que du coût et de la rareté des données expérimentales. Dans un contexte aussi complexe, il devient nécessaire d'avoir recours à des méthodes algorithmiques robustes, et en particulier à des heuristiques d'optimisation stochastique, car les méthodes classiques trouvent leurs limites. 
    Cet exposé est focalisé sur une approche développée dans l'équipe Malices de l'UMR GMPA en collaboration avec l'équipe AVIZ de l'INRIA, pour comprendre, organiser et structurer des données multi-dimensionnelles. L'exploration visuelle interactive, en formalisant l'exploration comme une tùche d'optimisation interactive prise en charge par un algorithme évolutionnaire, permet de proposer à des experts différents points de vues sur leurs jeux de données. Cette approche est utile à la fois au sein des premiÚres étapes de modélisation pour structurer les données brutes, mais aussi pour l'analyse du comportement de modÚles. 
     
  • Lundi 1er juin 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Marion Leclerc
    (INRA Micalis, Equipe PhylHom)
    Annulé
  • Lundi 18 mai 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Khashayar Pakdaman
    (U. Paris Diderot - CNRS, UMR 7592 Institut Jacques Monod, Biologie computationnelle et biomathématiques)
    Noise variability and synchronization

    Motivated by experiemental studies on the influence of noise on neuronal behavior, this presentation will review our work on the impact of stochastic variability on neuronal dynamics and synchronization of neuronal assemblies.

  • Lundi 13 avril 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Claire Rogel-Gaillard
    (INRA-AgroParisTech, UMR 1313 Génétique Animale et Biologie Intégrative)
    Sciences Animales Paris Saclay : pour qui, pourquoi, pour quoi faire?
  • Lundi 30 mars 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Guillaume Achaz
    (UPMC-CNRS UMR7138, IBPS ABI et CollĂšge de France SMILE)
    The strange case of the Standard Neutral Model of molecular evolution: when, what and why

    We will investigate the model that is commonly assumed throughout population genetics and more generally molecular evolution, the so-called Standard Neutral Model. I will start by describing some of its motivations, its underlying assumptions, its forward and backward perspectives, its usage and more importantly its limit. I will then illustrate the different points by biological examples and will argue that the use of a unique reference model may obscure our ability to apprehend correctly the evolution of life.

  • Lundi 16 mars 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Pierre RiviĂšre (1) et Olivier David (2)
    (RĂ©seau Semences Paysannes (1) et INRA UR1404 MaIAGE (2))
    Dispositifs expérimentaux pour la sélection décentralisée et participative sur le blé tendre

    Dans notre programme de sélection participative sur lé blé tendre, la sélection est décentralisée dans les environnements cibles and repose sur une collaboration étroite entre paysans, associations et équipes de recherche. Dans notre projet, chaque paysan conduit ses propres essais dans sa ferme et choisit les variétés qu'il souhaite semer. Cela génÚre des dispositifs expérimentaux trÚs déséquilibrés avec peu de degrés de liberté à la résiduelle dans chaque ferme et environ 95% de combinaisons variété x environnement manquante dans le réseau d'essais. Afin d'analyser ces données, nous avons développé deux modÚles hiérarchiques bayésiens afin de (1) réaliser des comparaisons de moyennes dans chaque ferme et (2) étudier les interactions variétés x environnements dans le réseau d'essais. Les deux modÚles apportent des résultats satisfaisants tant que le nombre de fermes est important dans le réseau et que chaque ferme a au moins une variété témoin répétée.

  • Lundi 2 mars 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Viet Chi Tran
    (CNRS - Université Lille 1, UMR 8524 Laboratoire Paul Painlevé)
    Un modÚle de propagation d'épidémie sur un graphe de configuration

    On considĂšre un modĂšle SIR, pour une maladie se propageant dans une population caractĂ©risĂ©e par une stucture sociale dĂ©crite par un graphe de configuration (Bollobas et Molloy-Reed). L'Ă©volution de l'Ă©pidĂ©mie peut-ĂȘtre rĂ©sumĂ©e par 3 Ă©quations Ă  valeurs mesures, d'oĂč l'on retrouve la description en 5 EDO proposĂ©e par Volz (2008). Nous expliquerons ensuite comment estimer les paramĂštres par ABC et comment faire une analyse de sensibilitĂ©.

  • Lundi 16 fĂ©vrier 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Sarah Lemler
    (UEVE - CNRS, UMR 8071 LAMME)
    Estimation pour les processus de comptage avec beaucoup de covariables

    Nous cherchons Ă  estimer l’intensitĂ© de sauts d’un processus de comptage en prĂ©sence d’un grand nombre de covariables. Nous proposons deux approches. D’abord, nous considĂ©rons une intensitĂ©
    non-paramĂ©trique et nous l’estimons par le meilleur modĂšle de Cox Ă©tant donnĂ© deux dictionnaires de fonctions. Le premier dictionnaire est utilisĂ© pour construire une approximation du logarithme
    du risque de base et le second pour approximer le risque relatif. Nous considérons une procédure Lasso, spécifique à la grande dimension, pour estimer simultanément les deux paramÚtres incon-
    nus du meilleur modĂšle de Cox approximant l’intensitĂ©. Nous prouvons des inĂ©galitĂ©s oracles non- asymptotiques pour l’estimateur Lasso obtenu.
    Dans une seconde partie, nous supposons que l’intensitĂ© satisfait un modĂšle de Cox. Nous proposons deux procĂ©dures en deux Ă©tapes pour estimer les paramĂštres inconnus du modĂšle de Cox. La premiĂšre Ă©tape est commune aux deux procĂ©dures, il s’agit d’estimer le paramĂštre de rĂ©gression en grande dimension via une procĂ©dure Lasso. Le risque de base est ensuite estimĂ© soit par sĂ©lection de modĂšles, soit par un estimateur Ă  noyau avec une fenĂȘtre choisie par la mĂ©thode de Goldenshluger et Lepski. Nous Ă©tablissons des inĂ©galitĂ©s oracles non-asymptotiques pour les deux estimateurs du risque de base ainsi obtenus. Nous menons une Ă©tude comparative de ces estimateurs sur des donnĂ©es simulĂ©es, et enfin, nous appliquons les procĂ©dures implĂ©mentĂ©es Ă  une base de donnĂ©es sur le cancer du sein.

  • Lundi 2 fĂ©vrier 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Charlotte Baey
    (Ecole Centrale Paris, laboratoire MAS)
    Modélisation de la variabilité inter-individuelle dans les modÚles de croissance de plantes

    Il existe une forte variabilitĂ© gĂ©nĂ©tique entre plantes, mĂȘme au sein de la mĂȘme variĂ©tĂ©, ce qui, combinĂ©e Ă  la variation locale des conditions climatiques dans le champ, peut conduire deux plantes voisines Ă  se dĂ©velopper de façon trĂšs diffĂ©rente. C’est l’une des raisons pour lesquelles les approches populationnelles dans les modĂšles de croissance de plantes suscitent un grand intĂ©rĂȘt. Nous proposons dans cette Ă©tude une extension du modĂšle individu-centrĂ© Greenlab Ă  l’échelle de la population dans le cas du colza, Ă  l’aide d’un modĂšle non linĂ©aire mixte. Deux variants stochastiques de l’algorithme EM (EspĂ©rance-Maximisation), le Monte-Carlo EM automatique (MCEM) et le SAEM seront comparĂ©s, en utilisant le fait que le modĂšle complet appartient Ă  la famille exponentielle.

  • Lundi 19 janvier 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Bogdan Mirauta
    (Biologie Computationnelle et Quantitative, UMR 7238 CNRS-UPMC; MaIAGE INRA)
    Transcriptome Analysis from High-Throughput Sequencing Count Data

    The most common RNA-Seq strategy consists of random shearing, amplification, and high-throughput sequencing, of the RNA fraction. Methods to analyze transcription level variations along the genome from the read count profiles generated by the is global RNA-Seq protocol are needed. We developed statistical approaches to estimate the local transcription levels and to identify transcript borders. The transcriptional landscape reconstruction relies on a state-space model to describe transcription level variations in terms of abrupt shifts and more progressive drifts. A new emission model is introduced to capture not only the read count variance inside a transcript but also its short-range autocorrelation and the fraction of positions with zero-counts. The estimation relies on a Sequential Monte Carlo algorithm, the Particle Gibbs.

  • Mardi 6 janvier 2015 - - Salle de rĂ©union 142, bĂątiment 210
    Rosemary Bailey
    (Queen Mary University of London, School of Mathematical Sciences)
    Designs for variety trials with very low replication

    In the early stages of testing new varieties, it is common that there are only small quantities of seed of many new varieties. In the UK (and some other countries with centuries of agriculture on the same land) variation within a field can be well represented by a division into blocks. Even when that is not the case, subsequent phases (such as testing for milling quality, or evaluation in a laboratory) have natural blocks, such as days or runs of a machine. I will discuss how to arrange the varieties in a block design when the average replication is less than two.

2014

  • Lundi 24 novembre 2014 - - Salle de rĂ©union 142, bĂątiment 210
    Maud Delattre
    (AgroParisTech-INRA, UMR 518 MIA)
    â–șTitre et rĂ©sumĂ© Ă  venir
  • Lundi 13 octobre 2014 - - Salle de rĂ©union 142, bĂątiment 210
    Simon Labarthe
    (INRA, MiaJ)
    â–ș Equations de rĂ©action diffusion en modĂ©lisation cardiaque et en dynamique de population.

    AprĂšs une brĂšve prĂ©sentation de modĂšles de dynamique de populations soumises au rĂ©chauffement climatique, j'aborderai de maniĂšre plus approfondie des questions de modĂ©lisation en cardiologie. J'illustrerai les mĂ©thodes et les interactions avec biologistes et cliniciens que j'ai pu mettre en Ɠuvre prĂ©cĂ©demment pour la modĂ©lisation des oreillettes et des veines pulmonaires. Les tissus auriculaires sont trĂšs fins : les modĂšles auriculaires surfaciques usuels tirent avantage de cette caractĂ©ristique.
    Cependant, des études cliniques ont montré que des événements électriques ont lieu dans l'épaisseur du tissu lors d'activité pathologique, ce que ne peuvent prendre en compte les modÚles surfaciques. Je présenterai un modÚle bisurfacique, de sa dérivation théorique par analyse asymptotique à son implémentation et à son utilisation pratique pour aborder des problématiques cliniques. Ces aspects applicatifs serviront à illustrer des aspects méthodologiques, plus susceptibles d'entrer en résonance avec les problématiques de recherche de l'unité.

  • Lundi 29 septembre 2014 - - Salle de rĂ©union 142, bĂątiment 210
    VĂ©ronique Cariou
    (ONIRIS, Nantes)
    â–ș Traitement de donnĂ©es mĂ©tabolomiques dans un contexte 3-voies. PrĂ©sentation de diffĂ©rentes approches en exploratoire versus en discrimination

    Dans le contexte de la chimiomĂ©trie, l’analyse des donnĂ©es mĂ©tabolomiques constitue un champ d’application en plein essor. La mĂ©tabolomique gĂ©nĂšre en effet de grands volumes de donnĂ©es pour lesquelles les techniques statistiques doivent ĂȘtre adaptĂ©es (volumĂ©trie, multi-colinĂ©aritĂ©, nombre d’individus trĂšs infĂ©rieur au nombre de variables, 
). Dans cet exposĂ©, nous nous intĂ©resserons Ă  une structure particuliĂšre de donnĂ©es mĂ©tabolomiques : celle des donnĂ©es trois voies. Ce type de donnĂ©es peut ĂȘtre directement issu des technologies employĂ©es (par exemple en chromatographie couplĂ©e Ă  la spectromĂ©trie de masse). La prise en compte d’une dimension temporelle (des prĂ©lĂšvements successifs pour un mĂȘme individu) engendre Ă©galement des donnĂ©es trois-voies, appelĂ©es aussi ternaires, tensorielles ou encore Ă  trois entrĂ©es. Nous prĂ©senterons tout d’abord ce type de donnĂ©es ainsi que les structures associĂ©es.
    En combinant une synthĂšse de travaux mĂ©thodologiques et des applications sur des jeux de donnĂ©es mĂ©tabolomiques, nous dĂ©velopperons ensuite diffĂ©rentes mĂ©thodes de traitement des donnĂ©es trois-voies. Nous nous placerons d’abord dans un contexte non supervisĂ© en montrant comment ces techniques gĂ©nĂ©ralisent l’analyse en composantes principales au cas des donnĂ©es trois-voies. Nous prĂ©senterons principalement deux modĂšles : Parafac et Tucker en dĂ©taillant Ă  la fois les critĂšres et les algorithmes. Dans un second temps, nous nous intĂ©resserons au cadre de la discrimination. Cette problĂ©matique est relativement rĂ©cente dans le cas de donnĂ©es trois voies. Parmi les travaux, nous dĂ©taillerons la NPLS proposĂ©e par R. Bro. Certaines limites de la NPLS dans le contexte des donnĂ©es mĂ©tabolomiques seront dĂ©gagĂ©es. Nous proposerons finalement une alternative relativement simple Ă  la NPLS, cette alternative reposant sur une approche exploratoire.

  • Vendredi 19 septembre 2014 - - AmphithĂ©Ăątre Jacques Poly, bĂątiment 440
    Plusieurs intervenants
    â–ș "Bucoliques, les MathĂ©matiques ?"

    Séminaire autour de l'apport des mathématiques en Agriculture, Alimentation et Environnement, à l'occasion du départ en retraite de Jean-Baptiste Denis.
    Pour plus de détail voir la page dédiée