# Difference between revisions of "Paul Lewis"

Line 30: | Line 30: | ||

==== C++ Bayesian Phylogenetics Software Tutorial ==== | ==== C++ Bayesian Phylogenetics Software Tutorial ==== | ||

− | <span style="color:red;font-weight:bold">New!</span> If you took John Huelsenbeck's morning C++ class, you may wish to continue learning by working through this [http://phylogeny.uconn.edu/phylogenetic-software-development-tutorial/|C++ Bayesian phylogenetics tutorial]. This tutorial starts by teaching you how to set up your Windows or Mac laptop for compiling C++ programs and then explaining each step until you have a functioning Bayesian MCMC sampler for inferring phylogenies or estimating parameters of evolutionary models. | + | <span style="color:red;font-weight:bold">New!</span> If you took John Huelsenbeck's morning C++ class, you may wish to continue learning by working through this [http://phylogeny.uconn.edu/phylogenetic-software-development-tutorial/ | C++ Bayesian phylogenetics tutorial]. This tutorial starts by teaching you how to set up your Windows or Mac laptop for compiling C++ programs and then explaining each step until you have a functioning Bayesian MCMC sampler for inferring phylogenies or estimating parameters of evolutionary models. |

==== The Chain Swapping Question ==== | ==== The Chain Swapping Question ==== |

## Revision as of 07:56, 28 July 2017

## Contents

### Time at MBL

I will be present from July 20 to July 30.

### Contact Info

E-mail: paul.lewis@uconn.edu

Address: Department of Ecology and Evolutionary Biology, University of Connecticut, 75 N. Eagleville Road, Unit 3043, Storrs, CT 06269-3043 U.S.A.

### Institutional home page

### Lecture Materials

#### Likelihood lecture (22 July 2017)

New! The link below has been updated with the 2017 version of the lecture.

- Likelihood lecture slides (89 pages, 3.5 MB, PDF format, uploaded 22 July 2017)

Here are a few slides that I use to explain how the Poisson probability formula arises (they help explain how the constant e sneaks into the JC69 transition probability formula):

- Poisson slides (7 pages, 516 KB, PDF format)

#### Bayesian lecture (23 July 2017)

New! The link below has been updated with the 2017 version of the lecture.

- Bayesian lecture slides (100 pages, 2.7 MB, PDF format, uploaded 23 July 2017)

MCMC Robot used in Bayesian lecture:

- free iOS or Windows app available from the MCMC Robot web site

New! The version I used in lecture is now available at http://phylogeny.uconn.edu/mcmc-robot/. This version still has a few glitches, but you should at least be able to replicate everything I did in the lecture without running into problems. Click outside an existing hill to move the starting point for a chain, and drag in an area outside a hill to create a new hill. Clicking on a hill to select it should allow you to delete it using the delete key. (It is currently easier to delete and recreate a hill than it is to change the diameter of an existing hill.)

#### C++ Bayesian Phylogenetics Software Tutorial

New! If you took John Huelsenbeck's morning C++ class, you may wish to continue learning by working through this | C++ Bayesian phylogenetics tutorial. This tutorial starts by teaching you how to set up your Windows or Mac laptop for compiling C++ programs and then explaining each step until you have a functioning Bayesian MCMC sampler for inferring phylogenies or estimating parameters of evolutionary models.

#### The Chain Swapping Question

During my demonstrations of the MCMC robot app, the question that is always asked is the following: Why does the swapping pattern of the cold chain appear concentrated between two of the three hills in the landscape; why is swapping not uniform and symmetrical between all three hills? The top figure below illustrates the pattern, with large magenta dots showing where each of the 3 chains are at the current time.

As you can see, swaps seem to avoid one pair of hills in the landscape. One key clue is also shown by this (top) figure: the 3 robots are not evenly distributed with 1 per hill. Instead, one hill currently has 2 robots, and another hills has just 1, and the third hill has 0 robots. While this configuration of robots persists, all swaps will be between the 2 hills that have robots, and a "rut" between these hills will develop. Eventually, one of the robots will manage to get over to the hill with no robots (note that this is hard because swapping does not help with this process), and swaps can then begin again to accrue to this third hill.

About 1/5 of the time, the 3 robots should all be on different hills (see the lower figure), and thus only 1/5 of the time should we expect to see symmetrical swapping.

Thus, we should expect ruts to develop, and it is not a problem because even in the top figure the cold chain robot is easily able to move between all three hills.