36-463/663: Hierarchical Linear Models
|
|
- Kathryn King
- 5 years ago
- Views:
Transcription
1 36-463/663: Hierarchical Linear Models Lmer model selection and residuals Brian Junker 132E Baker Hall 1 Outline The London Schools Data (again!) A nice random-intercepts, random-slopes model Residuals in MLM s Marginal residuals Conditional residuals Random effects residuals Variable selection Overall fit statistics Simulation-based methods (save for later!) AIC, BIC, DIC and all that An improved model for the London Schools data 2
2 The London Schools Data Student ( ) Gender (0=Female, 1=Male), per student VR = verbal reasoning level (High/Med/Low) LRT = London Reading test (at beginning of year) Y = end-of-year test School (1..38) School.gender (All.Boy, All.Girl, Mixed) School.denom (Other,CofE,RomCath,State) 3 London Schools Data The dotted lineis the pooledregression (ignoring schools) The solid linesare unpooledregressions (separate for each school) The solid lines look like a random sample of lines, with mean the solid line! Unpooled Pooled Note these plots were made with an earlier version of ggplot2. An alternative approach with the help of gridextra is shown in the R handout. 4
3 The London Schools Data This suggests a model like or, as a variance components model As an R model this would be Y ~ 1 + LRT + (1 + LRT school) 5 The London Schools Data > display(lmer.1 + <- lmer(y ~ 1 + LRT + + (1 + LRT school),data=school.frame)) Unpooled Pooled mlm alphas mlm betas coef.est coef.se (Intercept) LRT Error terms: Groups Name Std.Dev. Corr school (Int) 0.23 LRT Residual number of obs: 1978, groups: school, 38 AIC = , DIC = deviance =
4 Residuals In ordinary linear regression (and glm s) the residuals are easy to think about: E[y i ] = g -1 (X i β) r i = y i E[y i ] Multi-level models pose a couple of challenges 7 Residuals in Multi-Level Models Level 1 Level 2 Where are they? Level 1? Level 2? Some combination? What are they? The α s are random draws, so does the following make sense? E[y ij ] = α 0j + α 1j LRT ij?? r ij = y ij -E[y ij ]?? 8
5 Residuals in Multi-Level Models The variance components version of the model could be re-expressed in matrix form as i.e. Laird & Ware (1982, Biometrics) 9 Residuals in Multi-Level Models Given the Laird-Ware form, can formulate 3 different kinds of residuals: Marginal residuals: Conditional residuals: Random effects: In practice, estimate with, the MLE, and estimate with Although there is only one per group, there are as many s as there are observations. Nobre& Singer (2007, Biometrical Journal) 10
6 Residuals in Multi-Level Models Marginal residuals: Should be mean 0, but may show grouping structure May not be homoskedastic! Good for checking fixed effects, just like linear regr. Conditional residuals: Should be mean zero with no grouping structure Should be homoskedastic! Good for checking normality of ǫ, outliers Random effects: Will generally not be mean-zero May not be be homoskedastic! Good for checking normality of, outliers 11 Residuals in the London Schools Data > str(fixef(lmer.1)) > beta0 <- fixef(lmer.1)[1] > beta1 <- fixef(lmer.1)[2] > str(ranef(lmer.1)) > eta <- ranef(lmer.1)$school > attach(school.frame) > X <- cbind(1,lrt) > blocks <- lapply(split(x,school), + function(x){matrix(x,ncol=2)}) > J <- length(blocks) > n <- dim(school.frame)[1] > Z <- matrix(0,nrow=n,ncol=j*2) > row <- 1 > for (j in 1:J) { + col <- 2*j + nj <- dim(blocks[[j]])[1] + Z[row:(row+nj-1),c(col-1,col)] <- + blocks[[j]] + row <- row + nj + } > beta <- rbind(beta0,beta1) > # so beta is a column vector > eta <- c(t(eta)) > # so eta is a column vector > resid.marg <- Y - X%*%beta > resid.cond <- Y - X%*%beta - Z%*%eta > resid.reff <- Z%*%eta 12
7 Residuals in the London Schools Data Marginal residuals look pretty good Conditional residuals Marginal Residuals Conditional Residuals Index Index look pretty good Rand Effect residuals Random Effects look weird Index 13 Residuals in the London Schools Data Marginal residuals plotted by school No noticeable patterns Nice set of residuals 14
8 Residuals in the London Schools Data Conditional residuals plotted by school No noticeable patterns Nice set of residuals 15 Residuals in the London Schools Data Rand Effect residuals plotted by school The scaleof these residuals is smaller! A few schools show noticabledeviation from zero We do not expect mean-zero, but the BLUP estimates should cluster around a mean 16
9 Standardized Residuals We were looking at raw residuals You could standardize by dividing the raw residuals by the square root of the variance components, for example Quick and dirty: divide by their sample SD! If there are severe outliers, the sample SD will be inflated and this won t work well If we have standardized residuals, outlier detection is a bit easier. There are simulation-based methods for detecting outliers as well 17 Residuals: Practical Advice Looking at some residuals is better than looking at none. In many MLM s, marginal and conditional residuals can be used roughly as you would with ordinary linear regression It is worthwhile to plot residuals again the group/cluster indicators To identify and fix problems, plot residuals against other variables (within and/or across clusters), try transformations, etc residuals(lmer.1)gives you the conditional residuals! 18
10 The London Schools Data Variable Selection How can we improve the model? We have a bunch of other variables lying around: Unit-level (student): Gender, VR Group-level (school): School.denom, School.gender Which ones to include? Fixed effects or random effects? Interactions? Etc. 19 Comparing models Two basic tools Overall Fit Indices Likelihood ratio test / Deviance statistics AIC, BIC, DIC and all that Simulation-based Checks Fake data methods Posterior predictive checks 20
11 Overall fit indices Pro s: Quick-and-dirty (a drop of 2 is marginally interesting, a drop of 3 is interesting, a drop of 10 is a big deal) Tells you whether one model fits better overall than another Con s: Doesn t tell whether the best model really fits If a model doesn t fit, doesn t tell you what went wrong We will discuss these today 21 Simulation-based checks Con s The method is not automatic You have to choose a feature of the data or model that you want to check You have to program the simulation Pro s You can focus on particular features of the model or data that you are worried about You can see what went wrong with the model, and sometimes fix it We will discuss these later in the semester 22
12 Overall Fit Indices Basic tool: deviance of the model M Direct use: If M 0 is nested in M 1, then approx: Used to compare ordinary linear models, generalized linear models, etc. 23 What if models are not nested? The Chi-squared theory for D(M) is gone. Other considerations (prediction accuracy, Bayesian decis. theory, ) lead to penalized deviance measures: D P (M) = D(M) + P(#parameters) Three common cases (let k = #parameters) AIC = D(M) + 2 k BIC = D(M) + k log(n) DIC = mean sim (D(M)) + 2 k D D(M) + 2 k D k D is estimated#parameters 24
13 Why are we fudging about k vs k D? In ordinary linear models(and glm s), k = #X s. What about σ 2? what about over-dispersion parameter? In both M 1 and M 0, so df 1 df 0 is the same either way). In multilevel modelsit s not so clear: one-way ANOVA with J cells (df=j) fitting grand mean only (df=1) 1 k D J, depending on size of τ Variable Selection: Practical Advice Start with multilevel model that represents your initial guesses about group structure in the data Do variable selection on all the fixed effects first, using AIC or DIC (or if you prefer, BIC) AIC will result in bigger modelsthat predict better BIC will result in smaller models that interpret better DIC will result in models between AIC and BIC sizes Deviance/LRT only valid if models have same random effects, and are nested Then go back and use AIC or DIC (or BIC) to do selection on random effects 26
14 Digression: REML vs MLE In order to properly use AIC, BIC, or even LR tests, must calculate the true maximum log-likelihood. By default, lmer() calculates estimates known as REML estimates. Less biased estimates than true MLE s, especially for variance components Essentially, calculates fixed effect estimates by leastsquares, and calculates variance components from marginal residuals REML estimates do not maximize the true likelihood For AIC, BIC, LRT, etc., need to refit using MLEs! 27 Digression: REML vs MLE Can use lmer() or update() function to get MLE fit lmer.1 <-lmer(y ~ 1 + LRT + (1 + LRT school), data=school.frame, REML=F) lmer.1 <-update(lmer.1,. ~., REML=F) Can produce substantial differences in likelihood Use AIC(), BIC(), loglik() to extract these values directly from fitted model Not valid to compare models with these numbers OK to compare models with these numbers ANOVA() always refits using MLEs so that comparisons are valid 28
15 Back to London Schools Data > names(tmp) # main variabes in school.frame # [1] "Y" "LRT" "Gender" "School.gender" # [5] "School.denom" "VR" > lmer.2 <- update(lmer.1,. ~. + Gender) > anova(lmer.1,lmer.2) # refitting model(s) with ML (instead of REML) # Df AIC BIC loglik deviance Chisq Chi Df Pr(>Chisq) # lmer # lmer # --> AIC, BIC prefer lmer.2 > lmer.3 <- update(lmer.2,. ~. + School.gender) > anova(lmer.2,lmer.3) # refitting model(s) with ML (instead of REML) # Df AIC BIC loglik deviance Chisq Chi Df Pr(>Chisq) # lmer # lmer ** # --> AIC prefers lmer.3, BIC prefers lmer.2 Etc! 29 London Schools Data Tried Gender, School.gender, School.denom, and VR stepwise as fixed effects, found only Gender and VR reduce AIC by more than 3, so keep them. Gender and VR. Trying to convert Gender and VR to random effects does not improve AIC enough to keep them, so the final model we obtain is > lmer.5 <- lmer(y LRT + Gender + VR + + (1 + LRT school), data=school.frame)) 30
16 London Schools Data final model > display(lmer.5) lmer(formula = Y ~ LRT + Gender + VR + (1 + LRT school), data = school.frame) coef.est coef.se (Intercept) LRT Gender VRLow VRMed Error terms: Groups Name Std.Dev. Corr school (Intercept) 0.23 LRT Residual number of obs: 1978, groups: school, 38 AIC = , DIC = deviance = > anova(lmer.1,lmer.5) refitting model(s) with ML (instead of REML) Df AIC BIC loglik lmer lmer Some Automatic & Exact Methods There are a number of R packages that will do variable selection for lmer models, including: LMERConvenienceFunctionsautomates backwards selection of fixed effects and forward selection of random effects, using AIC, BIC, etc. fitlmer.fnc() is general-purpose function for this RLRsimprovides simulation-based exact likelihood ratio tests for random effects exactlrt() performs exact LRT test for true ML fits exactrlrt() performs exact LRT test for REML fits 32
17 Automated Variable Selection > library(lmerconveniencefunctions) # for fitlmer.fnc() function... # start with a "big fixed effects" model > lmer.10 <- lmer(y ~ LRT + VR + Gender + School.gender +School.denom + + (1+LRT school), data=school.frame) > bic_best.11 <- fitlmer.fnc(lmer.10, + ran.effects=c("(school.gender school)", + "(School.denom school)"),method="bic") fitlmer.fnc: > anova(lmer.5,lmer.10,lmer.11) 1. Backwards elimination of F.E s refitting model(s) with ML (instead of REML) 2. Forward selection of R.E. s Data: school.frame 3. Backwards elimination of F.E. s Models: lmer.5: Y ~ LRT + Gender + VR + (1 + LRT school) lmer.11: Y ~ LRT + VR + Gender + (1 + LRT school) lmer.10: Y ~ LRT + VR + Gender + School.gender + School.denom + (1 + LRT lmer.10: school) Df AIC BIC loglik deviance Chisq Chi Df Pr(>Chisq) lmer lmer <2e-16 *** lmer Exact Test of Random Effect.. library(rlrsim) m0 <- lmer(y ~ LRT + VR + Gender + (1 school), data=school.frame) lmer.11a <- lmer(y ~ LRT + VR + Gender + (1 school) + (0 + LRT school), data=school.frame) # need indep rand effects for RLRsim... lmer.lrt.only <- lmer(y ~ LRT + VR + Gender + (0 + LRT school), data=school.frame) formula(m0) # formula under H0: no random slopes for LRT formula(lmer.11a) # model under HA: yes random slopes for LRT formula(lmer.lrt.only) # model with *only* random slopes for LRT exactrlrt(lmer.lrt.only,lmer.11a,m0) # simulated finite sample distribution of RLRT. # # (p-value based on simulated values) # # data: # RLRT = , p-value =
18 Summary The London Schools Data (again!) A nice random-intercepts, random-slopes model Residuals in MLM s Marginal residuals Conditional residuals Random effects residuals Variable selection Overall fit statistics Simulation-based methods (save for later!) AIC, BIC, DIC and all that An improved model for the London Schools data 35
36-463/663: Multilevel & Hierarchical Models
36-463/663: Multilevel & Hierarchical Models (P)review: in-class midterm Brian Junker 132E Baker Hall brian@stat.cmu.edu 1 In-class midterm Closed book, closed notes, closed electronics (otherwise I have
More informationA brief introduction to mixed models
A brief introduction to mixed models University of Gothenburg Gothenburg April 6, 2017 Outline An introduction to mixed models based on a few examples: Definition of standard mixed models. Parameter estimation.
More informationIntroduction to Within-Person Analysis and RM ANOVA
Introduction to Within-Person Analysis and RM ANOVA Today s Class: From between-person to within-person ANOVAs for longitudinal data Variance model comparisons using 2 LL CLP 944: Lecture 3 1 The Two Sides
More informationOutline. Mixed models in R using the lme4 package Part 3: Longitudinal data. Sleep deprivation data. Simple longitudinal data
Outline Mixed models in R using the lme4 package Part 3: Longitudinal data Douglas Bates Longitudinal data: sleepstudy A model with random effects for intercept and slope University of Wisconsin - Madison
More informationOutline. Statistical inference for linear mixed models. One-way ANOVA in matrix-vector form
Outline Statistical inference for linear mixed models Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark general form of linear mixed models examples of analyses using linear mixed
More informationExploring Hierarchical Linear Mixed Models
Exploring Hierarchical Linear Mixed Models 1/49 Last time... A Greenhouse Experiment testing C:N Ratios Sam was testing how changing the C:N Ratio of soil affected plant leaf growth. He had 3 treatments.
More informationIntroduction and Background to Multilevel Analysis
Introduction and Background to Multilevel Analysis Dr. J. Kyle Roberts Southern Methodist University Simmons School of Education and Human Development Department of Teaching and Learning Background and
More informationMixed effects models
Mixed effects models The basic theory and application in R Mitchel van Loon Research Paper Business Analytics Mixed effects models The basic theory and application in R Author: Mitchel van Loon Research
More informationCorrelated Data: Linear Mixed Models with Random Intercepts
1 Correlated Data: Linear Mixed Models with Random Intercepts Mixed Effects Models This lecture introduces linear mixed effects models. Linear mixed models are a type of regression model, which generalise
More informationMutilevel Models: Pooled and Clustered Data
Mutilevel Models: Pooled and Clustered Data ICPSR, June 1-5, 2015 Tom Carsey carsey@unc.edu In the Beginning... We start with a simple linear model Where i subscripts individual observations. Y i = β 0
More informationA (Brief) Introduction to Crossed Random Effects Models for Repeated Measures Data
A (Brief) Introduction to Crossed Random Effects Models for Repeated Measures Data Today s Class: Review of concepts in multivariate data Introduction to random intercepts Crossed random effects models
More informationMixed models in R using the lme4 package Part 2: Longitudinal data, modeling interactions
Mixed models in R using the lme4 package Part 2: Longitudinal data, modeling interactions Douglas Bates Department of Statistics University of Wisconsin - Madison Madison January 11, 2011
More informationHierarchical Linear Models (HLM) Using R Package nlme. Interpretation. 2 = ( x 2) u 0j. e ij
Hierarchical Linear Models (HLM) Using R Package nlme Interpretation I. The Null Model Level 1 (student level) model is mathach ij = β 0j + e ij Level 2 (school level) model is β 0j = γ 00 + u 0j Combined
More informationHierarchical Generalized Linear Models. ERSH 8990 REMS Seminar on HLM Last Lecture!
Hierarchical Generalized Linear Models ERSH 8990 REMS Seminar on HLM Last Lecture! Hierarchical Generalized Linear Models Introduction to generalized models Models for binary outcomes Interpreting parameter
More informationIntroduction to the Analysis of Hierarchical and Longitudinal Data
Introduction to the Analysis of Hierarchical and Longitudinal Data Georges Monette, York University with Ye Sun SPIDA June 7, 2004 1 Graphical overview of selected concepts Nature of hierarchical models
More informationLongitudinal Data Analysis
Longitudinal Data Analysis Mike Allerhand This document has been produced for the CCACE short course: Longitudinal Data Analysis. No part of this document may be reproduced, in any form or by any means,
More informationIntroduction to Random Effects of Time and Model Estimation
Introduction to Random Effects of Time and Model Estimation Today s Class: The Big Picture Multilevel model notation Fixed vs. random effects of time Random intercept vs. random slope models How MLM =
More informationValue Added Modeling
Value Added Modeling Dr. J. Kyle Roberts Southern Methodist University Simmons School of Education and Human Development Department of Teaching and Learning Background for VAMs Recall from previous lectures
More informationA Re-Introduction to General Linear Models (GLM)
A Re-Introduction to General Linear Models (GLM) Today s Class: You do know the GLM Estimation (where the numbers in the output come from): From least squares to restricted maximum likelihood (REML) Reviewing
More informationStatistical Distribution Assumptions of General Linear Models
Statistical Distribution Assumptions of General Linear Models Applied Multilevel Models for Cross Sectional Data Lecture 4 ICPSR Summer Workshop University of Colorado Boulder Lecture 4: Statistical Distributions
More information10. Alternative case influence statistics
10. Alternative case influence statistics a. Alternative to D i : dffits i (and others) b. Alternative to studres i : externally-studentized residual c. Suggestion: use whatever is convenient with the
More informationReview of CLDP 944: Multilevel Models for Longitudinal Data
Review of CLDP 944: Multilevel Models for Longitudinal Data Topics: Review of general MLM concepts and terminology Model comparisons and significance testing Fixed and random effects of time Significance
More informationBIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation
BIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation Yujin Chung November 29th, 2016 Fall 2016 Yujin Chung Lec13: MLE Fall 2016 1/24 Previous Parametric tests Mean comparisons (normality assumption)
More informationFinal Review. Yang Feng. Yang Feng (Columbia University) Final Review 1 / 58
Final Review Yang Feng http://www.stat.columbia.edu/~yangfeng Yang Feng (Columbia University) Final Review 1 / 58 Outline 1 Multiple Linear Regression (Estimation, Inference) 2 Special Topics for Multiple
More informationRegression: Main Ideas Setting: Quantitative outcome with a quantitative explanatory variable. Example, cont.
TCELL 9/4/205 36-309/749 Experimental Design for Behavioral and Social Sciences Simple Regression Example Male black wheatear birds carry stones to the nest as a form of sexual display. Soler et al. wanted
More informationlme4 Luke Chang Last Revised July 16, Fitting Linear Mixed Models with a Varying Intercept
lme4 Luke Chang Last Revised July 16, 2010 1 Using lme4 1.1 Fitting Linear Mixed Models with a Varying Intercept We will now work through the same Ultimatum Game example from the regression section and
More information20. REML Estimation of Variance Components. Copyright c 2018 (Iowa State University) 20. Statistics / 36
20. REML Estimation of Variance Components Copyright c 2018 (Iowa State University) 20. Statistics 510 1 / 36 Consider the General Linear Model y = Xβ + ɛ, where ɛ N(0, Σ) and Σ is an n n positive definite
More information36-463/663Multilevel and Hierarchical Models
36-463/663Multilevel and Hierarchical Models From Bayes to MCMC to MLMs Brian Junker 132E Baker Hall brian@stat.cmu.edu 1 Outline Bayesian Statistics and MCMC Distribution of Skill Mastery in a Population
More informationCoping with Additional Sources of Variation: ANCOVA and Random Effects
Coping with Additional Sources of Variation: ANCOVA and Random Effects 1/49 More Noise in Experiments & Observations Your fixed coefficients are not always so fixed Continuous variation between samples
More informationRandom and Mixed Effects Models - Part II
Random and Mixed Effects Models - Part II Statistics 149 Spring 2006 Copyright 2006 by Mark E. Irwin Two-Factor Random Effects Model Example: Miles per Gallon (Neter, Kutner, Nachtsheim, & Wasserman, problem
More information36-309/749 Experimental Design for Behavioral and Social Sciences. Sep. 22, 2015 Lecture 4: Linear Regression
36-309/749 Experimental Design for Behavioral and Social Sciences Sep. 22, 2015 Lecture 4: Linear Regression TCELL Simple Regression Example Male black wheatear birds carry stones to the nest as a form
More informationover Time line for the means). Specifically, & covariances) just a fixed variance instead. PROC MIXED: to 1000 is default) list models with TYPE=VC */
CLP 944 Example 4 page 1 Within-Personn Fluctuation in Symptom Severity over Time These data come from a study of weekly fluctuation in psoriasis severity. There was no intervention and no real reason
More informationStat/F&W Ecol/Hort 572 Review Points Ané, Spring 2010
1 Linear models Y = Xβ + ɛ with ɛ N (0, σ 2 e) or Y N (Xβ, σ 2 e) where the model matrix X contains the information on predictors and β includes all coefficients (intercept, slope(s) etc.). 1. Number of
More informationWorkshop 9.3a: Randomized block designs
-1- Workshop 93a: Randomized block designs Murray Logan November 23, 16 Table of contents 1 Randomized Block (RCB) designs 1 2 Worked Examples 12 1 Randomized Block (RCB) designs 11 RCB design Simple Randomized
More informationMixed models in R using the lme4 package Part 7: Generalized linear mixed models
Mixed models in R using the lme4 package Part 7: Generalized linear mixed models Douglas Bates University of Wisconsin - Madison and R Development Core Team University of
More informationOutline. Mixed models in R using the lme4 package Part 5: Generalized linear mixed models. Parts of LMMs carried over to GLMMs
Outline Mixed models in R using the lme4 package Part 5: Generalized linear mixed models Douglas Bates University of Wisconsin - Madison and R Development Core Team UseR!2009,
More informationMixed models in R using the lme4 package Part 5: Generalized linear mixed models
Mixed models in R using the lme4 package Part 5: Generalized linear mixed models Douglas Bates Madison January 11, 2011 Contents 1 Definition 1 2 Links 2 3 Example 7 4 Model building 9 5 Conclusions 14
More informationStat 579: Generalized Linear Models and Extensions
Stat 579: Generalized Linear Models and Extensions Linear Mixed Models for Longitudinal Data Yan Lu April, 2018, week 14 1 / 64 Data structure and Model t1 t2 tn i 1st subject y 11 y 12 y 1n1 2nd subject
More informationOutline. Example and Model ANOVA table F tests Pairwise treatment comparisons with LSD Sample and subsample size determination
Outline 1 The traditional approach 2 The Mean Squares approach for the Completely randomized design (CRD) CRD and one-way ANOVA Variance components and the F test Inference about the intercept Sample vs.
More informationRandom and Mixed Effects Models - Part III
Random and Mixed Effects Models - Part III Statistics 149 Spring 2006 Copyright 2006 by Mark E. Irwin Quasi-F Tests When we get to more than two categorical factors, some times there are not nice F tests
More information36-720: Linear Mixed Models
36-720: Linear Mixed Models Brian Junker October 8, 2007 Review: Linear Mixed Models (LMM s) Bayesian Analogues Facilities in R Computational Notes Predictors and Residuals Examples [Related to Christensen
More information36-463/663: Multilevel & Hierarchical Models
36-463/663: Multilevel & Hierarchical Models Some Random Effects Configurations Brian Junker 132E Baker Hall brian@stat.cmu.edu 1 Outline Random Effect Configurations Most of our models so far: Level 1
More informationA Re-Introduction to General Linear Models
A Re-Introduction to General Linear Models Today s Class: Big picture overview Why we are using restricted maximum likelihood within MIXED instead of least squares within GLM Linear model interpretation
More information36-309/749 Experimental Design for Behavioral and Social Sciences. Dec 1, 2015 Lecture 11: Mixed Models (HLMs)
36-309/749 Experimental Design for Behavioral and Social Sciences Dec 1, 2015 Lecture 11: Mixed Models (HLMs) Independent Errors Assumption An error is the deviation of an individual observed outcome (DV)
More informationMultilevel Models in Matrix Form. Lecture 7 July 27, 2011 Advanced Multivariate Statistical Methods ICPSR Summer Session #2
Multilevel Models in Matrix Form Lecture 7 July 27, 2011 Advanced Multivariate Statistical Methods ICPSR Summer Session #2 Today s Lecture Linear models from a matrix perspective An example of how to do
More informationLecture 9 STK3100/4100
Lecture 9 STK3100/4100 27. October 2014 Plan for lecture: 1. Linear mixed models cont. Models accounting for time dependencies (Ch. 6.1) 2. Generalized linear mixed models (GLMM, Ch. 13.1-13.3) Examples
More informationStat 5303 (Oehlert): Randomized Complete Blocks 1
Stat 5303 (Oehlert): Randomized Complete Blocks 1 > library(stat5303libs);library(cfcdae);library(lme4) > immer Loc Var Y1 Y2 1 UF M 81.0 80.7 2 UF S 105.4 82.3 3 UF V 119.7 80.4 4 UF T 109.7 87.2 5 UF
More informationMixed models in R using the lme4 package Part 5: Generalized linear mixed models
Mixed models in R using the lme4 package Part 5: Generalized linear mixed models Douglas Bates 2011-03-16 Contents 1 Generalized Linear Mixed Models Generalized Linear Mixed Models When using linear mixed
More informationIntroducing Generalized Linear Models: Logistic Regression
Ron Heck, Summer 2012 Seminars 1 Multilevel Regression Models and Their Applications Seminar Introducing Generalized Linear Models: Logistic Regression The generalized linear model (GLM) represents and
More informationUsing R formulae to test for main effects in the presence of higher-order interactions
Using R formulae to test for main effects in the presence of higher-order interactions Roger Levy arxiv:1405.2094v2 [stat.me] 15 Jan 2018 January 16, 2018 Abstract Traditional analysis of variance (ANOVA)
More informationOutline for today. Two-way analysis of variance with random effects
Outline for today Two-way analysis of variance with random effects Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark Two-way ANOVA using orthogonal projections March 4, 2018 1 /
More informationRon Heck, Fall Week 8: Introducing Generalized Linear Models: Logistic Regression 1 (Replaces prior revision dated October 20, 2011)
Ron Heck, Fall 2011 1 EDEP 768E: Seminar in Multilevel Modeling rev. January 3, 2012 (see footnote) Week 8: Introducing Generalized Linear Models: Logistic Regression 1 (Replaces prior revision dated October
More informationA Handbook of Statistical Analyses Using R 2nd Edition. Brian S. Everitt and Torsten Hothorn
A Handbook of Statistical Analyses Using R 2nd Edition Brian S. Everitt and Torsten Hothorn CHAPTER 12 Analysing Longitudinal Data I: Computerised Delivery of Cognitive Behavioural Therapy Beat the Blues
More informationOutline for today. Maximum likelihood estimation. Computation with multivariate normal distributions. Multivariate normal distribution
Outline for today Maximum likelihood estimation Rasmus Waageetersen Deartment of Mathematics Aalborg University Denmark October 30, 2007 the multivariate normal distribution linear and linear mixed models
More informationDescribing Within-Person Change over Time
Describing Within-Person Change over Time Topics: Multilevel modeling notation and terminology Fixed and random effects of linear time Predicted variances and covariances from random slopes Dependency
More informationA Handbook of Statistical Analyses Using R 2nd Edition. Brian S. Everitt and Torsten Hothorn
A Handbook of Statistical Analyses Using R 2nd Edition Brian S. Everitt and Torsten Hothorn CHAPTER 12 Analysing Longitudinal Data I: Computerised Delivery of Cognitive Behavioural Therapy Beat the Blues
More informationSimple Linear Regression for the MPG Data
Simple Linear Regression for the MPG Data 2000 2500 3000 3500 15 20 25 30 35 40 45 Wgt MPG What do we do with the data? y i = MPG of i th car x i = Weight of i th car i =1,...,n n = Sample Size Exploratory
More information36-463/663: Multilevel & Hierarchical Models
36-463/663: Multilevel & Hierarchical Models Causal Inference Brian Junker 132E Baker Hall brian@stat.cmu.edu 1 Outline Causal Inference [G&H Ch 9] The Fundamental Problem Confounders, and how Controlled
More informationWeakly informative priors
Department of Statistics and Department of Political Science Columbia University 23 Apr 2014 Collaborators (in order of appearance): Gary King, Frederic Bois, Aleks Jakulin, Vince Dorie, Sophia Rabe-Hesketh,
More informationGeneralized Linear Models for Non-Normal Data
Generalized Linear Models for Non-Normal Data Today s Class: 3 parts of a generalized model Models for binary outcomes Complications for generalized multivariate or multilevel models SPLH 861: Lecture
More informationTopic 17 - Single Factor Analysis of Variance. Outline. One-way ANOVA. The Data / Notation. One way ANOVA Cell means model Factor effects model
Topic 17 - Single Factor Analysis of Variance - Fall 2013 One way ANOVA Cell means model Factor effects model Outline Topic 17 2 One-way ANOVA Response variable Y is continuous Explanatory variable is
More informationAn Introduction to Path Analysis
An Introduction to Path Analysis PRE 905: Multivariate Analysis Lecture 10: April 15, 2014 PRE 905: Lecture 10 Path Analysis Today s Lecture Path analysis starting with multivariate regression then arriving
More informationEPSY 905: Fundamentals of Multivariate Modeling Online Lecture #7
Introduction to Generalized Univariate Models: Models for Binary Outcomes EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #7 EPSY 905: Intro to Generalized In This Lecture A short review
More information13. October p. 1
Lecture 8 STK3100/4100 Linear mixed models 13. October 2014 Plan for lecture: 1. The lme function in the nlme library 2. Induced correlation structure 3. Marginal models 4. Estimation - ML and REML 5.
More informationAdditional Notes: Investigating a Random Slope. When we have fixed level-1 predictors at level 2 we show them like this:
Ron Heck, Summer 01 Seminars 1 Multilevel Regression Models and Their Applications Seminar Additional Notes: Investigating a Random Slope We can begin with Model 3 and add a Random slope parameter. If
More informationThe Simple Linear Regression Model
The Simple Linear Regression Model Lesson 3 Ryan Safner 1 1 Department of Economics Hood College ECON 480 - Econometrics Fall 2017 Ryan Safner (Hood College) ECON 480 - Lesson 3 Fall 2017 1 / 77 Bivariate
More informationIntroduction to Statistical modeling: handout for Math 489/583
Introduction to Statistical modeling: handout for Math 489/583 Statistical modeling occurs when we are trying to model some data using statistical tools. From the start, we recognize that no model is perfect
More information36-720: The Rasch Model
36-720: The Rasch Model Brian Junker October 15, 2007 Multivariate Binary Response Data Rasch Model Rasch Marginal Likelihood as a GLMM Rasch Marginal Likelihood as a Log-Linear Model Example For more
More informationWeakly informative priors
Department of Statistics and Department of Political Science Columbia University 21 Oct 2011 Collaborators (in order of appearance): Gary King, Frederic Bois, Aleks Jakulin, Vince Dorie, Sophia Rabe-Hesketh,
More informationTopic 20: Single Factor Analysis of Variance
Topic 20: Single Factor Analysis of Variance Outline Single factor Analysis of Variance One set of treatments Cell means model Factor effects model Link to linear regression using indicator explanatory
More informationCovariance Models (*) X i : (n i p) design matrix for fixed effects β : (p 1) regression coefficient for fixed effects
Covariance Models (*) Mixed Models Laird & Ware (1982) Y i = X i β + Z i b i + e i Y i : (n i 1) response vector X i : (n i p) design matrix for fixed effects β : (p 1) regression coefficient for fixed
More informationStatistics 262: Intermediate Biostatistics Model selection
Statistics 262: Intermediate Biostatistics Model selection Jonathan Taylor & Kristin Cobb Statistics 262: Intermediate Biostatistics p.1/?? Today s class Model selection. Strategies for model selection.
More informationStatistical Inference: The Marginal Model
Statistical Inference: The Marginal Model Edps/Psych/Stat 587 Carolyn J. Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Fall 2017 Outline Inference for fixed
More informationLinear Regression Models P8111
Linear Regression Models P8111 Lecture 25 Jeff Goldsmith April 26, 2016 1 of 37 Today s Lecture Logistic regression / GLMs Model framework Interpretation Estimation 2 of 37 Linear regression Course started
More informationA Brief and Friendly Introduction to Mixed-Effects Models in Linguistics
A Brief and Friendly Introduction to Mixed-Effects Models in Linguistics Cluster-specific parameters ( random effects ) Σb Parameters governing inter-cluster variability b1 b2 bm x11 x1n1 x21 x2n2 xm1
More informationMixed Model Theory, Part I
enote 4 1 enote 4 Mixed Model Theory, Part I enote 4 INDHOLD 2 Indhold 4 Mixed Model Theory, Part I 1 4.1 Design matrix for a systematic linear model.................. 2 4.2 The mixed model.................................
More informationSimple, Marginal, and Interaction Effects in General Linear Models
Simple, Marginal, and Interaction Effects in General Linear Models PRE 905: Multivariate Analysis Lecture 3 Today s Class Centering and Coding Predictors Interpreting Parameters in the Model for the Means
More informationMIXED MODELS FOR REPEATED (LONGITUDINAL) DATA PART 2 DAVID C. HOWELL 4/1/2010
MIXED MODELS FOR REPEATED (LONGITUDINAL) DATA PART 2 DAVID C. HOWELL 4/1/2010 Part 1 of this document can be found at http://www.uvm.edu/~dhowell/methods/supplements/mixed Models for Repeated Measures1.pdf
More informationIntroduction to Bayesian Statistics and Markov Chain Monte Carlo Estimation. EPSY 905: Multivariate Analysis Spring 2016 Lecture #10: April 6, 2016
Introduction to Bayesian Statistics and Markov Chain Monte Carlo Estimation EPSY 905: Multivariate Analysis Spring 2016 Lecture #10: April 6, 2016 EPSY 905: Intro to Bayesian and MCMC Today s Class An
More informationLinear Models (continued)
Linear Models (continued) Model Selection Introduction Most of our previous discussion has focused on the case where we have a data set and only one fitted model. Up until this point, we have discussed,
More informationAnalysis of binary repeated measures data with R
Analysis of binary repeated measures data with R Right-handed basketball players take right and left-handed shots from 3 locations in a different random order for each player. Hit or miss is recorded.
More informationStat 209 Lab: Linear Mixed Models in R This lab covers the Linear Mixed Models tutorial by John Fox. Lab prepared by Karen Kapur. ɛ i Normal(0, σ 2 )
Lab 2 STAT209 1/31/13 A complication in doing all this is that the package nlme (lme) is supplanted by the new and improved lme4 (lmer); both are widely used so I try to do both tracks in separate Rogosa
More informationSTAT 350. Assignment 4
STAT 350 Assignment 4 1. For the Mileage data in assignment 3 conduct a residual analysis and report your findings. I used the full model for this since my answers to assignment 3 suggested we needed the
More informationReview of Multiple Regression
Ronald H. Heck 1 Let s begin with a little review of multiple regression this week. Linear models [e.g., correlation, t-tests, analysis of variance (ANOVA), multiple regression, path analysis, multivariate
More informationModel Estimation Example
Ronald H. Heck 1 EDEP 606: Multivariate Methods (S2013) April 7, 2013 Model Estimation Example As we have moved through the course this semester, we have encountered the concept of model estimation. Discussions
More informationIntroduction to SAS proc mixed
Faculty of Health Sciences Introduction to SAS proc mixed Analysis of repeated measurements, 2017 Julie Forman Department of Biostatistics, University of Copenhagen Outline Data in wide and long format
More informationTime Invariant Predictors in Longitudinal Models
Time Invariant Predictors in Longitudinal Models Longitudinal Data Analysis Workshop Section 9 University of Georgia: Institute for Interdisciplinary Research in Education and Human Development Section
More informationClass Notes: Week 8. Probit versus Logit Link Functions and Count Data
Ronald Heck Class Notes: Week 8 1 Class Notes: Week 8 Probit versus Logit Link Functions and Count Data This week we ll take up a couple of issues. The first is working with a probit link function. While
More informationEstimation: Problems & Solutions
Estimation: Problems & Solutions Edps/Psych/Stat 587 Carolyn J. Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Fall 2017 Outline 1. Introduction: Estimation of
More informationA Handbook of Statistical Analyses Using R. Brian S. Everitt and Torsten Hothorn
A Handbook of Statistical Analyses Using R Brian S. Everitt and Torsten Hothorn CHAPTER 10 Analysing Longitudinal Data I: Computerised Delivery of Cognitive Behavioural Therapy Beat the Blues 10.1 Introduction
More informationRon Heck, Fall Week 3: Notes Building a Two-Level Model
Ron Heck, Fall 2011 1 EDEP 768E: Seminar on Multilevel Modeling rev. 9/6/2011@11:27pm Week 3: Notes Building a Two-Level Model We will build a model to explain student math achievement using student-level
More informationModel comparison. Patrick Breheny. March 28. Introduction Measures of predictive power Model selection
Model comparison Patrick Breheny March 28 Patrick Breheny BST 760: Advanced Regression 1/25 Wells in Bangladesh In this lecture and the next, we will consider a data set involving modeling the decisions
More informationOverdispersion Workshop in generalized linear models Uppsala, June 11-12, Outline. Overdispersion
Biostokastikum Overdispersion is not uncommon in practice. In fact, some would maintain that overdispersion is the norm in practice and nominal dispersion the exception McCullagh and Nelder (1989) Overdispersion
More informationUNIVERSITY OF TORONTO. Faculty of Arts and Science APRIL 2010 EXAMINATIONS STA 303 H1S / STA 1002 HS. Duration - 3 hours. Aids Allowed: Calculator
UNIVERSITY OF TORONTO Faculty of Arts and Science APRIL 2010 EXAMINATIONS STA 303 H1S / STA 1002 HS Duration - 3 hours Aids Allowed: Calculator LAST NAME: FIRST NAME: STUDENT NUMBER: There are 27 pages
More informationSimple, Marginal, and Interaction Effects in General Linear Models: Part 1
Simple, Marginal, and Interaction Effects in General Linear Models: Part 1 PSYC 943 (930): Fundamentals of Multivariate Modeling Lecture 2: August 24, 2012 PSYC 943: Lecture 2 Today s Class Centering and
More informationWU Weiterbildung. Linear Mixed Models
Linear Mixed Effects Models WU Weiterbildung SLIDE 1 Outline 1 Estimation: ML vs. REML 2 Special Models On Two Levels Mixed ANOVA Or Random ANOVA Random Intercept Model Random Coefficients Model Intercept-and-Slopes-as-Outcomes
More informationIntroduction to SAS proc mixed
Faculty of Health Sciences Introduction to SAS proc mixed Analysis of repeated measurements, 2017 Julie Forman Department of Biostatistics, University of Copenhagen 2 / 28 Preparing data for analysis The
More informationLecture 2: Linear Models. Bruce Walsh lecture notes Seattle SISG -Mixed Model Course version 23 June 2011
Lecture 2: Linear Models Bruce Walsh lecture notes Seattle SISG -Mixed Model Course version 23 June 2011 1 Quick Review of the Major Points The general linear model can be written as y = X! + e y = vector
More informationModel Selection in GLMs. (should be able to implement frequentist GLM analyses!) Today: standard frequentist methods for model selection
Model Selection in GLMs Last class: estimability/identifiability, analysis of deviance, standard errors & confidence intervals (should be able to implement frequentist GLM analyses!) Today: standard frequentist
More informationGeneralized Linear. Mixed Models. Methods and Applications. Modern Concepts, Walter W. Stroup. Texts in Statistical Science.
Texts in Statistical Science Generalized Linear Mixed Models Modern Concepts, Methods and Applications Walter W. Stroup CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint
More information