Random Coefficients Model Examples
|
|
- Shana Short
- 6 years ago
- Views:
Transcription
1 Random Coefficients Model Examples STAT:5201 Week 15 - Lecture 2 1 / 26
2 Each subject (or experimental unit) has multiple measurements (this could be over time, or it could be multiple measurements on a continuous variable x). Random effects are included in the model as to allow for a random intercept and a random slope, essentially allowing for a separate line for each subject. with b i = ( b0i b 1i Y ij = β 0 + β 1 x ij + b 0i + b 1i x ij + ɛ ) N ( [ 0 0, d11 d 12 d 21 d 22 b and ɛ independent from each other. Individual: E[Y ij b 0i, b 1i ] = (β 0 + b 0i ) + (β 1 + b 1i )x ij ]) and ɛ ij N(0, σ 2 ) Marginal: E[Y ij ] = β 0 + β 1 x ij 2 / 26
3 Adapted from John Fox example in Linear Mixed Models Appendix found in An R and S-PLUS Companion to Applied Regression. Eating disorders can be difficult to treat as many patients do not feel the need for treatment, when friends & family recognize the severity. Even after patients with eating disorders are hospitalized, they can continue behaviors that can be detrimental to their health. Here, data were collected recording the amount of exercise of 138 teenage girls hospitalized for eating disorders, and on a group of 93 control subjects. Variables subject: a factor with subject id codes. age: age in years. exercise: hours per week of exercise. group: factor indicating patient or control.. 3 / 26
4 We will consider the example in R and in SAS. > library(car) > data(blackmore) > head(blackmore) subject age exercise group patient patient patient patient patient patient > dim(blackmore) [1] > length(unique(blackmore$subject[blackmore$group=="patient"])) [1] 138 > length(unique(blackmore$subject[blackmore$group=="control"])) [1] 93 4 / 26
5 Fox transformed the response variable for numerous reasons (described in text) as log 2 (y + 5/60). > Blackmore$log.exercise <- log(blackmore$exercise + 5/60, 2) > attach(blackmore) Investigating the data with plots (in R). Use a random sample of 20 girls from each group for trend plotting. The groupeddata object from the nlme package is used to form the trellis plots. > library(nlme) > chosen.pat.ids=sample(unique(subject[group=="patient"]), 20) > chosen.pat.20=groupeddata(log.exercise ~ age subject, data=blackmore[is.element(subject,chosen.pat.ids),]) > chosen.con.ids=sample(unique(subject[group=="control"]), 20) > chosen.con.20=groupeddata(log.exercise ~ age subject, data=blackmore[is.element(subject,chosen.con.ids),]) 5 / 26
6 > print(plot(chosen.con.20, main="control Subjects", xlab="age",ylab="log2 Exercise", ylim=1.2*range(chosen.con.20$log.exercise, chosen.pat.20$log.exercise), layout=c(5,4),aspect=1), position=c(0, 0, 0.5, 1), more=t) > print(plot(chosen.pat.20, main="patients", xlab="age",ylab="log2 Exercise", ylim=1.2*range(chosen.con.20$log.exercise, chosen.pat.20$log.exercise), layout=c(5,4),aspect=1), position=c(0.5, 0, 1, 1)) 6 / 26
7 The groupeddata object is automatically plotted in order by average exercise. The subjects with the highest exercise values are in the top row, the subjects with the lowest exercise values are in the bottom row. Control Subjects Patients log2 Exercise log2 Exercise Age Age 7 / 26
8 Investigating the data with plots (in SAS). After I created subsetted data sets of patients from each group called control1 and patient1 (8 subjects per group), I used the PROC SGPANEL procedure to plot the individual trajectories. Here I ve asked for a linear regression line for each subject, but you can simply connect the observed points using the vline option instead of the reg option. proc sgpanel data=control1; title Control Subjects ; panelby subject/columns=4 rows=2; reg x=age y=log_exercise; rowaxis min=-4 max=4; colaxis values=(8, 10, 12, 14, 16); run; proc sgpanel data=patient1; <similar coding for the patient group as control group> 8 / 26
9 9 / 26
10 10 / 26
11 You can also plot the overlay of these individual lines using PROC SGPLOT... proc sgplot data=control1; title Subset of Control Subjects ; reg x=age y=log_exercise/group=subject; run; 11 / 26
12 Investigating the subject-specific parameter estimates (in R). Fox formally fits a linear regression to each subject (231 separately fit models) in order to investigate the variability and correlation in the slopes and intercept estimates from a graphical perspective. The predictor age is transformed to represent age after the start of the study or age-8. He points out that the random coefficients model (fitted to all the data) fits a unified model that considers slopes and intercepts as random effects, and in that case, the estimated random effects or û are estimated using BLUPs (best linear unbiased predictors). 12 / 26
13 Investigating the subject-specific parameter estimates (in R). For a model with independent random subject effects (i.e. just a random intercept, as in the gene expression line example from earlier), the BLUPs are actually shrinkage estimators and fall between the individual observed values and the overall mean values. Formally, BLUPs are estimated as û = ĜZ ˆΣ 1 (y X ˆβ) where Σ = var(y) = ZGZ + R. 13 / 26
14 Before moving to a unified mixed model, we consider truly fitting a separate line to each subject (so, not a random coefficients model). Again, the nlme package is utilized when employing the lmlist function: > pat.list=lmlist(log.exercise ~ I(age - 8) subject, subset = group=="patient", data=blackmore) > con.list=lmlist(log.exercise ~ I(age - 8) subject, subset = group=="control", data=blackmore) > pat.coef=coef(pat.list) > con.coef=coef(con.list) > par(mfrow=c(1,2)) > boxplot(pat.coef[,1], con.coef[,1], main="intercepts", names=c("patients","controls")) > boxplot(pat.coef[,2], con.coef[,2], main="slope", names=c("patients","controls")) 14 / 26
15 Intercepts Slope Patients Controls Patients Controls The intercept represents the level of exercise at the start of the study. As expected, there is a great deal of variation in both the intercepts and the slopes. The median intercepts are fairly similar for patients and controls, but there is somewhat more variation among patients. The slopes are higher on average for patients than for controls and the slopes tend to be positive (suggesting their exercise increases over time). 15 / 26
16 It makes sense to also plot the relationship between the estimated intercept and slope parameters. The dataellipse function is in the car library. > plot(c(-5,4),c(-1.2,1.2),xlab="intercept",ylab="slope",type="n", main="(individual) Estimates of slope and intercept") > points(con.coef[,1],con.coef[,2],col=1) > points(pat.coef[,1],pat.coef[,2],col=2) > abline(v=0) > abline(h=0) > legend(-4.5,-.7,c("controls","patients"),col=c(1,2),pch=c(1,1)) > dataellipse(con.coef[,1],con.coef[,2],levels=c(.5,.95),add=true, plot.points=false,col=1) > dataellipse(pat.coef[,1],pat.coef[,2],levels=c(.5,.95),add=true, plot.points=false,col=2) 16 / 26
17 (Individual) Estimates of slope and intercept slope Controls Patients intercept Recall that we are on the log-scale base 2 for our response, so y = 0 coincides with 1 hour of exercise a week. It looks like the two groups have a reasonably similar correlation structure for the slope and intercept. It also looks like the patients have a shifted distribution such that they tend to have higher slopes. 17 / 26
18 Fitting the random coefficients model in SAS This model allows for a random slope and random intercept for each subject (which are allowed to be correlated). The population-level mean structure allows for separate lines for each treatment group (control and patient). The predictor age is transformed to represent age after the start of the study or age-8. data Blackmore; set Blackmore; age_trans = age-8; run; proc mixed data=blackmore; class subject group; model log_exercise = group age_trans group*age_trans/solution ddfm=satterth; random intercept age_trans/subject=subject type=un gcorr; run; 18 / 26
19 The Mixed Procedure Dimensions Covariance Parameters 4 Columns in X 6 Columns in Z Per Subject 2 Subjects 231 Max Obs Per Subject 5 Estimated G Correlation Matrix Row Effect subject Col1 Col2 1 Intercept age_trans We see that the correlation between b 0i and b 1i is estimated to be negative (ρ = ). 19 / 26
20 The Mixed Procedure Covariance Parameter Estimates Standard Z Cov Parm Subject Estimate Error Value Pr Z UN(1,1) subject <.0001 UN(2,1) subject UN(2,2) subject Residual <.0001 We see that the correlation between b 0i and b 1i is estimated to be negative (ρ = ) and marginally significant with a p= / 26
21 Solution for Fixed Effects Standard Effect group Estimate Error DF t Value Pr > t Intercept <.0001 group control group patient age_trans <.0001 age_trans*group control <.0001 age_trans*group patient Type 3 Tests of Fixed Effects Num Den Effect DF DF F Value Pr > F group age_trans <.0001 age_trans*group <.0001 The groups do not have significantly different intercepts (average exercise values at start of study, at age 8), but they do have significantly different slopes with the patient group having a higher slope than the control group. 21 / 26
22 I can capture the estimated BLUPs or û = ĜZ ˆΣ 1 (y X ˆβ) where Σ = var(y) = ZGZ + R using the ODS output and the solution option in the random statement: ods output SolutionR=blups; proc mixed data=blackmore covtest; class subject group; model log_exercise = group age_trans group*age_trans/ddfm=satterth; random intercept age_trans/subject=subject type=un gcorr solution; run; /* Solution for the random effects are BLUPs*/ ods output close; 22 / 26
23 proc print data=blups (obs=10); run; StdErr Obs Effect subject Estimate Pred DF tvalue Probt 1 Intercept age_trans Intercept age_trans Intercept age_trans Intercept age_trans Intercept age_trans / 26
24 Below I ve plotted the estimated BLUPs for the random slopes against the estimated slopes from the separately fit regression lines (in absolute values). Absolute values of slopes BLUP of slope separately fit slope (individual regression) 24 / 26
25 Fitting the random coefficients model (in R) Using the lme function in the nlme package, we see the same estimates for the covariance parameters as in SAS: > lme.1=lme(log.exercise~i(age-8)*group, random=~i(age-8) subject, data=blackmore) > summary(lme.1) Linear mixed-effects model fit by REML Data: Blackmore Random effects: Formula: ~I(age - 8) subject Structure: General positive-definite, Log-Cholesky parametrization StdDev Corr (Intercept) (Intr) I(age - 8) Residual / 26
26 Fitting the random coefficients model (in R) Square the estimates to match the SAS estimates: Var(Intercept)= = Var(slope)= = Corr(Intercept, slope)= = Var(Residual)= = / 26
SAS Syntax and Output for Data Manipulation:
CLP 944 Example 5 page 1 Practice with Fixed and Random Effects of Time in Modeling Within-Person Change The models for this example come from Hoffman (2015) chapter 5. We will be examining the extent
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 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 informationRandom Coefficient Model (a.k.a. multilevel model) (Adapted from UCLA Statistical Computing Seminars)
STAT:5201 Applied Statistic II Random Coefficient Model (a.k.a. multilevel model) (Adapted from UCLA Statistical Computing Seminars) School math achievement scores The data file consists of 7185 students
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 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 informationAnswer to exercise: Blood pressure lowering drugs
Answer to exercise: Blood pressure lowering drugs The data set bloodpressure.txt contains data from a cross-over trial, involving three different formulations of a drug for lowering of blood pressure:
More informationBest Linear Unbiased Prediction (BLUP) of Random Effects in the Normal Linear Mixed Effects Model. *Modified notes from Dr. Dan Nettleton from ISU
Best Linear Unbiased Prediction (BLUP) of Random Effects in the Normal Linear Mixed Effects Model *Modified notes from Dr. Dan Nettleton from ISU Suppose intelligence quotients (IQs) for a population of
More informationBiostatistics Workshop Longitudinal Data Analysis. Session 4 GARRETT FITZMAURICE
Biostatistics Workshop 2008 Longitudinal Data Analysis Session 4 GARRETT FITZMAURICE Harvard University 1 LINEAR MIXED EFFECTS MODELS Motivating Example: Influence of Menarche on Changes in Body Fat Prospective
More informationSTAT 5200 Handout #23. Repeated Measures Example (Ch. 16)
Motivating Example: Glucose STAT 500 Handout #3 Repeated Measures Example (Ch. 16) An experiment is conducted to evaluate the effects of three diets on the serum glucose levels of human subjects. Twelve
More informationStatistical Analysis of Hierarchical Data. David Zucker Hebrew University, Jerusalem, Israel
Statistical Analysis of Hierarchical Data David Zucker Hebrew University, Jerusalem, Israel Unit 1 Linear Mixed Models 1 Examples of Hierarchical Data Hierarchical data = subunits within units Students
More informationMixed-Effects Models in R
Mixed-Effects Models in R An Appendix to An R Companion to Applied Regression, Second Edition John Fox & Sanford Weisberg last revision (DRAFT): November 6, 2014 Abstract 1 Introduction The normal linear
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 informationRandom Effects. Edps/Psych/Stat 587. Carolyn J. Anderson. Fall Department of Educational Psychology. university of illinois at urbana-champaign
Random Effects Edps/Psych/Stat 587 Carolyn J. Anderson Department of Educational Psychology I L L I N O I S university of illinois at urbana-champaign Fall 2012 Outline Introduction Empirical Bayes inference
More informationSAS Code for Data Manipulation: SPSS Code for Data Manipulation: STATA Code for Data Manipulation: Psyc 945 Example 1 page 1
Psyc 945 Example page Example : Unconditional Models for Change in Number Match 3 Response Time (complete data, syntax, and output available for SAS, SPSS, and STATA electronically) These data come from
More informationIES 612/STA 4-573/STA Winter 2008 Week 1--IES 612-STA STA doc
IES 612/STA 4-573/STA 4-576 Winter 2008 Week 1--IES 612-STA 4-573-STA 4-576.doc Review Notes: [OL] = Ott & Longnecker Statistical Methods and Data Analysis, 5 th edition. [Handouts based on notes prepared
More informationSubject-specific observed profiles of log(fev1) vs age First 50 subjects in Six Cities Study
Subject-specific observed profiles of log(fev1) vs age First 50 subjects in Six Cities Study 1.4 0.0-6 7 8 9 10 11 12 13 14 15 16 17 18 19 age Model 1: A simple broken stick model with knot at 14 fit with
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 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 informationLab 11. Multilevel Models. Description of Data
Lab 11 Multilevel Models Henian Chen, M.D., Ph.D. Description of Data MULTILEVEL.TXT is clustered data for 386 women distributed across 40 groups. ID: 386 women, id from 1 to 386, individual level (level
More informationMixed models with correlated measurement errors
Mixed models with correlated measurement errors Rasmus Waagepetersen October 9, 2018 Example from Department of Health Technology 25 subjects where exposed to electric pulses of 11 different durations
More informationSAS Syntax and Output for Data Manipulation: CLDP 944 Example 3a page 1
CLDP 944 Example 3a page 1 From Between-Person to Within-Person Models for Longitudinal Data The models for this example come from Hoffman (2015) chapter 3 example 3a. We will be examining the extent to
More informationST505/S697R: Fall Homework 2 Solution.
ST505/S69R: Fall 2012. Homework 2 Solution. 1. 1a; problem 1.22 Below is the summary information (edited) from the regression (using R output); code at end of solution as is code and output for SAS. a)
More informationRandomized Block Designs with Replicates
LMM 021 Randomized Block ANOVA with Replicates 1 ORIGIN := 0 Randomized Block Designs with Replicates prepared by Wm Stein Randomized Block Designs with Replicates extends the use of one or more random
More informationTesting Indirect Effects for Lower Level Mediation Models in SAS PROC MIXED
Testing Indirect Effects for Lower Level Mediation Models in SAS PROC MIXED Here we provide syntax for fitting the lower-level mediation model using the MIXED procedure in SAS as well as a sas macro, IndTest.sas
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 informationLinear Mixed Models. Appendix to An R and S-PLUS Companion to Applied Regression. John Fox. May 2002
Linear Mixed Models Appendix to An R and S-PLUS Companion to Applied Regression John Fox May 22 1 Introduction The normal linear model (described, for example, in Chapter 4 of the text), y i = β 1 x 1i
More informationRandom Effects. Edps/Psych/Stat 587. Carolyn J. Anderson. Fall Department of Educational Psychology. c Board of Trustees, University of Illinois
Random Effects Edps/Psych/Stat 587 Carolyn J. Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Fall 2017 Outline Introduction Empirical Bayes inference Henderson
More information22s:152 Applied Linear Regression. Returning to a continuous response variable Y...
22s:152 Applied Linear Regression Generalized Least Squares Returning to a continuous response variable Y... Ordinary Least Squares Estimation The classical models we have fit so far with a continuous
More informationSTAT 3A03 Applied Regression With SAS Fall 2017
STAT 3A03 Applied Regression With SAS Fall 2017 Assignment 2 Solution Set Q. 1 I will add subscripts relating to the question part to the parameters and their estimates as well as the errors and residuals.
More informationRandom Intercept Models
Random Intercept Models Edps/Psych/Soc 589 Carolyn J. Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Spring 2019 Outline A very simple case of a random intercept
More information22s:152 Applied Linear Regression. In matrix notation, we can write this model: Generalized Least Squares. Y = Xβ + ɛ with ɛ N n (0, Σ)
22s:152 Applied Linear Regression Generalized Least Squares Returning to a continuous response variable Y Ordinary Least Squares Estimation The classical models we have fit so far with a continuous response
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 informationDe-mystifying random effects models
De-mystifying random effects models Peter J Diggle Lecture 4, Leahurst, October 2012 Linear regression input variable x factor, covariate, explanatory variable,... output variable y response, end-point,
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 15 1 / 38 Data structure t1 t2 tn i 1st subject y 11 y 12 y 1n1 Experimental 2nd subject
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 informationFaculty of Health Sciences. Correlated data. Variance component models. Lene Theil Skovgaard & Julie Lyng Forman.
Faculty of Health Sciences Correlated data Variance component models Lene Theil Skovgaard & Julie Lyng Forman November 27, 2018 1 / 84 Overview One-way anova with random variation The rabbit example Hierarchical
More informationCorrelated data. Overview. Example: Swelling due to vaccine. Variance component models. Faculty of Health Sciences. Variance component models
Faculty of Health Sciences Overview Correlated data Variance component models One-way anova with random variation The rabbit example Hierarchical models with several levels Random regression Lene Theil
More informationLecture 1 Linear Regression with One Predictor Variable.p2
Lecture Linear Regression with One Predictor Variablep - Basics - Meaning of regression parameters p - β - the slope of the regression line -it indicates the change in mean of the probability distn of
More informationPractice with Interactions among Continuous Predictors in General Linear Models (as estimated using restricted maximum likelihood in SAS MIXED)
CLP 944 Example 2b page 1 Practice with Interactions among Continuous Predictors in General Linear Models (as estimated using restricted maximum likelihood in SAS MIXED) The models for this example come
More informationNon-independence due to Time Correlation (Chapter 14)
Non-independence due to Time Correlation (Chapter 14) When we model the mean structure with ordinary least squares, the mean structure explains the general trends in the data with respect to our dependent
More informationIntroduction to Mixed Models in R
Introduction to Mixed Models in R Galin Jones School of Statistics University of Minnesota http://www.stat.umn.edu/ galin March 2011 Second in a Series Sponsored by Quantitative Methods Collaborative.
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 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 informationTopic 25 - One-Way Random Effects Models. Outline. Random Effects vs Fixed Effects. Data for One-way Random Effects Model. One-way Random effects
Topic 5 - One-Way Random Effects Models One-way Random effects Outline Model Variance component estimation - Fall 013 Confidence intervals Topic 5 Random Effects vs Fixed Effects Consider factor with numerous
More informationCOMPREHENSIVE WRITTEN EXAMINATION, PAPER III FRIDAY AUGUST 26, 2005, 9:00 A.M. 1:00 P.M. STATISTICS 174 QUESTION
COMPREHENSIVE WRITTEN EXAMINATION, PAPER III FRIDAY AUGUST 26, 2005, 9:00 A.M. 1:00 P.M. STATISTICS 174 QUESTION Answer all parts. Closed book, calculators allowed. It is important to show all working,
More informationLecture 4 Multiple linear regression
Lecture 4 Multiple linear regression BIOST 515 January 15, 2004 Outline 1 Motivation for the multiple regression model Multiple regression in matrix notation Least squares estimation of model parameters
More informationANOVA Longitudinal Models for the Practice Effects Data: via GLM
Psyc 943 Lecture 25 page 1 ANOVA Longitudinal Models for the Practice Effects Data: via GLM Model 1. Saturated Means Model for Session, E-only Variances Model (BP) Variances Model: NO correlation, EQUAL
More informationVariance components and LMMs
Faculty of Health Sciences Variance components and LMMs Analysis of repeated measurements, 4th December 2014 Julie Lyng Forman & Lene Theil Skovgaard Department of Biostatistics, University of Copenhagen
More informationUNIVERSITY OF MASSACHUSETTS Department of Mathematics and Statistics Applied Statistics Friday, January 15, 2016
UNIVERSITY OF MASSACHUSETTS Department of Mathematics and Statistics Applied Statistics Friday, January 15, 2016 Work all problems. 60 points are needed to pass at the Masters Level and 75 to pass at the
More informationStatistics for exp. medical researchers Regression and Correlation
Faculty of Health Sciences Regression analysis Statistics for exp. medical researchers Regression and Correlation Lene Theil Skovgaard Sept. 28, 2015 Linear regression, Estimation and Testing Confidence
More informationDifferences of Least Squares Means
STAT:5201 Homework 9 Solutions 1. We have a model with two crossed random factors operator and machine. There are 4 operators, 8 machines, and 3 observations from each operator/machine combination. (a)
More informationVariance components and LMMs
Faculty of Health Sciences Topics for today Variance components and LMMs Analysis of repeated measurements, 4th December 04 Leftover from 8/: Rest of random regression example. New concepts for today:
More informationFaculty of Health Sciences. Correlated data. Variance component models. Lene Theil Skovgaard & Julie Lyng Forman.
Faculty of Health Sciences Correlated data Variance component models Lene Theil Skovgaard & Julie Lyng Forman November 28, 2017 1 / 96 Overview One-way anova with random variation The rabbit example Hierarchical
More informationCorrelated data. Overview. Variance component models. Terminology for correlated measurements. Faculty of Health Sciences. Variance component models
Faculty of Health Sciences Overview Correlated data Variance component models Lene Theil Skovgaard & Julie Lyng Forman November 28, 2017 One-way anova with random variation The rabbit example Hierarchical
More informationRepeated Measures Modeling With PROC MIXED E. Barry Moser, Louisiana State University, Baton Rouge, LA
Paper 188-29 Repeated Measures Modeling With PROC MIXED E. Barry Moser, Louisiana State University, Baton Rouge, LA ABSTRACT PROC MIXED provides a very flexible environment in which to model many types
More information36-402/608 Homework #10 Solutions 4/1
36-402/608 Homework #10 Solutions 4/1 1. Fixing Breakout 17 (60 points) You must use SAS for this problem! Modify the code in wallaby.sas to load the wallaby data and to create a new outcome in the form
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 informationCorrelated data. Repeated measurements over time. Typical set-up for repeated measurements. Traditional presentation of data
Faculty of Health Sciences Repeated measurements over time Correlated data NFA, May 22, 2014 Longitudinal measurements Julie Lyng Forman & Lene Theil Skovgaard Department of Biostatistics University of
More information17. Example SAS Commands for Analysis of a Classic Split-Plot Experiment 17. 1
17 Example SAS Commands for Analysis of a Classic SplitPlot Experiment 17 1 DELIMITED options nocenter nonumber nodate ls80; Format SCREEN OUTPUT proc import datafile"c:\data\simulatedsplitplotdatatxt"
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 informationAdvantages of Mixed-effects Regression Models (MRM; aka multilevel, hierarchical linear, linear mixed models) 1. MRM explicitly models individual
Advantages of Mixed-effects Regression Models (MRM; aka multilevel, hierarchical linear, linear mixed models) 1. MRM explicitly models individual change across time 2. MRM more flexible in terms of repeated
More informationLecture 6 Multiple Linear Regression, cont.
Lecture 6 Multiple Linear Regression, cont. BIOST 515 January 22, 2004 BIOST 515, Lecture 6 Testing general linear hypotheses Suppose we are interested in testing linear combinations of the regression
More informationMeasuring relationships among multiple responses
Measuring relationships among multiple responses Linear association (correlation, relatedness, shared information) between pair-wise responses is an important property used in almost all multivariate analyses.
More informationMultiple Predictor Variables: ANOVA
Multiple Predictor Variables: ANOVA 1/32 Linear Models with Many Predictors Multiple regression has many predictors BUT - so did 1-way ANOVA if treatments had 2 levels What if there are multiple treatment
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 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 informationUniversity of Minnesota Duluth
Trend analyses for species of concern: Analysis of CPUE data for walleye, cisco, and smallmouth bass 1970-2008 Dave Staples 1, Lucinda Johnson 2, Jennifer Olker 2, Dan Brenneman 2 1 Minnesota Department
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 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 information14 Multiple Linear Regression
B.Sc./Cert./M.Sc. Qualif. - Statistics: Theory and Practice 14 Multiple Linear Regression 14.1 The multiple linear regression model In simple linear regression, the response variable y is expressed in
More informationProperties of the least squares estimates
Properties of the least squares estimates 2019-01-18 Warmup Let a and b be scalar constants, and X be a scalar random variable. Fill in the blanks E ax + b) = Var ax + b) = Goal Recall that the least squares
More informationCorrelated data. Longitudinal data. Typical set-up for repeated measurements. Examples from literature, I. Faculty of Health Sciences
Faculty of Health Sciences Longitudinal data Correlated data Longitudinal measurements Outline Designs Models for the mean Covariance patterns Lene Theil Skovgaard November 27, 2015 Random regression Baseline
More informationAnalysis of variance and regression. May 13, 2008
Analysis of variance and regression May 13, 2008 Repeated measurements over time Presentation of data Traditional ways of analysis Variance component model (the dogs revisited) Random regression Baseline
More information6. Multiple regression - PROC GLM
Use of SAS - November 2016 6. Multiple regression - PROC GLM Karl Bang Christensen Department of Biostatistics, University of Copenhagen. http://biostat.ku.dk/~kach/sas2016/ kach@biostat.ku.dk, tel: 35327491
More informationCovariance Structure Approach to Within-Cases
Covariance Structure Approach to Within-Cases Remember how the data file grapefruit1.data looks: Store sales1 sales2 sales3 1 62.1 61.3 60.8 2 58.2 57.9 55.1 3 51.6 49.2 46.2 4 53.7 51.5 48.3 5 61.4 58.7
More informationMixed-Models. version 30 October 2011
Mixed-Models version 30 October 2011 Mixed models Mixed models estimate a vector! of fixed effects and one (or more) vectors u of random effects Both fixed and random effects models always include a vector
More informationSimple Linear Regression
Simple Linear Regression ST 370 Regression models are used to study the relationship of a response variable and one or more predictors. The response is also called the dependent variable, and the predictors
More informationModels for longitudinal data
Faculty of Health Sciences Contents Models for longitudinal data Analysis of repeated measurements, NFA 016 Julie Lyng Forman & Lene Theil Skovgaard Department of Biostatistics, University of Copenhagen
More informationGeneral Linear Model (Chapter 4)
General Linear Model (Chapter 4) Outcome variable is considered continuous Simple linear regression Scatterplots OLS is BLUE under basic assumptions MSE estimates residual variance testing regression coefficients
More informationLecture 3: Inference in SLR
Lecture 3: Inference in SLR STAT 51 Spring 011 Background Reading KNNL:.1.6 3-1 Topic Overview This topic will cover: Review of hypothesis testing Inference about 1 Inference about 0 Confidence Intervals
More informationSTAT 705 Chapters 23 and 24: Two factors, unequal sample sizes; multi-factor ANOVA
STAT 705 Chapters 23 and 24: Two factors, unequal sample sizes; multi-factor ANOVA Timothy Hanson Department of Statistics, University of South Carolina Stat 705: Data Analysis II 1 / 22 Balanced vs. unbalanced
More informationMixed-Model Estimation of genetic variances. Bruce Walsh lecture notes Uppsala EQG 2012 course version 28 Jan 2012
Mixed-Model Estimation of genetic variances Bruce Walsh lecture notes Uppsala EQG 01 course version 8 Jan 01 Estimation of Var(A) and Breeding Values in General Pedigrees The above designs (ANOVA, P-O
More informationPLS205 Lab 2 January 15, Laboratory Topic 3
PLS205 Lab 2 January 15, 2015 Laboratory Topic 3 General format of ANOVA in SAS Testing the assumption of homogeneity of variances by "/hovtest" by ANOVA of squared residuals Proc Power for ANOVA One-way
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 informationTwo-Variable Regression Model: The Problem of Estimation
Two-Variable Regression Model: The Problem of Estimation Introducing the Ordinary Least Squares Estimator Jamie Monogan University of Georgia Intermediate Political Methodology Jamie Monogan (UGA) Two-Variable
More informationLeast Squares Estimation
Least Squares Estimation Using the least squares estimator for β we can obtain predicted values and compute residuals: Ŷ = Z ˆβ = Z(Z Z) 1 Z Y ˆɛ = Y Ŷ = Y Z(Z Z) 1 Z Y = [I Z(Z Z) 1 Z ]Y. The usual decomposition
More informationVarians- og regressionsanalyse
Faculty of Health Sciences Varians- og regressionsanalyse Variance component models Lene Theil Skovgaard Department of Biostatistics Variance component models Definitions and motivation One-way anova with
More informationSerial Correlation. Edps/Psych/Stat 587. Carolyn J. Anderson. Fall Department of Educational Psychology
Serial Correlation Edps/Psych/Stat 587 Carolyn J. Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Fall 017 Model for Level 1 Residuals There are three sources
More informationTechnische Universität München. Zentrum Mathematik. Linear Mixed Models Applied to Bank Branch Deposit Data
Technische Universität München Zentrum Mathematik Linear Mixed Models Applied to Bank Branch Deposit Data Project by Eike Christian Brechmann Supervisor: Prof. Claudia Czado, Ph.D. Tutor: Dr. Mathias Hofmann
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 informationApplied Regression Analysis
Applied Regression Analysis Chapter 3 Multiple Linear Regression Hongcheng Li April, 6, 2013 Recall simple linear regression 1 Recall simple linear regression 2 Parameter Estimation 3 Interpretations of
More informationCorrelated data. Variance component models. Example: Evaluate vaccine. Traditional assumption so far. Faculty of Health Sciences
Faculty of Health Sciences Variance component models Definitions and motivation Correlated data Variance component models, I Lene Theil Skovgaard November 29, 2013 One-way anova with random variation The
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 informationSTAT 705 Generalized linear mixed models
STAT 705 Generalized linear mixed models Timothy Hanson Department of Statistics, University of South Carolina Stat 705: Data Analysis II 1 / 24 Generalized Linear Mixed Models We have considered random
More informationAnalysis of variance and regression. December 4, 2007
Analysis of variance and regression December 4, 2007 Variance component models Variance components One-way anova with random variation estimation interpretations Two-way anova with random variation Crossed
More informationSTATISTICS 479 Exam II (100 points)
Name STATISTICS 79 Exam II (1 points) 1. A SAS data set was created using the following input statement: Answer parts(a) to (e) below. input State $ City $ Pop199 Income Housing Electric; (a) () Give the
More informationSTAT3401: Advanced data analysis Week 10: Models for Clustered Longitudinal Data
STAT3401: Advanced data analysis Week 10: Models for Clustered Longitudinal Data Berwin Turlach School of Mathematics and Statistics Berwin.Turlach@gmail.com The University of Western Australia Models
More information#Alternatively we could fit a model where the rail values are levels of a factor with fixed effects
examples-lme.r Tue Nov 25 12:32:20 2008 1 library(nlme) # The following data shows the results of tests carried over 6 rails. The response # indicated the time needed for a an ultrasonic wave to travel
More informationRegression Models - Introduction
Regression Models - Introduction In regression models there are two types of variables that are studied: A dependent variable, Y, also called response variable. It is modeled as random. An independent
More information