# Markovian modelling and Bayesian learning, fall 2013

### 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. 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-50, Tuesday 12-14 and Thursday 12-14 in room B120. NB! No lectures on Nov 14. Course ends on Dec 10th.

### Exercises.

During weeks 45-49 there will be a weekly exercise session in room B120 on Thursdays 14-16. The teacher responsible for the exercise sessions is Elina Numminen (first.last@helsinki.fi). Solutions to each week's exercises must be presented to Elina in the sesssion or at least 1 hour beforehand by email.

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 May 1st 2014 to the lecturer. Home exam assignments are available here.

### Preliminary lecture diary

Week 44:

Tue

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).

Thu - no lecture

Week 47:

Continuous-time Markov chains (CTMCs) continued, 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:

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 (in press in the journal Statistical Applications in Genetics and Molecular Biology), paper2

Week 50: CTMS and HMMs continued, Variable Length Markov chains (see the article by Mächler & Buhlmann mentioned in the bibliography), sparse Markov chains. Finite mixture models and EM-algorithm,

### 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:

Baclawski, Kenneth. Introduction to probability with R. Chapman & Hall, 2008.

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.

Jordan, M.I. Graphical models. Statistical Science, 19, 140-155, 2004. Preprint available here

### Registration

Did you forget to register? What to do?

### Laskuharjoitukset

Ryhmä | Päivä | Aika | Paikka | Pitäjä |
---|---|---|---|---|

1. | Thursday | 14-16 | B120 | Elina Numminen (first dot last 'at' helsinki dot fi) |