Bayesian inference of biogeographical histories for hundreds of discrete areas
|
|
- Abner Gervase Walters
- 5 years ago
- Views:
Transcription
1 Bayesian inference of biogeographical histories for hundreds of discrete areas Michael Landis Nick Matzke Brian Moore John Huelsenbeck Evolu>on 06/23/13
2 Biogeography Every species has come into existence coincident both in space and 5me with a pre- exis5ng closely allied species. AR Wallace, 1855
3 Con>nental- scale biogeography (Mol Phylo Evol 2012) Supplemental Figure 1 a) b) ABCDEFGHI Abrocoma cinerea Abrocoma bennettii Octodontomys gliroides Octodon bridgesi Octodon lunatus Octodon degus Aconaemys porteri Aconaemys sagei Aconaemys fuscus Spalacopus cyanus Octomys mimax Pipanacoctomys aureus Tympanoctomys barrerae Ctenomys haigi Ctenomys boliviensis Ctenomys steinbachi Isothrix sinnamariensis Isothrix bistriata Isothrix barbarabrownae Makalata macrura Phyllomys brasiliensis Makalata didelphoides Phyllomys blainvillii Toromys grandis Echimys chrysurus Mesomys occultus Mesomys hispidus Lonchothrix emiliae Dactylomys dactylinus Dactylomys boliviensis Kannabateomys amblyonyx Proechimys roberti Proechimys longicaudatus Proechimys simonsi Proechimys quadruplicatus Hoplomys gymnurus Thrichomys apereoides Myocastor coypus Capromys pilorides Trinomys dimidiatus Trinomys iheringi Trinomys yonenagae Trinomys eliasi Trinomys paratus Trinomys setosus Clyomys laticeps Euryzygomatomys spinosus Dinomys branickii Lagostomus maximus Lagidium viscacia Chinchilla lanigera Sphiggurus melanurus Erethizon dorsatum Coendou bicolor Myoprocta acouchy Dasyprocta leporina Hydrochoerus hydrochaeris Kerodon rupestris Dolichotis patagonum Microcavia australis Cavia tschudii Cavia porcellus Cavia aperea Galea musteloides Cuniculus paca Cuniculus taczanowskii A B C E D H G F I Octodon degus Photo by José Cañas
4 Global- scale biogeography
5 For 8 areas
6 For 80 areas
7 For 800 areas
8 For 8 zillion areas
9 Why 9 areas? 86 occurrences used (Upham & PaYerson, 2012) a) b) A B C D E F G H I Abrocoma cinerea Abrocoma bennettii Octodontomys gliroides Octodon bridgesi Octodon lunatus Octodon degus Aconaemys porteri Aconaemys sagei Aconaemys fuscus Spalacopus cyanus Octomys mimax Pipanacoctomys aureus Tympanoctomys barrerae Ctenomys haigi Ctenomys boliviensis Ctenomys steinbachi Isothrix sinnamariensis Isothrix bistriata Isothrix barbarabrownae Makalata macrura Phyllomys brasiliensis Makalata didelphoides Phyllomys blainvillii Toromys grandis Echimys chrysurus Mesomys occultus Mesomys hispidus Lonchothrix emiliae Dactylomys dactylinus Dactylomys boliviensis Kannabateomys amblyonyx Proechimys roberti Proechimys longicaudatus Proechimys simonsi Proechimys quadruplicatus Hoplomys gymnurus Thrichomys apereoides Myocastor coypus Capromys pilorides Trinomys dimidiatus Trinomys iheringi Trinomys yonenagae Trinomys eliasi Trinomys paratus Trinomys setosus Clyomys laticeps Euryzygomatomys spinosus Dinomys branickii Lagostomus maximus Lagidium viscacia Chinchilla lanigera Sphiggurus melanurus Erethizon dorsatum Coendou bicolor Myoprocta acouchy Dasyprocta leporina Hydrochoerus hydrochaeris Kerodon rupestris Dolichotis patagonum Microcavia australis Cavia tschudii Cavia porcellus Cavia aperea Galea musteloides Cuniculus paca Cuniculus taczanowskii A C B D H E I G F 13,264 occurrences available (GBIF)
10 Transi>on between two ranges Founda>onal work: Ree et al. (Evolu5on 2005) Ree & Smith (Syst Biol 2008) Range Character Ancestral >me Observed & extant
11 Transi>on probability Q Instantaneous rate matrix Matrix exponen>a>on accounts for all intermediate events.
12 For few areas, no problem 3 areas
13 For more areas, Q explodes 3 areas 10 areas = Matrix exponen>a>on too slow for more than ten areas.
14 BayArea: Method for more areas Download BayArea: bayarea.googlecode.com Landis et al. (Syst Biol, in press)
15 BayArea: Method for more areas Inspired by Robinson et al. (Mol Biol Evol 2003) Landis et al. (Syst Biol, in press)
16 BayArea: Method for more areas Inspired by Robinson et al. (Mol Biol Evol 2003) Key concepts: 1. Sample biogeographic histories, H Landis et al. (Syst Biol, in press)
17 BayArea: Method for more areas Inspired by Robinson et al. (Mol Biol Evol 2003) Key concepts: 1. Sample biogeographic histories, H 2. Compute likelihood, L,H Landis et al. (Syst Biol, in press)
18 BayArea: Method for more areas Inspired by Robinson et al. (Mol Biol Evol 2003) Key concepts: 1. Sample biogeographic histories, L,H 2. Compute likelihood, 3. Approximate P (,H D) using H Markov chain Monte Carlo (MCMC) Landis et al. (Syst Biol, in press)
19 1. Sample biogeographic histories, H Nielsen (Syst Biol 2002) Landis et al. (Syst Biol, in press)
20 2. Compute likelihood, L,H Range evolu>on events from range : r = X r j re rt r i /r sum of rates leaving prob any event at >me prob next event is j L,H = product of event types & >mes over tree Landis et al. (Syst Biol, in press)
21 3. Approximate P (,H D) using MCMC high L,H low L,H P (,H D) Landis et al. (Syst Biol, in press)
22 Can we infer distance effects? Distance- dependent dispersal model Redistributes the rate of area gain Anywhere Nearby 0 ¼ ½ Collapses to independence model Simula>on: 600 areas, 50 replicates, 8 distances Landis et al. (Syst Biol, in press)
23 BayArea recovers true parameters Distance effects ¼ ½ Rate of area loss Rate of area gain Distance effects Landis et al. (Syst Biol, in press)
24 Bayes factors iden>fy distance effects 100% % of simula>ons supported 75% 5 50% 0 25% 5 Bayes factors support for distance model Favors M 0 Insubstantial Substantial Strong Very strong Decisive 0% ¼ ½ Landis et al. (Syst Biol, in press)
25 Malesian Rhododendron Vireya Re- analysis of Webb & Ree (2012) work 65 species, 20 areas Landis et al. (Syst Biol, in press)
26 Malesian Rhododendron Vireya Re- analysis of Webb & Ree (2012) work 65 species, 20 areas Wallace s Line Known dispersal barrier Vireya crossing? Landis et al. (Syst Biol, in press)
27 Malesian Rhododendron Vireya Re- analysis of Webb & Ree (2012) work 65 species, 20 areas Wallace s Line Known dispersal barrier Vireya crossing? Data from Brown et al. (J Biogeogr 2012) Webb & Ree (Chapter 8 in Bio5c Evolu5on and Environmental Change in Southeast Asia 2012) Landis et al. (Syst Biol, in press)
28 Distance mayers for Vireya dispersal Posterior Prior Rate of area loss Rate of area gain Distance Landis et al. (Syst Biol, in press)
29 B) Posterior of Branches: % of inferred range east of Wallace s Line W E Node maps: W Posterior probability ancestral of presence ranges per area E East of Wallace s Line West of Wallace s Line Wallace s Line & Lydekker s Line: 1 crossing West East Wallace s Line: 3+ crossings West East
30 Phylowood: biogeographic anima>ons
31 Future direc>ons Rate- modifiers for other traits/features Incorpora>ng on paleo- etc.- ical data Occupancy models to handle false absences Specia>on models (allopatry vs sympatry) Adding to RevBayes (easy to develop models)
32 Summary Allows hundreds of areas for analysis Joint posterior of parameters and ancestral ranges Simple distance- dependent dispersal model Efficient model tes>ng framework Open- source soqware available
33 Thanks! Ques>ons? BayArea Nick Matzke Brian Moore John Huelsenbeck Biogeography for many areas bayarea.google.code.com Phylowood Biogeographic anima>ons Trevor Bedford mlandis.github.com/phylowood Helpful folks Bas>en Boussau Tracy Heath Josh Schraiber Sebas>an Höhna Contact twiyer.com/landismj
34
35 Extra slides
36 Malesian paleogeography We assume constant geography, but E 110 E 120 E 130 E Volcanoes Highlands Land Lakes Shallow sea Carbonate platforms Deep sea Trenches a 100 E 110 E 120 E 130 E b 0 10 S 10 S 20 S 60 Mya Paleocene 40 Mya Late Eocene 20 S c d S 10 S 20 S 30 Mya Middle Oligocene 20 Mya Early Miocene 20 S e f S 10 S 20 S 10 Mya Late Miocene 5 Mya Early Pliocene 100 E 110 E 120 E 130 E 100 E 110 E 120 E 130 E Lohman et al. (2011) 20 S
37 Vireya results Assume Vireya root age is 55 Mya Ancestral range posterior Joint WL and LL crossing once ~40 Mya All other WL crossings < 15 Mya Plausible biogeographical scenario Single long distance dispersal event around 40 Mya As Sundi and Sahul Shelf converge, repeated short dispersals
38 Dispersal- ex>nc>on model R (a) Y i,y j = 8 0 if Y j,a =0 >< 1 if Y j,a =1 0 if Y i and Y j di er at more than one area >: 0 if Y j =(0, 0,...,0) iid, Jukes- Cantor, forbids ex>nc>on
39 Rate- modified dispersal model 8 0 if Y j,a =0 R (a) Y i,y j = >< 1 (Y i,y j,a, ) if Y j,a =1 0 if Y i and Y j di er at more than one area >: 0 if Y j =(0, 0,...,0) Per- area rate of gain depends on current biogeographical range.
40 Compute likelihood, L,H
41 Distance- dependent dispersal model 2, the numerator). (Y i =(1, 1, 0, 0)! Y j =(1, 1, 0, 1),a=4, )= d(g 1,G 4 ) + d(g 2,G 4 ) {z } Rate- modifier d(g 1,G 3 ) + d(g 2,G 3 ) {z } + d(g 1,G 4 ) + d(g 2,G 4 ) {z } Normaliza>on
42 BayArea recovers rate of area gain Posterior rate of area gain True rate True distance effects Landis et al. (Syst Biol, in press)
43 BayArea recovers rate of area loss Posterior rate of area loss True distance effects True rate Landis et al. (Syst Biol, in press)
44 BayArea recovers distance effects Posterior of distance effects True value True distance effects Landis et al. (Syst Biol, in press)
A primer on phylogenetic biogeography and DEC models. March 13, 2017 Michael Landis Bodega Bay Workshop Sunny California
A primer on phylogenetic biogeography and DEC models March 13, 2017 Michael Landis Bodega Bay Workshop Sunny California Epidemiology Air travel communities H3N2 Influenza virus samples 2000 2007 CE Flu,
More informationTransitions between Andean and Amazonian centers of endemism in the radiation of some arboreal rodents
Transitions between Andean and Amazonian centers of endemism in the radiation of some arboreal rodents Upham et al. Upham et al. BMC Evolutionary Biology 2013, 13:191 Upham et al. BMC Evolutionary Biology
More informationUnexpected phylogenetic relationships of the painted tree rat Callistomys pictus (Rodentia: Echimyidae)
Copyright jul-set 2014 do(s) autor(es). Publicado pela ESFA [on line] http://www.naturezaonline.com.br Loss AC, Moura RT, Leite YLR (2014) Unexpected phylogenetic relationships of the painted tree rat
More informationBULLAR VOLUME IN SUBTERRANEAN AND SURFACE-DWELLING CAVIOMORPH RODENTS
Journal of Mammalogy, 84(1):185 189, 2003 BULLAR VOLUME IN SUBTERRANEAN AND SURFACE-DWELLING CAVIOMORPH RODENTS CRISTIAN ERIC SCHLEICH* AND ALDO IVAN VASSALLO Departamento de Biología, Facultad de Ciencias
More informationThe history of South American octodontoid rodents and its contribution to evolutionary generalisations
5 The history of South American octodontoid rodents and its contribution to evolutionary generalisations diego h. verzi, cecilia c. morgan and a. itatí olivares Introduction The peculiar New World hystricomorph
More informationPhylogenetic Inference using RevBayes
Phylogenetic Inference using RevBayes Model section using Bayes factors Sebastian Höhna 1 Overview This tutorial demonstrates some general principles of Bayesian model comparison, which is based on estimating
More informationIntroduction ORIGINAL PAPER. Alicia Álvarez 1,3,4 & Michelle Arnal 2,3
J Mammal Evol (2015) 22:519 533 DOI 10.1007/s10914-015-9291-z ORIGINAL PAPER First Approach to the Paleobiology of Extinct Prospaniomys (Rodentia, Hystricognathi, Octodontoidea) Through Head Muscle Reconstruction
More informationUsing phylogenetics to estimate species divergence times... Basics and basic issues for Bayesian inference of divergence times (plus some digression)
Using phylogenetics to estimate species divergence times... More accurately... Basics and basic issues for Bayesian inference of divergence times (plus some digression) "A comparison of the structures
More informationConcepts and Methods in Molecular Divergence Time Estimation
Concepts and Methods in Molecular Divergence Time Estimation 26 November 2012 Prashant P. Sharma American Museum of Natural History Overview 1. Why do we date trees? 2. The molecular clock 3. Local clocks
More informationCarolina M Voloch 1, Julio F Vilela 2, Leticia Loss-Oliveira 2 and Carlos G Schrago 3*
Voloch et al. BMC Research Notes 2013, 6:160 RESEARCH ARTICLE Open Access Phylogeny and chronology of the major lineages of New World hystricognath rodents: insights on the biogeography of the Eocene/Oligocene
More informationPhylogenetics: Bayesian Phylogenetic Analysis. COMP Spring 2015 Luay Nakhleh, Rice University
Phylogenetics: Bayesian Phylogenetic Analysis COMP 571 - Spring 2015 Luay Nakhleh, Rice University Bayes Rule P(X = x Y = y) = P(X = x, Y = y) P(Y = y) = P(X = x)p(y = y X = x) P x P(X = x 0 )P(Y = y X
More informationWho was Bayes? Bayesian Phylogenetics. What is Bayes Theorem?
Who was Bayes? Bayesian Phylogenetics Bret Larget Departments of Botany and of Statistics University of Wisconsin Madison October 6, 2011 The Reverand Thomas Bayes was born in London in 1702. He was the
More informationBayesian Phylogenetics
Bayesian Phylogenetics Bret Larget Departments of Botany and of Statistics University of Wisconsin Madison October 6, 2011 Bayesian Phylogenetics 1 / 27 Who was Bayes? The Reverand Thomas Bayes was born
More informationG. B. DIAZ,* R. A. OJEDA* and E. L. REZENDE
Functional Ecology 2006 Renal morphology, phylogenetic history and desert Blackwell Publishing Ltd adaptation of South American hystricognath rodents G. B. DIAZ,* R. A. OJEDA* and E. L. REZENDE *Grupo
More informationArtificial Intelligence
ICS461 Fall 2010 Nancy E. Reed nreed@hawaii.edu 1 Lecture #14B Outline Inference in Bayesian Networks Exact inference by enumeration Exact inference by variable elimination Approximate inference by stochastic
More informationInfer relationships among three species: Outgroup:
Infer relationships among three species: Outgroup: Three possible trees (topologies): A C B A B C Model probability 1.0 Prior distribution Data (observations) probability 1.0 Posterior distribution Bayes
More informationBayesian Inference using Markov Chain Monte Carlo in Phylogenetic Studies
Bayesian Inference using Markov Chain Monte Carlo in Phylogenetic Studies 1 What is phylogeny? Essay written for the course in Markov Chains 2004 Torbjörn Karfunkel Phylogeny is the evolutionary development
More informationMarkov Chain Monte Carlo
Markov Chain Monte Carlo Recall: To compute the expectation E ( h(y ) ) we use the approximation E(h(Y )) 1 n n h(y ) t=1 with Y (1),..., Y (n) h(y). Thus our aim is to sample Y (1),..., Y (n) from f(y).
More informationAnnouncements. Inference. Mid-term. Inference by Enumeration. Reminder: Alarm Network. Introduction to Artificial Intelligence. V22.
Introduction to Artificial Intelligence V22.0472-001 Fall 2009 Lecture 15: Bayes Nets 3 Midterms graded Assignment 2 graded Announcements Rob Fergus Dept of Computer Science, Courant Institute, NYU Slides
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee and Andrew O. Finley 2 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry & Department
More information(5) Multi-parameter models - Gibbs sampling. ST440/540: Applied Bayesian Analysis
Summarizing a posterior Given the data and prior the posterior is determined Summarizing the posterior gives parameter estimates, intervals, and hypothesis tests Most of these computations are integrals
More informationHistorical biogeography models with dispersal probability as a func7on of distance(s)
Historical biogeography models with dispersal probability as a func7on of distance(s) Nicholas J. Matzke, Postdoctoral Fellow, NIMBioS (Na6onal Ins6tute of Mathema6cal, www.nimbios.org) SSE symposium:
More informationSome of these slides have been borrowed from Dr. Paul Lewis, Dr. Joe Felsenstein. Thanks!
Some of these slides have been borrowed from Dr. Paul Lewis, Dr. Joe Felsenstein. Thanks! Paul has many great tools for teaching phylogenetics at his web site: http://hydrodictyon.eeb.uconn.edu/people/plewis
More informationDiscrete & continuous characters: The threshold model
Discrete & continuous characters: The threshold model Discrete & continuous characters: the threshold model So far we have discussed continuous & discrete character models separately for estimating ancestral
More informationMarkov Networks.
Markov Networks www.biostat.wisc.edu/~dpage/cs760/ Goals for the lecture you should understand the following concepts Markov network syntax Markov network semantics Potential functions Partition function
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry & Department
More informationPhylogenetic Inference using RevBayes
Phylogenetic Inference using RevBayes Substitution Models Sebastian Höhna 1 Overview This tutorial demonstrates how to set up and perform analyses using common nucleotide substitution models. The substitution
More informationEstimating Evolutionary Trees. Phylogenetic Methods
Estimating Evolutionary Trees v if the data are consistent with infinite sites then all methods should yield the same tree v it gets more complicated when there is homoplasy, i.e., parallel or convergent
More informationMonte Carlo in Bayesian Statistics
Monte Carlo in Bayesian Statistics Matthew Thomas SAMBa - University of Bath m.l.thomas@bath.ac.uk December 4, 2014 Matthew Thomas (SAMBa) Monte Carlo in Bayesian Statistics December 4, 2014 1 / 16 Overview
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 7 Approximate
More informationThe Importance of Systema2cs and Rassenkreis. Reading: Willi Hennig
The Importance of Systema2cs and Rassenkreis Reading: Throughout the class so far we seen that the distribu2on of an organism is the result of its biological history as well as geologic and clima4c history
More informationBayesian Estimation of Input Output Tables for Russia
Bayesian Estimation of Input Output Tables for Russia Oleg Lugovoy (EDF, RANE) Andrey Polbin (RANE) Vladimir Potashnikov (RANE) WIOD Conference April 24, 2012 Groningen Outline Motivation Objectives Bayesian
More informationHow should we go about modeling this? Model parameters? Time Substitution rate Can we observe time or subst. rate? What can we observe?
How should we go about modeling this? gorilla GAAGTCCTTGAGAAATAAACTGCACACACTGG orangutan GGACTCCTTGAGAAATAAACTGCACACACTGG Model parameters? Time Substitution rate Can we observe time or subst. rate? What
More informationCS 188: Artificial Intelligence. Bayes Nets
CS 188: Artificial Intelligence Probabilistic Inference: Enumeration, Variable Elimination, Sampling Pieter Abbeel UC Berkeley Many slides over this course adapted from Dan Klein, Stuart Russell, Andrew
More informationBayesian phylogenetics. the one true tree? Bayesian phylogenetics
Bayesian phylogenetics the one true tree? the methods we ve learned so far try to get a single tree that best describes the data however, they admit that they don t search everywhere, and that it is difficult
More informationSTA414/2104 Statistical Methods for Machine Learning II
STA414/2104 Statistical Methods for Machine Learning II Murat A. Erdogdu & David Duvenaud Department of Computer Science Department of Statistical Sciences Lecture 3 Slide credits: Russ Salakhutdinov Announcements
More informationSupplemental Information Likelihood-based inference in isolation-by-distance models using the spatial distribution of low-frequency alleles
Supplemental Information Likelihood-based inference in isolation-by-distance models using the spatial distribution of low-frequency alleles John Novembre and Montgomery Slatkin Supplementary Methods To
More informationBayesian Networks. instructor: Matteo Pozzi. x 1. x 2. x 3 x 4. x 5. x 6. x 7. x 8. x 9. Lec : Urban Systems Modeling
12735: Urban Systems Modeling Lec. 09 Bayesian Networks instructor: Matteo Pozzi x 1 x 2 x 3 x 4 x 5 x 6 x 7 x 8 x 9 1 outline example of applications how to shape a problem as a BN complexity of the inference
More informationBayesian inference & Markov chain Monte Carlo. Note 1: Many slides for this lecture were kindly provided by Paul Lewis and Mark Holder
Bayesian inference & Markov chain Monte Carlo Note 1: Many slides for this lecture were kindly provided by Paul Lewis and Mark Holder Note 2: Paul Lewis has written nice software for demonstrating Markov
More informationInferring Complex DNA Substitution Processes on Phylogenies Using Uniformization and Data Augmentation
Syst Biol 55(2):259 269, 2006 Copyright c Society of Systematic Biologists ISSN: 1063-5157 print / 1076-836X online DOI: 101080/10635150500541599 Inferring Complex DNA Substitution Processes on Phylogenies
More informationPitfalls of Heterogeneous Processes for Phylogenetic Reconstruction
Pitfalls of Heterogeneous Processes for Phylogenetic Reconstruction Daniel Štefankovič Eric Vigoda June 30, 2006 Department of Computer Science, University of Rochester, Rochester, NY 14627, and Comenius
More informationBayesian Inference. Chapter 1. Introduction and basic concepts
Bayesian Inference Chapter 1. Introduction and basic concepts M. Concepción Ausín Department of Statistics Universidad Carlos III de Madrid Master in Business Administration and Quantitative Methods Master
More informationBayesian support is larger than bootstrap support in phylogenetic inference: a mathematical argument
Bayesian support is larger than bootstrap support in phylogenetic inference: a mathematical argument Tom Britton Bodil Svennblad Per Erixon Bengt Oxelman June 20, 2007 Abstract In phylogenetic inference
More informationDown by the Bayes, where the Watermelons Grow
Down by the Bayes, where the Watermelons Grow A Bayesian example using SAS SUAVe: Victoria SAS User Group Meeting November 21, 2017 Peter K. Ott, M.Sc., P.Stat. Strategic Analysis 1 Outline 1. Motivating
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 3 Linear
More informationApproximate Bayesian Computation: a simulation based approach to inference
Approximate Bayesian Computation: a simulation based approach to inference Richard Wilkinson Simon Tavaré 2 Department of Probability and Statistics University of Sheffield 2 Department of Applied Mathematics
More informationSea level and vertical motion of continents from dynamic earth models since the Late Cretaceous
Datashare 45: Sea level and vertical motion of continents from dynamic earth models since the Late Cretaceous Sonja Spasojevic and Michael Gurnis AAPG Bulletin, v. 96, no. 11 (November 2012), pp. 2037
More informationBMI/CS 776 Lecture 4. Colin Dewey
BMI/CS 776 Lecture 4 Colin Dewey 2007.02.01 Outline Common nucleotide substitution models Directed graphical models Ancestral sequence inference Poisson process continuous Markov process X t0 X t1 X t2
More informationAlgorithmic Methods Well-defined methodology Tree reconstruction those that are well-defined enough to be carried out by a computer. Felsenstein 2004,
Tracing the Evolution of Numerical Phylogenetics: History, Philosophy, and Significance Adam W. Ferguson Phylogenetic Systematics 26 January 2009 Inferring Phylogenies Historical endeavor Darwin- 1837
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee University of Minnesota July 20th, 2008 1 Bayesian Principles Classical statistics: model parameters are fixed and unknown. A Bayesian thinks of parameters
More informationBayes Nets III: Inference
1 Hal Daumé III (me@hal3.name) Bayes Nets III: Inference Hal Daumé III Computer Science University of Maryland me@hal3.name CS 421: Introduction to Artificial Intelligence 10 Apr 2012 Many slides courtesy
More informationCS 343: Artificial Intelligence
CS 343: Artificial Intelligence Bayes Nets: Sampling Prof. Scott Niekum The University of Texas at Austin [These slides based on those of Dan Klein and Pieter Abbeel for CS188 Intro to AI at UC Berkeley.
More informationBayesian Networks BY: MOHAMAD ALSABBAGH
Bayesian Networks BY: MOHAMAD ALSABBAGH Outlines Introduction Bayes Rule Bayesian Networks (BN) Representation Size of a Bayesian Network Inference via BN BN Learning Dynamic BN Introduction Conditional
More informationSTAT 499/962 Topics in Statistics Bayesian Inference and Decision Theory Jan 2018, Handout 01
STAT 499/962 Topics in Statistics Bayesian Inference and Decision Theory Jan 2018, Handout 01 Nasser Sadeghkhani a.sadeghkhani@queensu.ca There are two main schools to statistical inference: 1-frequentist
More informationRobust demographic inference from genomic and SNP data
Robust demographic inference from genomic and SNP data Laurent Excoffier Isabelle Duperret, Emilia Huerta-Sanchez, Matthieu Foll, Vitor Sousa, Isabel Alves Computational and Molecular Population Genetics
More informationBayesian Networks in Educational Assessment
Bayesian Networks in Educational Assessment Estimating Parameters with MCMC Bayesian Inference: Expanding Our Context Roy Levy Arizona State University Roy.Levy@asu.edu 2017 Roy Levy MCMC 1 MCMC 2 Posterior
More informationConsistent Downscaling of Seismic Inversions to Cornerpoint Flow Models SPE
Consistent Downscaling of Seismic Inversions to Cornerpoint Flow Models SPE 103268 Subhash Kalla LSU Christopher D. White LSU James S. Gunning CSIRO Michael E. Glinsky BHP-Billiton Contents Method overview
More informationBayesian Inference and MCMC
Bayesian Inference and MCMC Aryan Arbabi Partly based on MCMC slides from CSC412 Fall 2018 1 / 18 Bayesian Inference - Motivation Consider we have a data set D = {x 1,..., x n }. E.g each x i can be the
More informationReconstruire le passé biologique modèles, méthodes, performances, limites
Reconstruire le passé biologique modèles, méthodes, performances, limites Olivier Gascuel Centre de Bioinformatique, Biostatistique et Biologie Intégrative C3BI USR 3756 Institut Pasteur & CNRS Reconstruire
More informationIntroduction to Bayesian Learning
Course Information Introduction Introduction to Bayesian Learning Davide Bacciu Dipartimento di Informatica Università di Pisa bacciu@di.unipi.it Apprendimento Automatico: Fondamenti - A.A. 2016/2017 Outline
More informationProbabilistic Graphical Networks: Definitions and Basic Results
This document gives a cursory overview of Probabilistic Graphical Networks. The material has been gleaned from different sources. I make no claim to original authorship of this material. Bayesian Graphical
More informationApproximate Inference
Approximate Inference Simulation has a name: sampling Sampling is a hot topic in machine learning, and it s really simple Basic idea: Draw N samples from a sampling distribution S Compute an approximate
More informationBayesian Statistics: An Introduction
: An Introduction Tim Frasier Copyright Tim Frasier This work is licensed under the Creative Commons Attribution 4.0 International license. Click here for more information. Outline 1. Bayesian statistics,
More informationBayes Nets: Sampling
Bayes Nets: Sampling [These slides were created by Dan Klein and Pieter Abbeel for CS188 Intro to AI at UC Berkeley. All CS188 materials are available at http://ai.berkeley.edu.] Approximate Inference:
More informationBiol 206/306 Advanced Biostatistics Lab 12 Bayesian Inference
Biol 206/306 Advanced Biostatistics Lab 12 Bayesian Inference By Philip J. Bergmann 0. Laboratory Objectives 1. Learn what Bayes Theorem and Bayesian Inference are 2. Reinforce the properties of Bayesian
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods Tomas McKelvey and Lennart Svensson Signal Processing Group Department of Signals and Systems Chalmers University of Technology, Sweden November 26, 2012 Today s learning
More informationOutline. 1 Introduction. 2 Problem and solution. 3 Bayesian tracking model of group debris. 4 Simulation results. 5 Conclusions
Outline 1 Introduction 2 Problem and solution 3 Bayesian tracking model of group debris 4 Simulation results 5 Conclusions Problem The limited capability of the radar can t always satisfy the detection
More informationMCMC and Gibbs Sampling. Kayhan Batmanghelich
MCMC and Gibbs Sampling Kayhan Batmanghelich 1 Approaches to inference l Exact inference algorithms l l l The elimination algorithm Message-passing algorithm (sum-product, belief propagation) The junction
More informationSTA 4273H: Sta-s-cal Machine Learning
STA 4273H: Sta-s-cal Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 2 In our
More informationLecture 6 Phylogenetic Inference
Lecture 6 Phylogenetic Inference From Darwin s notebook in 1837 Charles Darwin Willi Hennig From The Origin in 1859 Cladistics Phylogenetic inference Willi Hennig, Cladistics 1. Clade, Monophyletic group,
More informationMachine Learning
Machine Learning 10-701 Tom M. Mitchell Machine Learning Department Carnegie Mellon University January 13, 2011 Today: The Big Picture Overfitting Review: probability Readings: Decision trees, overfiting
More informationBayesian Statistics. State University of New York at Buffalo. From the SelectedWorks of Joseph Lucke. Joseph F. Lucke
State University of New York at Buffalo From the SelectedWorks of Joseph Lucke 2009 Bayesian Statistics Joseph F. Lucke Available at: https://works.bepress.com/joseph_lucke/6/ Bayesian Statistics Joseph
More informationPitfalls of Heterogeneous Processes for Phylogenetic Reconstruction
Pitfalls of Heterogeneous Processes for Phylogenetic Reconstruction DANIEL ŠTEFANKOVIČ 1 AND ERIC VIGODA 2 1 Department of Computer Science, University of Rochester, Rochester, New York 14627, USA; and
More informationJournal of Arid Environments
Journal of Arid Environments 73 (2009) 173 181 Contents lists available at ScienceDirect Journal of Arid Environments journal homepage: www.elsevier.com/locate/jaridenv The mammals of the Monte Desert
More informationBayesian Inference. Anders Gorm Pedersen. Molecular Evolution Group Center for Biological Sequence Analysis Technical University of Denmark (DTU)
Bayesian Inference Anders Gorm Pedersen Molecular Evolution Group Center for Biological Sequence Analysis Technical University of Denmark (DTU) Background: Conditional probability A P (B A) = A,B P (A,
More informationBIG4: Biosystematics, informatics and genomics of the big 4 insect groups- training tomorrow s researchers and entrepreneurs
BIG4: Biosystematics, informatics and genomics of the big 4 insect groups- training tomorrow s researchers and entrepreneurs Kick-Off Meeting 14-18 September 2015 Copenhagen, Denmark This project has received
More informationRisk Estimation and Uncertainty Quantification by Markov Chain Monte Carlo Methods
Risk Estimation and Uncertainty Quantification by Markov Chain Monte Carlo Methods Konstantin Zuev Institute for Risk and Uncertainty University of Liverpool http://www.liv.ac.uk/risk-and-uncertainty/staff/k-zuev/
More informationProbabilistic assessment of regional climate change: a Bayesian approach to combining. predictions from multi-model ensembles
Probabilistic assessment of regional climate change: a Bayesian approach to combining predictions from multi-model ensembles Claudia Tebaldi (ESIG/NCAR), Richard L. Smith (UNC-Chapel Hill) Doug Nychka
More informationAnnouncements. CS 188: Artificial Intelligence Fall Causality? Example: Traffic. Topology Limits Distributions. Example: Reverse Traffic
CS 188: Artificial Intelligence Fall 2008 Lecture 16: Bayes Nets III 10/23/2008 Announcements Midterms graded, up on glookup, back Tuesday W4 also graded, back in sections / box Past homeworks in return
More informationBiol 206/306 Advanced Biostatistics Lab 12 Bayesian Inference Fall 2016
Biol 206/306 Advanced Biostatistics Lab 12 Bayesian Inference Fall 2016 By Philip J. Bergmann 0. Laboratory Objectives 1. Learn what Bayes Theorem and Bayesian Inference are 2. Reinforce the properties
More informationSession: Probabilistic methods in avalanche hazard assessment. Session moderator: Marco Uzielli Norwegian Geotechnical Institute
Session: Probabilistic methods in avalanche hazard assessment Session moderator: Marco Uzielli Norwegian Geotechnical Institute Session plan Title Presenter Affiliation A probabilistic framework for avalanche
More informationAn Introduction to Reversible Jump MCMC for Bayesian Networks, with Application
An Introduction to Reversible Jump MCMC for Bayesian Networks, with Application, CleverSet, Inc. STARMAP/DAMARS Conference Page 1 The research described in this presentation has been funded by the U.S.
More informationBayesian Inference. Chapter 4: Regression and Hierarchical Models
Bayesian Inference Chapter 4: Regression and Hierarchical Models Conchi Ausín and Mike Wiper Department of Statistics Universidad Carlos III de Madrid Advanced Statistics and Data Mining Summer School
More informationAn example to illustrate frequentist and Bayesian approches
Frequentist_Bayesian_Eample An eample to illustrate frequentist and Bayesian approches This is a trivial eample that illustrates the fundamentally different points of view of the frequentist and Bayesian
More informationA Bayesian Approach to Phylogenetics
A Bayesian Approach to Phylogenetics Niklas Wahlberg Based largely on slides by Paul Lewis (www.eeb.uconn.edu) An Introduction to Bayesian Phylogenetics Bayesian inference in general Markov chain Monte
More informationUnfolding techniques for neutron spectrometry
Uncertainty Assessment in Computational Dosimetry: A comparison of Approaches Unfolding techniques for neutron spectrometry Physikalisch-Technische Bundesanstalt Braunschweig, Germany Contents of the talk
More informationA Probabilistic Framework for solving Inverse Problems. Lambros S. Katafygiotis, Ph.D.
A Probabilistic Framework for solving Inverse Problems Lambros S. Katafygiotis, Ph.D. OUTLINE Introduction to basic concepts of Bayesian Statistics Inverse Problems in Civil Engineering Probabilistic Model
More informationMarkov Networks. l Like Bayes Nets. l Graphical model that describes joint probability distribution using tables (AKA potentials)
Markov Networks l Like Bayes Nets l Graphical model that describes joint probability distribution using tables (AKA potentials) l Nodes are random variables l Labels are outcomes over the variables Markov
More informationOutline. CSE 573: Artificial Intelligence Autumn Agent. Partial Observability. Markov Decision Process (MDP) 10/31/2012
CSE 573: Artificial Intelligence Autumn 2012 Reasoning about Uncertainty & Hidden Markov Models Daniel Weld Many slides adapted from Dan Klein, Stuart Russell, Andrew Moore & Luke Zettlemoyer 1 Outline
More informationStatistical Inference
Statistical Inference Robert L. Wolpert Institute of Statistics and Decision Sciences Duke University, Durham, NC, USA Spring, 2006 1. DeGroot 1973 In (DeGroot 1973), Morrie DeGroot considers testing the
More informationMCMC: Markov Chain Monte Carlo
I529: Machine Learning in Bioinformatics (Spring 2013) MCMC: Markov Chain Monte Carlo Yuzhen Ye School of Informatics and Computing Indiana University, Bloomington Spring 2013 Contents Review of Markov
More informationRapid evolution of the cerebellum in humans and other great apes
Rapid evolution of the cerebellum in humans and other great apes Article Accepted Version Barton, R. A. and Venditti, C. (2014) Rapid evolution of the cerebellum in humans and other great apes. Current
More informationCSL302/612 Artificial Intelligence End-Semester Exam 120 Minutes
CSL302/612 Artificial Intelligence End-Semester Exam 120 Minutes Name: Roll Number: Please read the following instructions carefully Ø Calculators are allowed. However, laptops or mobile phones are not
More informationEvolution of morphological adaptations for digging in living and extinct ctenomyid and octodontid rodents
Biological Journal of the Linnean Society, 2008, 95, 267 283. With 8 figures Evolution of morphological adaptations for digging in living and extinct ctenomyid and octodontid rodents ENRIQUE P. LESSA 1
More informationCS 188: Artificial Intelligence Spring Announcements
CS 188: Artificial Intelligence Spring 2011 Lecture 16: Bayes Nets IV Inference 3/28/2011 Pieter Abbeel UC Berkeley Many slides over this course adapted from Dan Klein, Stuart Russell, Andrew Moore Announcements
More informationMarkov Networks. l Like Bayes Nets. l Graph model that describes joint probability distribution using tables (AKA potentials)
Markov Networks l Like Bayes Nets l Graph model that describes joint probability distribution using tables (AKA potentials) l Nodes are random variables l Labels are outcomes over the variables Markov
More informationBayesian inference for multivariate extreme value distributions
Bayesian inference for multivariate extreme value distributions Sebastian Engelke Clément Dombry, Marco Oesting Toronto, Fields Institute, May 4th, 2016 Main motivation For a parametric model Z F θ of
More informationEE562 ARTIFICIAL INTELLIGENCE FOR ENGINEERS
EE562 ARTIFICIAL INTELLIGENCE FOR ENGINEERS Lecture 16, 6/1/2005 University of Washington, Department of Electrical Engineering Spring 2005 Instructor: Professor Jeff A. Bilmes Uncertainty & Bayesian Networks
More informationSome slides from Carlos Guestrin, Luke Zettlemoyer & K Gajos 2
Logistics CSE 446: Point Estimation Winter 2012 PS2 out shortly Dan Weld Some slides from Carlos Guestrin, Luke Zettlemoyer & K Gajos 2 Last Time Random variables, distributions Marginal, joint & conditional
More informationThe University of Auckland Applied Mathematics Bayesian Methods for Inverse Problems : why and how Colin Fox Tiangang Cui, Mike O Sullivan (Auckland),
The University of Auckland Applied Mathematics Bayesian Methods for Inverse Problems : why and how Colin Fox Tiangang Cui, Mike O Sullivan (Auckland), Geoff Nicholls (Statistics, Oxford) fox@math.auckland.ac.nz
More information