Analysis of Longitudinal Data: Comparison Between PROC GLM and PROC MIXED. Maribeth Johnson Medical College of Georgia Augusta, GA
|
|
- Helen Townsend
- 6 years ago
- Views:
Transcription
1 Analysis of Longitudinal Data: Comparison Between PROC GLM and PROC MIXED Maribeth Johnson Medical College of Georgia Augusta, GA
2 Overview Introduction to longitudinal data Describe the data for examples (WT and HR) Analysis using PROC GLM and the RANDOM statement Analysis using PROC GLM and the REPEATED statement Analysis using PROC MIXED Compare error V-C estimates from each Use data with no missing values Repeat using data with missing values Conclusions and Implications
3 Longitudinal Data Multiple measurements of a response variable On the same experimental unit Made over a period of time Experimental units are considered as random effects
4 Correlated Data Quantify the inter-relationships of the measurements If variances and covariances are constant then the relationship can be thought of as ICC reliability Assumption of compound symmetry can be tested against other variance structures using PROC MIXED
5 Correlated Data The ability to model different covariance structures provides an opportunity to investigate and possibly quantify the tracking of k measurements as well as provide the basis for other relationships y ~ N( μ, Σ) where Σ 1 1 = M 1 κ 1 M κ L L O L 1κ κ M k
6 Data Texas site in the studies of child activity and nutrition (SCAN) program - TXSCAN Children aged 3 or 4 years at the start
7 Data Structure Number of children Number of records WT HR 1 of of of of
8 Data Structure WT HR Year N Mean SD N Mean SD
9 Using GLM With The Random Statement Univariate GLM proc glm; class child year; model depvar = child year; random child; Mixed model univariate analysis of variance Random statement prints table of expected MS
10 Using GLM With The Random Statement Univariate GLM The GLM Procedure Source Type III Expected Mean Square CHILD Var(Error) + 4 Var(CHILD) YEAR Var(Error) + Q(YEAR)
11 Using GLM With The Random Statement Univariate GLM E(MS) Child = Var(Error) + k Var(CHILD) k is the average number of observations per child CHILD is treated as a fixed effect Variance components calculated using Method of Moments Var(CHILD) = MS(CHILD) - MS(Error) k ICC = VAR(CHILD) VAR(CHILD) + VAR(Error)
12 Using GLM With The Random Statement Univariate GLM Reliability of a single year of measurement The individual components of variation are maximum-likelihood estimators when the data are balanced
13 Using GLM With The Repeated Statement MV GLM Can t handle missing data Repeated measurements must appear in a multivariate mode in the dataset, one observation for each experimental unit Allows for the multivariate test of the assumption of compound symmetry (CS) If the CS assumption is rejected it can t help in determining the correct underlying covariance structure
14 Using GLM With The Repeated Statement MV GLM proc glm; model depvar1-depvar4 = /nouni; repeated year / printe; NOUNI option suppresses the univariate analyses of each year PRINTE option outputs the partial correlations computed from residuals after fitting the between-subjects model
15 Test for Sphericity Tests whether a set of orthonormal contrasts of the repeated measures variables are independent and have equal variances, i.e. are the data compound symmetric Significance tells you that this condition is not met The estimate of the correlation between measures from the univariate GLM analysis is not valid
16 Test for Sphericity Correlations from the PRINTE option of this analysis are identical to the correlations computed from the RCORR option in the REPEATED statement in PROC MIXED when TYPE=UN is specified as the covariance structure REPEATED statements perform different functions
17 Using PROC MIXED Data is in the same format as GLM using the Random statement (i.e. Univariate GLM) Multiple observations for each experimental unit
18 Standard Linear Model y = Xβ + ε y X is vector of observed data is the design matrix β is unkown vector ο f fixed effects ε is unknown random error vector, iid N ~ (0, )
19 Mixed Model y= Xβ + Zγ + ε Z is a design matrix γ is unkown vector ο f random effects ε is unknown random error vector, no longer required to be iid
20 Using PROC MIXED proc mixed; class child year; model depvar = year; repeated year / subject=child r rcorr type=cov-structure ; REPEATED statement models the covariance structures in R, the variance-covariance matrix of the vector of errors If no REPEATED statement is specified, R is assumed to be equal to I.
21 Using PROC MIXED proc mixed; class child year; model depvar = year; repeated year / subject=child r rcorr type=cov-structure ; SUBJECT=CHILD option is the mechanism for block diagonalizing R
22 Using PROC MIXED proc mixed; class child year; model depvar = year; repeated year / subject=child r type=cov-structure ; rcorr R option of the REPEATED statement requests that the first block of the R matrix be printed
23 Using PROC MIXED proc mixed; class child year; model depvar = year; repeated year / subject=child r rcorr type=cov-structure ; RCORR options prints the correlation matrix corresponding to R
24 Using PROC MIXED proc mixed; class child year; model depvar = year; repeated year / subject=child r rcorr type=cov-structure ; TYPE= option is what determines the V-C structure
25 Covariance Structures Compound symmetric (CS): Most specific structure Variance within years is constant Common correlation between years Two parameters estimated Assumption of Univariate GLM estimates
26 Covariance Structures Heterogeneous compound symmetric (CSH): Common correlation Different variances along the diagonal Five parameters estimated for these data ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ
27 Covariance Structures First-order autoregressive (AR(1)): Variance within years is constant Estimate of the autoregressive parameter Correlations between years separated by the same amount of time are the same (ρ m ) Two parameters estimated ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ
28 Covariance Structures Heterogeneous first-order autoregressive (ARH(1)): Different variances along the diagonal Estimate of the autoregressive parameter Correlations between years separated by the same amount of time are the same (ρ m ) Five parameters estimated for these data ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ ρ
29 Covariance Structures Unstructured (UN): Estimates of all four variances and six covariances for these data All of the correlations between years may be different Identical to those from the PRINTE option of the REPEATED statement in GLM
30 Determine The Preferred Model Likelihood Ratio Test (LRT) One model is a submodel of another Compute - times the difference between their residual log likelihoods (-RLL) Chi-square distribution with degrees of freedom equal to the difference in the number of parameters for the two models Models are preferred where the -RLL is smaller
31 Determine The Preferred Model Akaike s Information Criterion (AIC) Schwarz s Bayesian Criterion (BIC) Model that has the smallest value is the preferred model BIC penalizes models with more covariance parameters more than AIC
32 Determine The Preferred Model Fit Statistics - Res Log Likelihood AIC (smaller is better) AICC (smaller is better) BIC (smaller is better)
33 Example WT - Balanced Data WT Year N Mean SD
34 V-C and Correlations WT - Balanced Data (N=98) Type = CS: (Variance and covariances in top line, correlations below) Type = CSH: (Variances on diagonal, covariances above, correlations below) parameters estimated from the CS structure are identical to those calculated using the Univariate GLM mean square estimates
35 V-C and Correlations WT - Balanced Data (N=98) Type = CS: (Variance and covariances in top line, correlations below) Type = CSH: (Variances on diagonal, covariances above, correlations below) The sphericity test from the MV GLM was significant
36 V-C and Correlations WT - Balanced Data (N=98) Type = CS: (Variance and covariances in top line, correlations below) Type = CSH: (Variances on diagonal, covariances above, correlations below) The correlation from the CSH structure is identical to that calculated using the Univariate GLM analysis of WT standardized within year
37 Model Comparisons WT - Balanced Data Type Parameters -RLL Comparison Model Chi-square / df CS CSH CS 88 / 3 * *p<.0005 Heterogeneous variances along the diagonal provide a significantly better fit for the CS model
38 V-C and Correlations WT - Balanced Data (N=98) Type = AR(1): (Variance and covariances in top line, correlations below) Type = ARH(1): (Variances on diagonal, covariances above, correlations below) The correlations do not fall off as quickly when variances are allowed to differ
39 Model Comparisons WT - Balanced Data Type Parameters -RLL Comparison Model Chi-square / df AR(1) ARH(1) 5 19 AR(1) 6 / 3 * *p<.0005 Heterogeneous autoregressive covariance structure (ARH(1)) provides a better fit than the assumption that within year variances are equal (AR(1))
40 V-C and Correlations WT - Balanced Data (N=98) Type = CSH: Type = ARH(1): Type = UN: UN produced very similar results to ARH(1) UN identical to the correlations from the Multivariate GLM analysis
41 Model Comparisons WT - Balanced Data Type Parameters -RLL Comparison Model Chi-square / df CSH ARH(1) UN CSH 86 / 5 * UN ARH(1) 8 / 5 *p<.0005 UN provides a significantly better fit to the data than CSH ARH(1) appears to provide the best fit for these data Improper models may underestimate the correlation between adjacent measurements
42 Example WT - Unbalanced WT Year N Mean SD
43 V-C and Correlations WT - Unbalanced Data (N=57) Type = CS: (Variance and covariances in top line, correlations below) Type = CSH: (Variances on diagonal, covariances above, correlations below) Estimates from Univariate GLM no longer ML estimators, not identical to CS but similar Var: 1. vs 11.5 Cov: 10.1 vs 9.5
44 V-C and Correlations WT - Unbalanced Data (N=57) Type = CSH: Type = ARH(1): Type = UN: UN produced very similar results to ARH(1)
45 V-C and Correlations WT Preferred Models Type = ARH(1): Balanced data Type = ARH(1): Unbalanced data Correlation estimates are similar Anywhere from 37 to 145 more children measured in any one year
46 Example HR - Balanced Data HR Year N Mean SD
47 V-C and Correlations HR Balanced Data (N=60) ANOVA estimates and those using the CS structure are identical (r=0.55) Common correlation estimate is similar using CSH (r=0.56) Heterogenous variances do not provide a better fit for these data The MV GLM s test of sphericity was significant so CS is not the correct structure
48 V-C and Correlations HR - Balanced Data (N=60) Type = UN: (Variances on diagonal, covariances above, correlations below) The correlation between successive years is lower when the children are young
49 Model Comparisons WT - Balanced Data Type Parameters -RLL Comparison Model Chi-square / df CS AR(1) UN CS 17 / 8 * UN AR(1) 1 / 8 * *p<.05 UN provides a significantly better fit to the data than CS UN provides a significantly better fit to the data than AR(1)
50 V-C and Correlations HR Preferred Models Type = UN: Balanced data Type = UN: Unbalanced data Correlation estimates are similar Variance smaller in third year There are from 50 to 189 more children in unbalanced analysis
51 V-C and Correlations HR Preferred Models Type = UN: Balanced data Type = UN: Unbalanced data Correlation higher between year 3 and 4 As children age there may be higher correlation between adjacent years
52 Conclusions Univariate GLM V-C calculations are identical to MIXED estimates using type=cs when there is no missing data (ICC) MV GLM V-C estimates are identical to MIXED estimates using type=un If sphericity assumption violated, MV GLM cannot determine best fit
53 Conclusions The ability to model a broader class of V-C structures yield results that make more sense in the context of the problem
54 Conclusions WT variability increases with age Correlation higher between adjacent years ARH(1) appears useful in explaining how WT tracks in very young children
55 Conclusions The UN nature of the relationship between HR measurements aids in assessment of the quality of this measurement
56 Conclusions PROC MIXED gives an opportunity to understand and quantify the true inter-relationships between repeated measurements
57 Implications Inferences for fixed effects may be impacted by the poor choice of error covariance matrix Those impacted will be the repeated variable and its interactions
58 Implications - Example Group Effect of TM intervention on resting SBP N Obs Variable Time N Mean Std Dev TM 30 SUPSBP_A SUPSBP_C SUPSBP_D Control 3 SUPSBP_A SUPSBP_C SUPSBP_D
59 Implications - Example Effect of TM intervention on resting SBP TM or Control=TM TM or Control=Control sbp 160 sbp time time
60 Implications - Example Univariate GLM with Random statement proc glm data=long; class group time family; model sbp=group family(group) time group*time; random family(group) / test; lsmeans group*time/stderr; run; quit;
61 Implications - Example Univariate GLM with Random statement Source DF Type III SS Mean Square F Value Pr > F GROUP*time Error: MS(Error)
62 Implications - Example Multivariate GLM with Repeated statement proc glm data=wide; class group; model supsbp_a supsbp_c supsbp_d = group /nouni; repeated time 3 (0 48) / printe; lsmeans group/ stderr; run; quit;
63 Implications - Example Multivariate GLM with Repeated statement Univariate Tests of Hypotheses for Within Subject Effects Adj Pr > F Source DF Type III SS Mean Square F Value Pr > F G-G H-F time*group Error(time)
64 Implications - Example Multivariate GLM with Repeated statement Multivariate Tests - identical results for Pillai s Trace, Hotelling-Lawley Trace, and Roy s Greatest Root MANOVA Test Criteria and Exact F Statistics for the Hypothesis of no time*group Effect H = Type III SSCP Matrix for time*group E = Error SSCP Matrix S=1 M=0 N=8.5 Statistic Value F Value Num DF Den DF Pr > F Wilks' Lambda
65 Implications - Example Multivariate GLM with Repeated statement Partial Correlation Coefficients from the Error SSCP Matrix / Prob > r DF = 60 SUPSBP_A SUPSBP_C SUPSBP_D SUPSBP_A < <.0001 SUPSBP_C < <.0001 SUPSBP_D < < Sphericity Tests Variables DF Mauchly's Criterion Chi-Square Pr > ChiSq Transformed Variates <.0001 Orthogonal Components
66 Implications - Example PROC MIXED with Repeated statement proc mixed data=long; class group time family; model sbp= group time group*time/ddfm=kenwardroger; repeated time / subject= family(group) type=un r rcorr; lsmeans group*time; run; quit;
67 Implications - Example PROC MIXED with Repeated statement Type=UN and DDFM=kenwardroger Type 3 Tests of Fixed Effects Effect Num DF Den DF F Value Pr > F GROUP*time Estimated R Correlation Matrix for FAMILY(GROUP) 1003 Control Row Col1 Col Col
68 Implications - Example PROC MIXED with Repeated statement Type=CS Type 3 Tests of Fixed Effects Effect Num DF Den DF F Value Pr > F GROUP*time
69 Implications - Example PROC MIXED with Repeated statement Least Squares Means type=un Effect Group time Estimate Standard Error GROUP*time Control GROUP*time Control GROUP*time Control GROUP*time TM GROUP*time TM Least Squares Means - type=cs Estimate Standard Error GROUP*time TM
70 Implications - Example PROC MIXED with Repeated statement Diff -RLL = 9. with 4 df, p= Both AIC smaller for UN, BIC is not Fit Statistics UN 6 parameters estimated - Res Log Likelihood Fit Statistics CS parameters estimated - Res Log Likelihood AIC (smaller is better) AIC (smaller is better) AICC (smaller is better) AICC (smaller is better) BIC (smaller is better) BIC (smaller is better)
71 Implications Inferences for fixed effects may be impacted by the poor choice of error covariance matrix Those impacted will be the repeated variable and its interactions It is suggested that if Number of subjects - Number of treatment groups - Number of repeated measures > 30 then UN will usually suffice
72 Implications MV GLM results are identical to PROC MIXED results when ddfm=kenwardrogers is used
Analysis of Longitudinal Data: Comparison between PROC GLM and PROC MIXED.
Analysis of Longitudinal Data: Comparison between PROC GLM and PROC MIXED. Maribeth Johnson, Medical College of Georgia, Augusta, GA ABSTRACT Longitudinal data refers to datasets with multiple measurements
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 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 informationRepeated Measures Part 2: Cartoon data
Repeated Measures Part 2: Cartoon data /*********************** cartoonglm.sas ******************/ options linesize=79 noovp formdlim='_'; title 'Cartoon Data: STA442/1008 F 2005'; proc format; /* value
More informationWITHIN-PARTICIPANT EXPERIMENTAL DESIGNS
1 WITHIN-PARTICIPANT EXPERIMENTAL DESIGNS I. Single-factor designs: the model is: yij i j ij ij where: yij score for person j under treatment level i (i = 1,..., I; j = 1,..., n) overall mean βi treatment
More informationChapter 9. Multivariate and Within-cases Analysis. 9.1 Multivariate Analysis of Variance
Chapter 9 Multivariate and Within-cases Analysis 9.1 Multivariate Analysis of Variance Multivariate means more than one response variable at once. Why do it? Primarily because if you do parallel analyses
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 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 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 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 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 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 informationRepeated Measures ANOVA Multivariate ANOVA and Their Relationship to Linear Mixed Models
Repeated Measures ANOVA Multivariate ANOVA and Their Relationship to Linear Mixed Models EPSY 905: Multivariate Analysis Spring 2016 Lecture #12 April 20, 2016 EPSY 905: RM ANOVA, MANOVA, and Mixed Models
More informationSAS 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 information4.1 Computing section Example: Bivariate measurements on plants Post hoc analysis... 7
Master of Applied Statistics ST116: Chemometrics and Multivariate Statistical data Analysis Per Bruun Brockhoff Module 4: Computing 4.1 Computing section.................................. 1 4.1.1 Example:
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 informationOther hypotheses of interest (cont d)
Other hypotheses of interest (cont d) In addition to the simple null hypothesis of no treatment effects, we might wish to test other hypothesis of the general form (examples follow): H 0 : C k g β g p
More informationChapter 5: Multivariate Analysis and Repeated Measures
Chapter 5: Multivariate Analysis and Repeated Measures Multivariate -- More than one dependent variable at once. Why do it? Primarily because if you do parallel analyses on lots of outcome measures, the
More informationRepeated Measures Data
Repeated Measures Data Mixed Models Lecture Notes By Dr. Hanford page 1 Data where subjects are measured repeatedly over time - predetermined intervals (weekly) - uncontrolled variable intervals between
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 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 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 informationInternational Journal of Current Research in Biosciences and Plant Biology ISSN: Volume 2 Number 5 (May-2015) pp
Original Research Article International Journal of Current Research in Biosciences and Plant Biology ISSN: 349-00 Volume Number (May-01) pp. -19 www.ijcrbp.com Graphical Approaches to Support Mixed Model
More informationUV Absorbance by Fish Slime
Data Set 1: UV Absorbance by Fish Slime Statistical Setting This handout describes a repeated-measures ANOVA, with two crossed amongsubjects factors and repeated-measures on a third (within-subjects) factor.
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 informationThis is a Randomized Block Design (RBD) with a single factor treatment arrangement (2 levels) which are fixed.
EXST3201 Chapter 13c Geaghan Fall 2005: Page 1 Linear Models Y ij = µ + βi + τ j + βτij + εijk This is a Randomized Block Design (RBD) with a single factor treatment arrangement (2 levels) which are fixed.
More informationAnalysis of Repeated Measures Data of Iraqi Awassi Lambs Using Mixed Model
American Journal of Applied Scientific Research 01; 1(): 1-6 Published online November 1, 01 (http://www.sciencepublishinggroup.com/j/ajasr) doi: 10.64/j.ajasr.00.13 Analysis of Repeated Measures Data
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 informationApplication of Ghosh, Grizzle and Sen s Nonparametric Methods in. Longitudinal Studies Using SAS PROC GLM
Application of Ghosh, Grizzle and Sen s Nonparametric Methods in Longitudinal Studies Using SAS PROC GLM Chan Zeng and Gary O. Zerbe Department of Preventive Medicine and Biometrics University of Colorado
More informationBiostatistics 301A. Repeated measurement analysis (mixed models)
B a s i c S t a t i s t i c s F o r D o c t o r s Singapore Med J 2004 Vol 45(10) : 456 CME Article Biostatistics 301A. Repeated measurement analysis (mixed models) Y H Chan Faculty of Medicine National
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 informationGLM Repeated Measures
GLM Repeated Measures Notation The GLM (general linear model) procedure provides analysis of variance when the same measurement or measurements are made several times on each subject or case (repeated
More informationT. Mark Beasley One-Way Repeated Measures ANOVA handout
T. Mark Beasley One-Way Repeated Measures ANOVA handout Profile Analysis Example In the One-Way Repeated Measures ANOVA, two factors represent separate sources of variance. Their interaction presents an
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 informationANOVA approaches to Repeated Measures. repeated measures MANOVA (chapter 3)
ANOVA approaches to Repeated Measures univariate repeated-measures ANOVA (chapter 2) repeated measures MANOVA (chapter 3) Assumptions Interval measurement and normally distributed errors (homogeneous across
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 informationLecture 4. Random Effects in Completely Randomized Design
Lecture 4. Random Effects in Completely Randomized Design Montgomery: 3.9, 13.1 and 13.7 1 Lecture 4 Page 1 Random Effects vs Fixed Effects Consider factor with numerous possible levels Want to draw inference
More informationOdor attraction CRD Page 1
Odor attraction CRD Page 1 dm'log;clear;output;clear'; options ps=512 ls=99 nocenter nodate nonumber nolabel FORMCHAR=" ---- + ---+= -/\*"; ODS LISTING; *** Table 23.2 ********************************************;
More informationAnalyzing the Behavior of Rats by Repeated Measurements
Georgia State University ScholarWorks @ Georgia State University Mathematics Theses Department of Mathematics and Statistics 5-3-007 Analyzing the Behavior of Rats by Repeated Measurements Kenita A. Hall
More informationMultivariate analysis of variance and covariance
Introduction Multivariate analysis of variance and covariance Univariate ANOVA: have observations from several groups, numerical dependent variable. Ask whether dependent variable has same mean for each
More informationMIXED MODELS FOR REPEATED (LONGITUDINAL) DATA PART 1 DAVID C. HOWELL 4/26/2010
MIXED MODELS FOR REPEATED (LONGITUDINAL) DATA PART 1 DAVID C. HOWELL 4/26/2010 FOR THE SECOND PART OF THIS DOCUMENT GO TO www.uvm.edu/~dhowell/methods/supplements/mixed Models Repeated/Mixed Models for
More informationChapter 7, continued: MANOVA
Chapter 7, continued: MANOVA The Multivariate Analysis of Variance (MANOVA) technique extends Hotelling T 2 test that compares two mean vectors to the setting in which there are m 2 groups. We wish to
More informationdm'log;clear;output;clear'; options ps=512 ls=99 nocenter nodate nonumber nolabel FORMCHAR=" = -/\<>*"; ODS LISTING;
dm'log;clear;output;clear'; options ps=512 ls=99 nocenter nodate nonumber nolabel FORMCHAR=" ---- + ---+= -/\*"; ODS LISTING; *** Table 23.2 ********************************************; *** Moore, David
More informationPOWER ANALYSIS TO DETERMINE THE IMPORTANCE OF COVARIANCE STRUCTURE CHOICE IN MIXED MODEL REPEATED MEASURES ANOVA
POWER ANALYSIS TO DETERMINE THE IMPORTANCE OF COVARIANCE STRUCTURE CHOICE IN MIXED MODEL REPEATED MEASURES ANOVA A Thesis Submitted to the Graduate Faculty of the North Dakota State University of Agriculture
More informationDynamic Determination of Mixed Model Covariance Structures. in Double-blind Clinical Trials. Matthew Davis - Omnicare Clinical Research
PharmaSUG2010 - Paper SP12 Dynamic Determination of Mixed Model Covariance Structures in Double-blind Clinical Trials Matthew Davis - Omnicare Clinical Research Abstract With the computing power of SAS
More informationDescribing Within-Person Fluctuation over Time using Alternative Covariance Structures
Describing Within-Person Fluctuation over Time using Alternative Covariance Structures Today s Class: The Big Picture ACS models using the R matrix only Introducing the G, Z, and V matrices ACS models
More informationGeneral Linear Model. Notes Output Created Comments Input. 19-Dec :09:44
GET ILE='G:\lare\Data\Accuracy_Mixed.sav'. DATASET NAME DataSet WINDOW=RONT. GLM Jigsaw Decision BY CMCTools /WSACTOR= Polynomial /METHOD=SSTYPE(3) /PLOT=PROILE(CMCTools*) /EMMEANS=TABLES(CMCTools) COMPARE
More information1998, Gregory Carey Repeated Measures ANOVA - 1. REPEATED MEASURES ANOVA (incomplete)
1998, Gregory Carey Repeated Measures ANOVA - 1 REPEATED MEASURES ANOVA (incomplete) Repeated measures ANOVA (RM) is a specific type of MANOVA. When the within group covariance matrix has a special form,
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 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 informationTopic 12. The Split-plot Design and its Relatives (continued) Repeated Measures
12.1 Topic 12. The Split-plot Design and its Relatives (continued) Repeated Measures 12.9 Repeated measures analysis Sometimes researchers make multiple measurements on the same experimental unit. We have
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 informationPrepared by: Prof. Dr Bahaman Abu Samah Department of Professional Development and Continuing Education Faculty of Educational Studies Universiti
Prepared by: Prof. Dr Bahaman Abu Samah Department of Professional Development and Continuing Education Faculty of Educational Studies Universiti Putra Malaysia Serdang Use in experiment, quasi-experiment
More informationTopic 23: Diagnostics and Remedies
Topic 23: Diagnostics and Remedies Outline Diagnostics residual checks ANOVA remedial measures Diagnostics Overview We will take the diagnostics and remedial measures that we learned for regression and
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 informationApplied Multivariate Analysis
Department of Mathematics and Statistics, University of Vaasa, Finland Spring 2017 Discriminant Analysis Background 1 Discriminant analysis Background General Setup for the Discriminant Analysis Descriptive
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 informationStep 2: Select Analyze, Mixed Models, and Linear.
Example 1a. 20 employees were given a mood questionnaire on Monday, Wednesday and again on Friday. The data will be first be analyzed using a Covariance Pattern model. Step 1: Copy Example1.sav data file
More informationGeneral Linear Model
GLM V1 V2 V3 V4 V5 V11 V12 V13 V14 V15 /WSFACTOR=placeholders 2 Polynomial target 5 Polynomial /METHOD=SSTYPE(3) /EMMEANS=TABLES(OVERALL) /EMMEANS=TABLES(placeholders) COMPARE ADJ(SIDAK) /EMMEANS=TABLES(target)
More informationM M Cross-Over Designs
Chapter 568 Cross-Over Designs Introduction This module calculates the power for an x cross-over design in which each subject receives a sequence of treatments and is measured at periods (or time points).
More informationStevens 2. Aufl. S Multivariate Tests c
Stevens 2. Aufl. S. 200 General Linear Model Between-Subjects Factors 1,00 2,00 3,00 N 11 11 11 Effect a. Exact statistic Pillai's Trace Wilks' Lambda Hotelling's Trace Roy's Largest Root Pillai's Trace
More informationMultivariate Analysis of Variance
Chapter 15 Multivariate Analysis of Variance Jolicouer and Mosimann studied the relationship between the size and shape of painted turtles. The table below gives the length, width, and height (all in mm)
More informationMULTIVARIATE ANALYSIS OF VARIANCE
MULTIVARIATE ANALYSIS OF VARIANCE RAJENDER PARSAD AND L.M. BHAR Indian Agricultural Statistics Research Institute Library Avenue, New Delhi - 0 0 lmb@iasri.res.in. Introduction In many agricultural experiments,
More informationOutline. Topic 20 - Diagnostics and Remedies. Residuals. Overview. Diagnostics Plots Residual checks Formal Tests. STAT Fall 2013
Topic 20 - Diagnostics and Remedies - Fall 2013 Diagnostics Plots Residual checks Formal Tests Remedial Measures Outline Topic 20 2 General assumptions Overview Normally distributed error terms Independent
More informationMultivariate Linear Models
Multivariate Linear Models Stanley Sawyer Washington University November 7, 2001 1. Introduction. Suppose that we have n observations, each of which has d components. For example, we may have d measurements
More informationGraphical Procedures, SAS' PROC MIXED, and Tests of Repeated Measures Effects. H.J. Keselman University of Manitoba
1 Graphical Procedures, SAS' PROC MIXED, and Tests of Repeated Measures Effects by H.J. Keselman University of Manitoba James Algina University of Florida and Rhonda K. Kowalchuk University of Manitoba
More informationMultivariate Statistical Analysis
Multivariate Statistical Analysis Fall 2011 C. L. Williams, Ph.D. Lecture 17 for Applied Multivariate Analysis Outline Multivariate Analysis of Variance 1 Multivariate Analysis of Variance The hypotheses:
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 informationReview of Unconditional Multilevel Models for Longitudinal Data
Review of Unconditional Multilevel Models for Longitudinal Data Topics: Course (and MLM) overview Concepts in longitudinal multilevel modeling Model comparisons and significance testing Describing within-person
More informationAn Introduction to Multilevel Models. PSYC 943 (930): Fundamentals of Multivariate Modeling Lecture 25: December 7, 2012
An Introduction to Multilevel Models PSYC 943 (930): Fundamentals of Multivariate Modeling Lecture 25: December 7, 2012 Today s Class Concepts in Longitudinal Modeling Between-Person vs. +Within-Person
More informationunadjusted model for baseline cholesterol 22:31 Monday, April 19,
unadjusted model for baseline cholesterol 22:31 Monday, April 19, 2004 1 Class Level Information Class Levels Values TRETGRP 3 3 4 5 SEX 2 0 1 Number of observations 916 unadjusted model for baseline cholesterol
More informationMANOVA is an extension of the univariate ANOVA as it involves more than one Dependent Variable (DV). The following are assumptions for using MANOVA:
MULTIVARIATE ANALYSIS OF VARIANCE MANOVA is an extension of the univariate ANOVA as it involves more than one Dependent Variable (DV). The following are assumptions for using MANOVA: 1. Cell sizes : o
More informationExample 7b: Generalized Models for Ordinal Longitudinal Data using SAS GLIMMIX, STATA MEOLOGIT, and MPLUS (last proportional odds model only)
CLDP945 Example 7b page 1 Example 7b: Generalized Models for Ordinal Longitudinal Data using SAS GLIMMIX, STATA MEOLOGIT, and MPLUS (last proportional odds model only) This example comes from real data
More informationSTAT 501 EXAM I NAME Spring 1999
STAT 501 EXAM I NAME Spring 1999 Instructions: You may use only your calculator and the attached tables and formula sheet. You can detach the tables and formula sheet from the rest of this exam. Show your
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 informationLongitudinal data: simple univariate methods of analysis
Longitudinal data: simple univariate methods of analysis Danish version by Henrik Stryhn, June 1996 Department of Mathematcis and Physics, KVL Translation (and rewriting) by Ib Skovgaard, March 1998 (the
More informationRepeated-Measures ANOVA in SPSS Correct data formatting for a repeated-measures ANOVA in SPSS involves having a single line of data for each
Repeated-Measures ANOVA in SPSS Correct data formatting for a repeated-measures ANOVA in SPSS involves having a single line of data for each participant, with the repeated measures entered as separate
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 informationMixed Effects Models
Mixed Effects Models What is the effect of X on Y What is the effect of an independent variable on the dependent variable Independent variables are fixed factors. We want to measure their effect Random
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 informationRANDOM and REPEATED statements - How to Use Them to Model the Covariance Structure in Proc Mixed. Charlie Liu, Dachuang Cao, Peiqi Chen, Tony Zagar
Paper S02-2007 RANDOM and REPEATED statements - How to Use Them to Model the Covariance Structure in Proc Mixed Charlie Liu, Dachuang Cao, Peiqi Chen, Tony Zagar Eli Lilly & Company, Indianapolis, IN ABSTRACT
More informationLecture 10: Experiments with Random Effects
Lecture 10: Experiments with Random Effects Montgomery, Chapter 13 1 Lecture 10 Page 1 Example 1 A textile company weaves a fabric on a large number of looms. It would like the looms to be homogeneous
More informationApplied Multivariate Statistical Modeling Prof. J. Maiti Department of Industrial Engineering and Management Indian Institute of Technology, Kharagpur
Applied Multivariate Statistical Modeling Prof. J. Maiti Department of Industrial Engineering and Management Indian Institute of Technology, Kharagpur Lecture - 29 Multivariate Linear Regression- Model
More informationAccounting for Correlation in the Analysis of Randomized Controlled Trials with Multiple Layers of Clustering
Duquesne University Duquesne Scholarship Collection Electronic Theses and Dissertations Spring 2016 Accounting for Correlation in the Analysis of Randomized Controlled Trials with Multiple Layers of Clustering
More informationMultivariate Tests. Mauchly's Test of Sphericity
General Model Within-Sujects Factors Dependent Variale IDLS IDLF IDHS IDHF IDHCLS IDHCLF Descriptive Statistics IDLS IDLF IDHS IDHF IDHCLS IDHCLF Mean Std. Deviation N.0.70.0.0..8..88.8...97 Multivariate
More informationReview of Multilevel Models for Longitudinal Data
Review of Multilevel Models for Longitudinal Data Topics: Concepts in longitudinal multilevel modeling Describing within-person fluctuation using ACS models Describing within-person change using random
More informationApplied Multivariate and Longitudinal Data Analysis
Applied Multivariate and Longitudinal Data Analysis Chapter 2: Inference about the mean vector(s) II Ana-Maria Staicu SAS Hall 5220; 919-515-0644; astaicu@ncsu.edu 1 1 Compare Means from More Than Two
More informationDescribing Change over Time: Adding Linear Trends
Describing Change over Time: Adding Linear Trends Longitudinal Data Analysis Workshop Section 7 University of Georgia: Institute for Interdisciplinary Research in Education and Human Development Section
More informationGregory Carey, 1998 Regression & Path Analysis - 1 MULTIPLE REGRESSION AND PATH ANALYSIS
Gregory Carey, 1998 Regression & Path Analysis - 1 MULTIPLE REGRESSION AND PATH ANALYSIS Introduction Path analysis and multiple regression go hand in hand (almost). Also, it is easier to learn about multivariate
More information5.3 Three-Stage Nested Design Example
5.3 Three-Stage Nested Design Example A researcher designs an experiment to study the of a metal alloy. A three-stage nested design was conducted that included Two alloy chemistry compositions. Three ovens
More informationDescriptive Statistics
*following creates z scores for the ydacl statedp traitdp and rads vars. *specifically adding the /SAVE subcommand to descriptives will create z. *scores for whatever variables are in the command. DESCRIPTIVES
More informationLinear Mixed Models with Repeated Effects
1 Linear Mixed Models with Repeated Effects Introduction and Examples Using SAS/STAT Software Jerry W. Davis, University of Georgia, Griffin Campus. Introduction Repeated measures refer to measurements
More informationIncomplete Block Designs
Incomplete Block Designs Recall: in randomized complete block design, each of a treatments was used once within each of b blocks. In some situations, it will not be possible to use each of a treatments
More informationRepeated Measures Design. Advertising Sales Example
STAT:5201 Anaylsis/Applied Statistic II Repeated Measures Design Advertising Sales Example A company is interested in comparing the success of two different advertising campaigns. It has 10 test markets,
More informationThe SAS System 18:28 Saturday, March 10, Plot of Canonical Variables Identified by Cluster
The SAS System 18:28 Saturday, March 10, 2018 1 The FASTCLUS Procedure Replace=FULL Radius=0 Maxclusters=2 Maxiter=10 Converge=0.02 Initial Seeds Cluster SepalLength SepalWidth PetalLength PetalWidth 1
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 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 informationLecture 6: Single-classification multivariate ANOVA (k-group( MANOVA)
Lecture 6: Single-classification multivariate ANOVA (k-group( MANOVA) Rationale and MANOVA test statistics underlying principles MANOVA assumptions Univariate ANOVA Planned and unplanned Multivariate ANOVA
More informationST 732, Midterm Solutions Spring 2019
ST 732, Midterm Solutions Spring 2019 Please sign the following pledge certifying that the work on this test is your own: I have neither given nor received aid on this test. Signature: Printed Name: There
More informationNeuendorf MANOVA /MANCOVA. Model: X1 (Factor A) X2 (Factor B) X1 x X2 (Interaction) Y4. Like ANOVA/ANCOVA:
1 Neuendorf MANOVA /MANCOVA Model: X1 (Factor A) X2 (Factor B) X1 x X2 (Interaction) Y1 Y2 Y3 Y4 Like ANOVA/ANCOVA: 1. Assumes equal variance (equal covariance matrices) across cells (groups defined by
More information