The textbook for the course is

- Madsen, H. and Thyregod, P. (2011) Introduction to General and Generalized Linear Models, Chapman and Hall/CRC, parts of chapter 5 and 6.

You may find it useful to also consult Faraway (2006) Extending the linear model with R, Chapman and Hall/CRC, chapter 2 and 8, which is a useful reference regarding more practical aspects of fitting mixed models. Pages 41-60 in Demidenko (2013) Mixed Models: Theory and Applications with R, Second Edition, give technical background on MLE and RMLE. Pages 9-49 in West et al (2014) Linear Mixed Models: A Practical Guide Using Statistical Software, Second Edition, also give background on linear mixed models and MLE.

The original reference regarding analysis of variance for orthogonal designs is Tjur, Tue (1984) Analysis of variance models in orthogonal designs, *International Statistical Review*, **52**, 33-65.

The course will be evaluated by active participation in solving exercises and presenting solutions of exercises. Mat8/Mat8tek/Mat8øk students must complete an additional miniproject related to Bayesian inference.

Some data sets used in the course. The orthodontic data set "Orthodont" is a part of the nlme R package. The wheezing data "Ohio" is a part of the geepack R package.

During the course you will work on and present several mini-projects:

**First miniproject**: ANOVA for split-plot
experiments - theory and practice.

**Second miniproject:** Linear mixed
model with AR(1) errors - parts 1-3. Along with the project description
comes some code
and more code.

The following will be updated during the course (*remember to refresh to get latest version of handouts*).

**1. Feb 2 8.15-12** Handouts (with exercises). Read M & T page 157-168.

**2. Feb 13 12.30-16** Handouts (with exercises). Read M & T page 169-171 and 175-182.

**3. Feb 16 8.15-12** Handouts (with exercises). We first consider practical implementation of ML and REML for linear mixed models. Secondly we consider a geometric approach to one-way anova. We also discuss the so-called hierarchical principle - try this R-code to get a better understanding of the importance of the hierarchical principle and the implications of different parametrizations. A few basic facts about orthogonal projections.

**4. Feb 20 12.30-16** We finish the treatment of one-way ANOVA and continue with two-way ANOVA for balanced designs and some examples of data analyses. Handouts (with exercises).

**5. Feb 23 8.15-12** We consider inference for anova-models Handouts (with exercises). Some R-code and data for the gene-expression example.

**6. Feb 27 12.30-16** We consider asymptotic inference for general linear mixed
models. Handouts. Some R-code
for a simulation study and parametric bootstrap. We also conduct a midway evaluation.

**7 - 10. (weeks 9-11)** Self-study. Groups work on and complete mini-projects.

**11. Mar 20 12.30-16** Groups present mini-projects.

**12. Mar 23 8.15-12** Prediction. Handouts. Read M & T page 171-174 and 182-186.

**13. Mar 27 12.30-16** Conditional independence, unnormalized densities and the Kalman-filter. Slides and tutorial.

**14. April 3 12.30-16** Bayesian inference. Slides. Read sections 6.1, 6.2 and 6.4 in Chapter 6 in M & T.

**15. April 10 12.30-16** Bayesian
inference. Slides.

**16. April 20 8.15-12** Presentation of last exercises and presentation of mini-projects for Mat8-oek-tek.

**Third miniproject:** Topics in Bayesian
statistics: choose either large sample methods in Bayesian inference (first 5 pages)
or choice of priors for variances in
hierarchical models. Here is a short introduction to the two options. Uddybende spørgsmål til 3. miniprojekt.

**17. May 17 8.15-12** Presentation of miniproject 3 (G5-112).

Last modified