Multilevel Modeling Using R

David Rogosa Sequoia 224, rag{AT}stat{DOT}stanford{DOT}edu Office hours:

Course web page: http://rogosateaching.com/stat196/

To see full course materials from Spring 2015 go here

From explorecourses

STATS 196A (EDUC 401D): Multilevel Modeling Using R Multilevel data analysis examples using R. Topics include: two-level nested data, growth curve modeling, generalized linear models for counts and categorical data, nonlinear models, three-level analyses. For more information, see course website: http://rogosateaching.com/stat196/ Terms: Spr | Units: 1 | Grading: Satisfactory/No Credit Instructors: Rogosa, D. (PI) STATS 196A | 1 units | Class # 47104 | Section 01 | Grading: Satisfactory/No Credit | WKS | Students enrolled: 3 03/28/2016 - 06/01/2016 Wed 3:30 PM - 5:20 PM at Sequoia Hall 200 with Rogosa, D. (PI) Instructors: Rogosa, D. (PI) Notes: Class meets on April 6, April 13, April 20, April 27, May 11.

For the 1-unit enrollment in this short course, students are expected to engage in the four presentation class sessions, and for the fifth session each student makes a short (~10 min) presentation of a relevant data analysis they have conducted.*Preliminary* Course ScheduleFive (2hr) mtgs W 3:30 - 5:20 on April 6 13 20 27, May 11 Sequoia Hall 200 Weeks 1 - 4. a. Introduction: Basic analyses for two-level nested data, normal models (UK Exam data) b. Additional two-level (normal) models: experimental designs (Dyestuff), longitudinal data (growth curves, sleepstudy), observational data (High School and Beyond) c. Generalized linear mixed models for counts and categorical outcomes d. Three-level analyses (nested data and longitudinal data) e. Specialized applications: regression diagnostics, power calculations and design, ecological inference, survival analysis, nonlinear functional forms, mediation analysis, propensity scores and matching, imputation, item response theory Week 5. Student presentations of multilevel data analyses

The Registrar does have a form (no-fee) for faculty, staff, post-docs: Application for Auditor or Permit to Attend (PTA) Status

Many of the example presented in this short course are described in Examples from Multilevel Software Comparative Reviews Douglas Bates. Code version of MlmSoftRev

R-package containing mlmRev data examples. Bates talk on mlmRev Additional examples in core package lme4

Another set of examples: lmer for SAS PROC MIXED Users Douglas Bates Department of Statistics University of Wisconsin Madison Data sets from SAS System for Mixed Models

Overviews and additional examples from Doug Bates:

lme4: Mixed-effects modeling with R February 17, 2010 Springer (chapters). An merged version of Bates book: lme4: Mixed-effects modeling with R January 11, 2010

R Journal intro Fitting linear mixed models in R Using the lme4 package Douglas Bates (pp.27-30)

Collection of all Doug Bates lme4 talks Fitting linear mixed-effects models using lme4,

HSB and growth curve examples in John Fox lme tutorial

Current version of R is version 3.2.4 (Very Secure Dishes) released 2016-03-10..(Rogosa is running 3.2.2) . For references and software: The R Project for Statistical Computing Closest download mirror is Berkeley

A new text (potentially) provides more infrastructure for this short course, but sadly it has many shortcomings. This text has free access at Stanford via crcnetbase.com Multilevel Modeling Using R

http://www.crcpress.com/product/isbn/9781466515857 Chapman & Hall/CRC Statistics in the Social and Behavioral Sciences Published:June 23, 2014 by CRC Press

Journal of Statistical Software Book Review

Book website, including data

1. Introductory Example. Goldstein Exam Data.

a. Introductory descriptive approaches for gender gap analysis (Smart First Year Student analyses using lmList, additional plots).

b. Various lmer analyses for gender gap.

Rogosa R-session basic plots models used

Stat 209 Gender gap data analysis. scanned class handout

2. Growth Curve Modeling exercise, Brain Volume Data Analysis. analyses from "Variation in longitudinal trajectories of regional brain volumes of healthy men and women (ages 10 to 85 years) measured with atlas-based parcellation of MRI" cartoon plot of Lateral Ventricles data; actual data plot of Lateral Ventricles data; development of lmer (mixed effect) growth models

1. Two-level nested (normal) data recap. Brief overview of HSB (High School and Beyond) analysis (from Stat209): plots and model.

a full single Bryk dataset (longform) (abbreviated) Rogosa R-session Bryk data plots, Rogosa R-session

A nice teaching document from Indiana that does HSB data with every known statistical package (including

2. Growth Curve models and analysis. Sleepstudy example (from Stat222).

Sleepstudy class handout, pdf scan [more Doug Bates Slides (pdf pages 8-28) ] Individual plots (frame-by-frame) Plot of straight-line fits

3. Data from designed experiments.

a. Dyestuff data, Bates book, Chapter 1 (sec 1.2, 1.3) Rogosa Dyestuff session

b. Penicillin data (also Pastes, ratbrain), Bates book, Chapter 2

4. Generalized Linear Mixed Models: counts and proportions.

nice overview: Generalized Linear Mixed Models from Encyclopedia of Statistics in Behavioral Science.

Dichotomous outcomes,

a. Respiratory clinical trial from HSAUR. lmList does logistic, introducing glmer lmList, glmer for respiration data (placebo group)

b. Contraception (Bangladesh) use from Bates review Rogosa R-session

c. Test scores (pass/fail outcome) from Ch 8, Multilevel Modeling Using R. Rogosa R-session

1. one-way anova, E(MS) recap

2. P-values, tests add-ons to lmer. Faraway text addendums: Inferential Methods for Linear Mixed Models

3. Residual plots, add-on regression diagnostics: packages

Rogosa session with Exam data (week 1) (ascii) resulting plots

4. Matrix Formulation for Mixed Effects Models (growth curves and nested data).

5. GLMM continued, dichotomous outcome

Contraception (Bangladesh) use from Bates review Rogosa R-session

6. Count outcome, GLMM poisson models.

Count data: Contagious bovine pleuropneumonia, data(cbpp) in lme4. Rogosa R-session herd plots

Another count data example from mlmRev package,

data(Mmmec): Malignant Melanoma Mortality in the European Community associated with the impact of UV radiation exposure. Rogosa glmer session

Measured outcomes

a. Achieve data from

b. example from mlmRev package, data(Chem97): Scores on A-level Chemistry in 1997. Rogosa R-session

Count outcome.

c. data(Mmmec): Malignant Melanoma Mortality in the European Community associated with the impact of UV radiation exposure.

A. Nonlinear Models.

B. Extensions of lme4 modeling:

Package

C. Multilevel Survival Analysis.

Mixed effects (Frailty) survival models. Package

Stat222 coxme handout plain text version

D. Mediation Analysis Package

E. Propensity Score and Matching Methods. Package

F. Ecological Inference R-package

G. Multiple Imputation. Package

H. Power calculations and Design;

Package pamm Title Power Analysis for Random Effects in Mixed Models

Package simr Title Power Analysis for Generalised Linear Mixed Models by Simulation.

longpower power and sample size for linear models of longitudinal data

Package MultiRR Title Bias, Precision, and Power for Multi-Level Random Regression

Student presentations, multilevel analysis examples.