# Markovian modelling and Bayesian learning, fall 2015

### Lecturer

### Scope

5 sp.

### Type

Advanced studies

### Prerequisites

Basic calculus, linear algebra, introductory course on probability and statistical inference are absolutely necessary. First course level knowledge on algebra, probability and inference will be recommendable for many parts of the course. Short course synopsis:

Part 1. Basic properties of discrete-time Markov chains (DTMCs), irreducibility, ergodicity, stationarity, invariant distributions,higher order Markov chains

Part 2. Statistical inference for DTMCs, maximum likelihood estimation, Bayesian estimation, inference about the order of a DTMC, full likelihood, model averaging, applications to clustering of DNA sequences

Part 3. Continuous-time Markov chains (CTMCs), basic properties, waiting time distributions, matrix forward equations, generator, absolute vs relative time, maximum likelihood estimation, unidentifiability of root position, applications to DNA and amino acid sequence evolution

Part 4. Hidden Markov models (HMMs). Basic properties, inference tasks related to smoothing, filtering and prediction, Derin's algorithm, factorization and recursion for HMM calculations, applications to classification and modeling of DNA sequences

Part 5. Markov random fields and pseudolikelihood inference, variable-order DTMCs (VOMs/VLMCs), sparse higher-order Markov chains (SMCs), inference for VOM/VLMC/SMC models, applications to sequence prediction and classification

### Lectures

Weeks 44-49, Tuesday 12-14 and Thursday 12-14 in Exactum room B120. Course end on December 4th.

Anonymous online chat for the course is available at http://presemo.helsinki.fi/mmbl2015

### Exercises

During weeks 45-49 there will be a weekly exercise session in Exactum room C129 on Thursday 14-16. The teacher responsible for the exercise sessions is Alberto Pessia (he has 'at' helsinki dot fi adress). Solutions to each week's exercises must be presented in the sesssion or at least 1 hour beforehand by email to the responsible teacher.

Exercises for week 45 are available here

Exercises for week 46 are available here

Exercises for week 47 are available here

Exercises for week 48 are available here

Exercises for week 49 are available here

### Exams

To gain the credits from this course, it is necessary to do at least 50% of the exercises and a home exam. Additional solved exercises will yield bonus points for the grade. The home exam will consist of a number of larger assignments that must be returned by March 1st 2016 to the lecturer. Home exam assignments are available here.

### Preliminary lecture diary

Week 44:

Tue

Log-linear expansions of Markov random field probabilities considered during the lecture are shown in detail in this document, three examples on why Markovian type models are useful: modeling the Ebola outbreak, HMM and tracking of objects (examples starting at 49:30), compression of DNA sequence data (backup of the paper if don't have access rights), some slides about modeling neutral evolution, Teaser trailer, Eye-opener on conditional probabilities and Bayes' theorem, basic properties of Markov chains. This excerpt from the HMM book by T. Koski is mainly used during the lecture and also this short excerpt on periodicity from the book of Isaacson & Madsen, Markov chains. For further illustrations and mathematical details on Markov chains, see the link to Sirl and Norris in Bibliography.

Thu

Basic properties of Markov chains continued. To get going with the basics of simulating Markov chains, you might find these Matlab codes useful.

Week 45:

Tue

Properties of Markov chains continued. Derivation of (4.10) on p 164 in the HMM book. Basics of ML and Bayesian learning, see this excerpt from the HMM book by T. Koski.

Thu

Statistical learning for DTMC's, see this excerpt from the HMM book by T. Koski. Also, this appendix from the HMM book is useful for refreshing details on various distributions.

Week 46:

Tue

Bayesian estimation of DTMC parameters, Continuous-time Markov chains (see the e-book by Koski).

Week 47:

Slides for inference about intractable CTMC models, Continuous-time Markov chains (CTMCs) continued, for an application of birth-death processes to modeling the 2014 Ebola epidemic see this paper, here is also an application of CTMCs to infectious disease epidemiology

basic properties of hidden Markov models, see: Ch. 10,Ch. 12,Ch. 13,Ch. 14 from the HMM book. An example of using HMM in classification. A biological example of the use of HMM is here. HMMs are also relevant for a multitude of engineering applications, such as dynamic tracking, an excellent technical review of this field by Arnaud Doucet is here, another excellent review by Cyrill Stachniss and an excellent short introduction by Bryan Minor is here.

Weeks 48-49:

Slides for introduction to pseudolikelihood inference for Markov random fields

A primer on Occham's razor and Bayesian model comparison for Markov chains, Information-theoretic book by D MacKay where Ch 28 contains a detailed explanation of the Occham's razor principle and Bayesian model comparison, Bayesian learning of the order of a DTMC, see also the use of information-theoretic criteria for choosing the model dimension as explained in this nice review paper, use of Markov chains to clustering of DNA sequences in metagenomics applications, paper1 (published in the journal Statistical Applications in Genetics and Molecular Biology), paper2

CTMS and HMMs continued, Variable Length Markov chains (see the article by Mächler & Buhlmann mentioned in the bibliography), Mixture Transition Markov models by Adrian Raftery, sparse Markov chains. Finite mixture models and EM-algorithm, Markov random fields: basics and learning, see this tutorial and this recent method based on Bayesian pseudolikelihood.

### Bibliography

Various references will be used during the course. The lecture diary will also include links to some additional materials. Parts of the following books will be considered:

Timo Koski. Hidden Markov models for bioinformatics. Kluwer, 2001.

Timo Koski & John M. Noble. Bayesian networks: An introduction. Wiley, 2009.

Timo Koski. Lectures at RNI on Probabilistic Models and Inference for Phylogenetics. Free e-book available here.

In addition, we will consider a number of articles & tutorials (articles not directly linked here are generally available form JSTOR collection or are otherwise online):

Braun, J.V. & Muller, H-G. Statistical methods for DNA sequence segmentation. Statistical Science, 13, 142-162, 1998.

Sirl, D. Markov Chains: An Introduction/Review. pdf.

Norris, J. Markov chains. CUP, see online resource.

Gu, L. Notes on Dirichlet distribution with relatives. This document provides a concise recapitulation of some of the central formulas that are needed in the exercises and assignments when doing Bayesian learning. More comprehensive derivations can be found in several books on Bayesian modeling, e.g. in Koski & Noble (2009), which is listed above.

Mächler, M. & Buhlmann, P. Variable length Markov chains: Methodology, computing and software. Journal of Computational and Graphical Statistics 13, 435-455, 2004. Preprint available here

Kass, R.E. & Raftery, A.E. Bayes factors. Journal of the American Statistical Association, 90, 773-795, 1995.

## Registration

Did you forget to register? What to do?

## Exercises

### Assignments

- Set 1
- Set 2

### Exercise classes

Group | Day | Time | Room | Instructor |
---|---|---|---|---|

1. |

## Course feedback

Course feedback can be given at any point during the course. Click here.