Advantages of Mixed-effects Regression Models (MRM; aka multilevel, hierarchical linear, linear mixed models) 1. MRM explicitly models individual
|
|
- Branden Newman
- 5 years ago
- Views:
Transcription
1 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 measures (a) need not have same number of obs per subject (b) time can be continuous, rather than a fixed set of points 3. Flexible specification of the covariance structure among repeated measures methods for testing specific determinants of this structure 4. MRM can be extended to higher-level models repeated observations within individuals within clusters 5. Generalizations for non-normal data 1
2 2-level model for longitudinal data y i n i 1 = X i n i p β p 1 + Z i n i r υ i r 1 + ε i n i 1 y i = n i 1 response vector for individual i i = 1... N individuals j = 1... n i observations for individual i X i = n i p design matrix for the fixed effects β = p 1 vector of unknown fixed parameters Z i = n i r design matrix for the random effects υ i = r 1 vector of unknown random effects N (0, Σ υ ) ε i = n i 1 residual vector N (0, σ 2 I ni ) 2
3 Random-intercepts Model each subject is parallel to their group trend y = T ime + Grp + (Grp T ime) + Subj + Error y ij = β 0 + β 1 T ij + β 2 G i + β 3 (G i T ij ) + υ 0i + ε ij υ 0i N (0, σ 2 υ) ε ij N (0, σ 2 ) 3
4 Random Intercepts and Trend Model subjects deviate in terms of both intercept & slope y = T ime + Grp + (G T ) + Subj + (S T ) + Error y ij = β 0 + β 1 T ij + β 2 G i + β 3 (G i T ij ) + υ 0i + υ 1i T ij + ε ij υ 0i υ 1i N 0 0, σ 2 υ 0 σ υ0 υ 1 σ υ0 υ 1 σ 2 υ 1 ε ij N (0, σ 2 ) 4
5 Within-Unit / Between-Unit representation Within-subjects model - level 1 (j = 1,..., n i ) y ij = b 0i + b 1i X 1ij b p1i X p1ij + ε ij Between-subjects model - level 2 (i = 1,..., N) slopes as outcomes model b 0i = β 0 + β 0(2) x i + υ 0i b 1i = β 1 + β 1(2) x i + υ 1i. =... b p1i = β p1 + β p1(2) x i β = β 0 β 1... β p1 β 0(2) β 1(2)... β p1(2) intercept level-1 level-2 cross-level 5
6 Matrix form of model for individual i y i1 y i2... y ini y i n i 1 = 1 T ime i1 Group i Grp i T i1 1 T ime i2 Group i Grp i T i T ime ini Group i Grp i T ini X i n i p β 0 β 1 β 2 β 3 β p T ime i1 1 T ime i T ime ini Z i n i r υ 0i υ 1i υ i r 1 + ε i1 ε i2... ε ini ε i n i 1 Time might be years or months, and could differ for each subject 6
7 The conditional variance-covariance matrix is now of the form: Σy i = Z i Σ υ Z i + σ 2 I ni For example, with r = 2, n = 3, and Z i = the conditional variance-covariance Σy i = σ 2 I ni + σ 2 υ 0 σ 2 υ 0 + σ υ0 υ 1 σ 2 υ 0 + 2σ υ0 υ 1 σ 2 υ 0 + σ υ0 υ 1 σ 2 υ 0 + 2σ υ0 υ 1 + σ 2 υ 1 σ 2 υ 0 + 3σ υ0 υ 1 + 2σ 2 υ 1 σ 2 υ 0 + 2σ υ0 υ 1 σ 2 υ 0 + 3σ υ0 υ 1 + 2σ 2 υ 1 σ 2 υ 0 + 4σ υ0 υ 1 + 4σ 2 υ 1 variances and covariances change across time More general models allow autocorrelated errors, ε i N (0, σ 2 Ω i ), where Ω might represent AR or MA process 7
8 Estimation - EM algorithm opposite process of I do cocaine so I can work more, so I can do more cocaine, so I can work more, etc., Effect of increasing cocaine use Cocaine Work Health do cocaine work more declines do more cocaine work more declines more do even more cocaine work even more declines even more do a ton of cocaine always working death Effect of EM estimation of parameters M-Step (ML) E-Step (EB) Estimation starting values β, σ 2, Σ υ estimate υ i Σ υ yi improves re-estimate β, σ 2, Σ υ re-estimate υ i Σ υ yi improves more re-re-estimate β, σ 2, Σ υ re-re-estimate υ i Σ υ yi improves even more RE-estimate β, σ 2, Σ υ RE-estimate υ i Σ υ yi convergence EM is better than cocaine since EM leads to convergence and not death 8
9 EM solution - random intercepts model E-step (expectation - Expected A Posteriori or Empirical Bayes) 1 n i υ i = ρ ni n y i ij x ijβ j=1 n i σ 2 υ y i = σ 2 υ(1 ρ ni n i ) where ρ ni n i = M-step (maximization - Maximium Likelihood ) ˆβ = N i X ix i 1 N i X i(y i 1 i υ i ) n i r 1 + (n i 1)r and r = σ2 υ σ 2 υ + σ 2 ˆσ 2 υ = 1 N N i υ i 2 + σ 2 υ y i ˆσ 2 = ( N i n i) 1 N i (y i X i ˆβ 1i υ i ) (y i X i ˆβ 1i υ i ) + n i σ 2 υ y i provide starting values for β, σ 2 υ, and σ 2 perform E-step, perform M-step, repeat early and often (until convergence) 9
10 Example: Drug Plasma Levels and Clinical Response Riesby and associates (Riesby et al., 1977) examined the relationship between Imipramine (IMI) and Desipramine (DMI) plasma levels and clinical response in 66 depressed inpatients (37 endogenous and 29 non-endogenous) Drug-Washout day0 day7 day14 day21 day28 day35 wk 0 wk 1 wk 2 wk 3 wk 4 wk 5 Hamilton Depression HD 1 HD 2 HD 3 HD 4 HD 5 HD 6 Diagnosis Dx IMI IMI 3 IMI 4 IMI 5 IMI 6 DMI DMI 3 DMI 4 DMI 5 DMI 6 n
11 outcome variable Hamilton Depression Scores (HD) independent variables Dx, IMI and DMI Dx - endogenous (=1) or non-endogenous (=0) IM I (imipramine) drug-plasma levels (µg/l) antidepressant given 225 mg/day, weeks 3-6 DM I (desipramine) drug-plasma levels (µg/l) metabolite of imipramine 11
12 Descriptive Statistics Observed HDRS Means, n, and sd Washout wk 0 wk 1 wk 2 wk 3 wk 4 wk 5 Endog n Non-Endog n pooled sd
13 Correlations: n = 46 and 46 n 66 wk 0 wk 1 wk 2 wk 3 wk 4 wk 5 week week week week week week
14 increasing variance across time general linear decline over time 14
15 Examination of HD across all weeks HD i1 HD i2... HD ini y i n i 1 = 1 W EEK i1 1 W EEK i W EEK ini X i n i p β 0 β 1 β p W EEK i1 1 W EEK i W EEK ini υ 0i υ 1i + ε i1 ε i2... ε ini Z i υ i ε i n i r r 1 n i 1 where max(n i ) = 6, and X i = Z i =
16 Within-subjects and between-subjects components Within-subjects model HD ij = b 0i + b 1i T ime ij + RESID ij y ij = b 0i + b 1i x ij + ε ij b 0i b 1i = week 0 HD level for patient i = weekly change in HD for patient i i = patients j = 1... n i observations (max = 6) for patient i Between-subjects models b 0i = β 0 + υ 0i b 1i = β 1 + υ 1i β 0 β 1 υ 0i υ 1i = average week 0 HD level = average HD weekly improvement = individual deviation from average intercept = individual deviation from average improvement 16
17 parameter ML estimate se z p < β β σ 2 υ σ υ0 υ σ 2 υ σ log L = χ 2 2 = 66.1, p <.0001 for H 0 : σ υ0 υ 1 = συ 2 1 = 0 σ υ0 υ 1 as corr between intercept and slope = Wald tests are dubious for variance parameters, likelihood-ratio tests are preferred (though divide p-value by 2) Wald z-statistics sometimes expressed as χ 2 1 (by squaring z-value) 17
18 Observed and estimated means (= X ˆβ) wk 0 wk 1 wk 2 wk 3 wk 4 wk 5 n obs est
19 Obs. (pairwise) and est. variance-covariance matrix Σy = ˆΣy = Z ˆΣ υ Z + ˆσ 2 I Z = = ˆΣυ = note: from random-int model: ˆσ 2 υ = and ˆσ2 =
20 Empirical Bayes estimates of Subject Trends 20
21 Empirical Bayes estimates of Subject Trends 21
22 Examination of HD across all weeks by diagnosis HD i1 HD i2... HD ini y i n i 1 = 1 W EEK i1 Dx i Dx i W k i1 1 W EEK i2 Dx i Dx i W k i W EEK ini Dx i Dx i W k ini X i n i p β 0 β 1 β 2 β 3 β p W EEK i1 1 W EEK i W EEK ini υ 0i υ 1i + ε i1 ε i2... ε ini Z i υ i ε i n i r r 1 n i 1 where max(n i ) = 6, Z i = , Dx i = 0 for NE 1 for E 22
23 Within-subjects and between-subjects components Within-subjects model HD ij = b 0i + b 1i T ime ij + RESID ij b 0i b 1i = week 0 HD level for patient i = weekly change in HD for patient i Between-subjects models b 0i = β 0 + β 2 Dx i + υ 0i b 1i = β 1 + β 3 Dx i + υ 1i β 0 = average week 0 HD level for NE patients (Dx i = 0) β 1 = average HD weekly improvement for NE patients (Dx i = 0) β 2 = average week 0 HD difference for E patients β 3 = average HD weekly improvement difference for endogenous patients = individual deviation from average intercept = individual deviation from average improvement υ 0i υ 1i 23
24 parameter ML estimate se z p < NE int β NE slope β E int diff β E slope diff β σ 2 υ σ υ0 υ σ 2 υ σ log L = χ 2 2 = 4.1, p ns, compared to model with β 2 = β 3 = 0 σ β0 β 1 as corr between intercept and slope =
25 Riesby data - model fit by diagnosis 25
26 Examination of HD across all weeks - quadratic trend HD i1 HD i2... HD ini y i n i 1 = 1 W EEK i1 W EEKi1 2 1 W EEK i2 W EEKi W EEK ini W EEKin 2 i X i n i p β 0 β 1 β 2 β p W EEK i1 W EEKi1 2 1 W EEK i2 W EEKi W EEK ini W EEKin 2 i υ 0i υ 1i υ 2i + ε i1 ε i2... ε ini Z i υ i ε i n i r r 1 n i 1 where max(n i ) = 6, and X i = Z i =
27 Within-subjects and between-subjects components Within-subjects model HD ij y ij = b 0i + b 1i T ime ij + b 2i T ime 2 ij + RESID ij = b 0i + b 1i x ij + b 2i x 2 ij + ε ij b 0i b 1i b 2i = week 0 HD level for patient i = weekly linear change in HD for patient i = weekly quadratic change in HD for patient i Between-subjects models b 0i = β 0 + υ 0i b 1i = β 1 + υ 1i b 2i = β 2 + υ 2i β 0 β 1 β 2 υ 0i υ 1i υ 2i = average week 0 HD level = average HD weekly linear change = average HD weekly quadratic change = individual deviation from average intercept = individual deviation from average linear change = individual deviation from average quadratic change 27
28 parameter ML estimate se z p < β β β συ σ υ0 υ συ σ υ0 υ σ υ1 υ συ σ log L = χ 2 4 = 11.4, p < 0.025, compared to model with β 2 = σ 2 υ 2 = σ υ0 υ 2 = σ υ1 υ 2 = 0 χ 2 3 = 11.0, p < 0.02, compared to model with σ2 υ 2 = σ υ0 υ 2 = σ υ1 υ 2 = 0 σ υ1 υ 2 as corr between linear and quadratic terms =
29 Average linear and individual quadratic trends 29
30 Observed (pairwise) and estimated variance-covariance matrix Σy = ˆΣy = Z ˆΣ υ Z + ˆσ 2 I = where Z = ˆΣ υ =
31 Empirical Bayes estimates of Subject Trends 31
32 SAS PROC MIXED formulation y n 1 = X n p β p 1 + Z n Sr υ Sr 1 + ε n 1 i = 1... S individuals n = ( S i n i ) total number of observations y = vector obtained by stacking y i vectors X = matrix obtained by stacking X i matrices υ = vector obtained by stacking υ i vectors ε = vector obtained by stacking ε i vectors Z = block diagonal matrix with Z i on diagonal Z = Z Z Z S 32
33 V ar υ ε = G 0 0 R V ar(y) = Z G Z + R n n n Sr Sr Sr Sr n n n G = Σ υ Σ υ Σ υs Can model variance/covariance of y in terms of: G - random-effects only (R = σ 2 I n ) R - variance/covariance structures unstructured, or AR(1), or Toeplitz (banded) G and R - random-effects with correlated errors 33
34 Example 4a: Analysis of Riesby dataset using MRM. This example has a few different PROC MIXED specifications, and includes a grouping variable and curvilinear effect of time. (SAS code and output) hedeker/riesbym.txt 34
35 SAS MIXED code - RIESBYM.SAS TITLE1 analysis of riesby data - hdrs scores across time ; DATA one; INFILE c:\mixdemo\riesby.dat ; INPUT id hamd intcpt week endog endweek ; PROC FORMAT; VALUE endog 0= nonendog 1= endog ; VALUE week 0= week 0 1= week 1 2= week 2 3= week 3 4= week 4 5= week 5 ; PROC MIXED METHOD=ML COVTEST; CLASS id; MODEL hamd = week /SOLUTION; RANDOM INTERCEPT /SUB=id TYPE=UN G; TITLE2 random intercepts model: compound symmetry structure ; PROC MIXED METHOD=ML COVTEST; CLASS id; MODEL hamd = week /SOLUTION; RANDOM INTERCEPT week /SUB=id TYPE=UN G GCORR; TITLE2 random trend model ; 35
36 PROC MIXED METHOD=ML COVTEST; CLASS id; MODEL hamd = week endog endweek /SOLUTION; RANDOM INTERCEPT week /SUB=id TYPE=UN G GCORR; TITLE2 random trend model with group effects ; PROC MIXED METHOD=ML COVTEST; CLASS id; MODEL hamd = week week*week /SOLUTION; RANDOM INTERCEPT week week*week /SUB=id TYPE=UN G GCORR; TITLE2 random quadratic trend model ; RUN; 36
37 Riesby.dat - data from a few subjects
38 Example 4b: Analysis of Riesby dataset. This handout shows how empirical Bayes estimates can be output to a dataset in order to calculate estimated individual scores at all timepoints. (SAS code and output) hedeker/riesbym2.txt 38
39 SAS MIXED code - RIESBYM2.SAS TITLE1 analysis of riesby data - empirical bayes estimates ; DATA one; INFILE c:\mixdemo\riesby.dat ; INPUT id hamd intcpt week endog endweek ; PROC FORMAT; VALUE endog 0= nonendog 1= endog ; VALUE week 0= week 0 1= week 1 2= week 2 3= week 3 4= week 4 5= week 5 ; PROC MIXED METHOD=ML; CLASS id; MODEL hamd = week /SOLUTION; RANDOM INTERCEPT week /SUB=id TYPE=UN G S; ODS LISTING EXCLUDE SOLUTIONR; ODS OUTPUT SOLUTIONR=randest; TITLE2 random trend model ; /* print out the estimated random effects dataset */ PROC PRINT DATA=randest; 39
40 /* get a printout of the data in multivariate form */ PROC SORT DATA=one; BY id; DATA t0;set one; IF week=0; hamd 0 = hamd; DATA t1;set one; IF week=1; hamd 1 = hamd; DATA t2;set one; IF week=2; hamd 2 = hamd; DATA t3;set one; IF week=3; hamd 3 = hamd; DATA t4;set one; IF week=4; hamd 4 = hamd; DATA t5;set one; IF week=5; hamd 5 = hamd; DATA comp (KEEP=id hamd 0-hamd 5); MERGE t0 t1 t2 t3 t4 t5; BY id; PROC PRINT DATA=comp; VAR id hamd 0-hamd 5; 40
41 /* extract the intercepts and slopes for each person */ /* and compute the estimated hamd values across time */ PROC SORT DATA=randest; BY id; DATA randest2 (KEEP=id intdev slopedev int slope hdest 0-hdest 5); ARRAY y(2) intdev slopedev; DO par = 1 TO 2; SET randest; BY id; y(par) = ESTIMATE; IF par = 2 THEN DO; int = intdev; slope = slopedev; hdest 0 = int; hdest 1 = int + slope; hdest 2 = int + 2*slope; hdest 3 = int + 3*slope; hdest 4 = int + 4*slope; hdest 5 = int + 5*slope; END; IF LAST.id THEN RETURN; END; 41
42 PROC PRINT DATA=randest2; VAR id hdest 0-hdest 5; PROC PLOT DATA=randest2; PLOT intdev * slopedev; PLOT int * slope; TITLE2 plot of individual intercepts versus slopes ; RUN; 42
43 Time-varying Covariates - WS and BS effects Section in Hedeker & Gibbons (2006), Longitudinal Data Analysis, Wiley. 43
44 Examination of HD across 4 weeks by plasma drug-levels HD i1 HD i2... HD ini y i n i 1 = 1 W EEK i1 lnimi i1 lndmi i1 1 W EEK i2 lnimi i2 lndmi i W EEK ini lnimi ini lndmi ini X i n i p β 0 β 1 β 2 β 3 β p W EEK i1 1 W EEK i W EEK ini υ 0i υ 1i + ε i1 ε i2... ε ini Z i υ i ε i n i r r 1 n i 1 where max(n i ) = 4, and Z i =
45 Within-subjects and between-subjects components Within-subjects model HD ij = b 0i + b 1i T ij + b 2i ln IMI ij + b 3i ln DMI ij + Res ij b 0i = week 2 HD level for patient i with both ln IMI and ln DMI = 0 b 1i = weekly change in HD for patient i b 2i = change in HD due to ln IMI = change in HD due to ln DMI b 3i Between-subjects models b 0i = β 0 + υ 0i b 1i = β 1 + υ 1i b 2i = β 2 b 3i = β 3 45
46 β 0 β 1 β 2 β 3 υ 0i υ 1i = average week 2 HD level for drug-free patients = average HD weekly improvement = average HD difference for unit change in ln IMI = average HD difference for unit change in ln DMI = individual intercept deviation from model = individual slope deviation from model Here, week 2 is the actual study week (i.e., one week after the drug washout period), which is coded as 0 in this analysis of the last four study timepoints 46
47 parameter ML estimate se z p < int β slope β ln IMI β ln DMI β σ 2 υ σ υ0 υ σ 2 υ σ log L = σ υ0 υ 1 as corr between intercept and slope =
48 parameter estimate se p < HD total score intercept β slope β Baseline HD β ln IMI β ln DMI β συ σ υ0 υ συ σ HD change from baseline intercept β ns slope β ln IMI β ns ln DMI β συ σ υ0 υ συ σ
49 Correlation between HD scores and plasma levels (ln units) HD total score week 2 week 3 week 4 week 5 IMI DMI HD change from baseline week 2 week 3 week 4 week 5 IMI DMI p <
50 Model with time-varying covariate X ij Within-subjects model Y ij = b 0i + b 1i T ij + b 2i X ij + E ij Between-subjects models b 0i = β 0 + υ 0i b 1i = β 1 + υ 1i b 2i = β 2 Is the effect of X ij purely within-subjects? What about X ij = X ij + X i X i = X i + (X ij X i ) X i is between-subjects component of X X ij X i is within-subjects component of X 50
51 Time-varying covariate effects: purely between-subjects 51
52 Time-varying covariate effects: purely within-subjects 52
53 Model with decomposition of time-varying covariate X ij Within-subjects model Y ij = b 0i + b 1i T ij + b 2i (X ij X i ) + E ij Between-subjects models b 0i = β 0 + β BS Xi + υ 0i b 1i = β 1 + υ 1i b 2i = β W S Notice, effect of X is now β BS Xi + β W S (X ij X i ) β BS = effect of Xi on Ȳi BS or cross-sectional β W S = effect of (X ij X i ) on (Y ij Ȳi) WS or longitudinal 53
54 Model with only X ij assumes equal BS and WS effects (β BS = β W S ) suppose β BS = β W S = β, then in the model with decomposition, the effect of X ij = β Xi + β (X ij X i ) = β X ij precisely what the model with only X ij assumes Equal WS and BS effects of X ij? can be a dubious assumption needs to be tested (by comparing two models via LR test) there is no guarantee that β BS and β W S even agree on sign 54
55 Time-varying covariate effects: opposite sign WS and BS effects 55
56 parameter estimate se p < assuming BS=WS drug effects intercept ns slope ln IMI ns ln DMI deviance = relaxing BS=WS drug effects intercept ns slope ln IMI BS ns ln DMI BS ln IMI WS ns ln DMI WS ns deviance = X 2 2 = = 3 Accept H 0 : β BS = β W S 56
57 Example 4c: Analysis of Riesby dataset. This handout has the analysis considering the time-varying drug plasma levels, separating the within-subjects from the between-subjects effects for these time-varying covariates. (SAS code and output) hedeker/riesbsws.txt 57
58 SAS MIXED code - RIESBSWS.SAS TITLE1 partitioning BS and WS effects of drug levels ; DATA one; INFILE c:\mixdemo\riesbyt4.dat ; INPUT id hamdelt intcpt week sex endog lnimi lndmi ; PROC SORT; BY id; PROC MEANS NOPRINT; CLASS id; VAR lnimi lndmi; OUTPUT OUT = two MEAN = mlnimi mlndmi; DATA three; MERGE one two; BY id; lnidev = lnimi - mlnimi; lnddev = lndmi - mlndmi; PROC MIXED METHOD=ML COVTEST; CLASS id; MODEL hamdelt = week lnimi lndmi /SOLUTION; RANDOM INTERCEPT week /SUB=id TYPE=UN G GCORR; TITLE2 assuming bs=ws drug effects ; PROC MIXED METHOD=ML COVTEST; CLASS id; MODEL hamdelt = week mlnimi mlndmi lnidev lnddev /SOLUTION; RANDOM INTERCEPT week /SUB=id TYPE=UN G GCORR; TITLE2 relaxing bs=ws drug effects ; 58
Missing Data in Longitudinal Studies: Mixed-effects Pattern-Mixture and Selection Models
Missing Data in Longitudinal Studies: Mixed-effects Pattern-Mixture and Selection Models Hedeker D & Gibbons RD (1997). Application of random-effects pattern-mixture models for missing data in longitudinal
More informationMixed Models for Longitudinal Binary Outcomes. Don Hedeker Department of Public Health Sciences University of Chicago.
Mixed Models for Longitudinal Binary Outcomes Don Hedeker Department of Public Health Sciences University of Chicago hedeker@uchicago.edu https://hedeker-sites.uchicago.edu/ Hedeker, D. (2005). Generalized
More informationMixed-effects Polynomial Regression Models chapter 5
Mixed-effects Polynomial Regression Models chapter 5 1 Figure 5.1 Various curvilinear models: (a) decelerating positive slope; (b) accelerating positive slope; (c) decelerating negative slope; (d) accelerating
More informationLongitudinal Data Analysis via Linear Mixed Model
Longitudinal Data Analysis via Linear Mixed Model Edps/Psych/Stat 587 Carolyn J. Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Fall 2017 Outline Introduction
More informationMixed-Effects Pattern-Mixture Models for Incomplete Longitudinal Data. Don Hedeker University of Illinois at Chicago
Mixed-Effects Pattern-Mixture Models for Incomplete Longitudinal Data Don Hedeker University of Illinois at Chicago This work was supported by National Institute of Mental Health Contract N44MH32056. 1
More informationSection IV MODELS FOR MULTILEVEL DATA
Section IV MODELS FOR MULTILEVEL DATA Chapter 12 An Introduction to Growth Modeling Donald Hedeker 12.1. Introduction Longitudinal studies are increasingly common in social sciences research. In these
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 informationGEE for Longitudinal Data - Chapter 8
GEE for Longitudinal Data - Chapter 8 GEE: generalized estimating equations (Liang & Zeger, 1986; Zeger & Liang, 1986) extension of GLM to longitudinal data analysis using quasi-likelihood estimation method
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 informationFixed effects results...32
1 MODELS FOR CONTINUOUS OUTCOMES...7 1.1 MODELS BASED ON A SUBSET OF THE NESARC DATA...7 1.1.1 The data...7 1.1.1.1 Importing the data and defining variable types...8 1.1.1.2 Exploring the data...12 Univariate
More informationJournal of Statistical Software
JSS Journal of Statistical Software March 2013, Volume 52, Issue 12. http://www.jstatsoft.org/ MIXREGLS: A Program for Mixed-Effects Location Scale Analysis Donald Hedeker University of Illinois at Chicago
More informationMixed Models for Longitudinal Ordinal and Nominal Outcomes
Mixed Models for Longitudinal Ordinal and Nominal Outcomes Don Hedeker Department of Public Health Sciences Biological Sciences Division University of Chicago hedeker@uchicago.edu Hedeker, D. (2008). Multilevel
More informationAn Application of a Mixed-Effects Location Scale Model for Analysis of Ecological Momentary Assessment (EMA) Data
An Application of a Mixed-Effects Location Scale Model for Analysis of Ecological Momentary Assessment (EMA) Data Don Hedeker, Robin Mermelstein, & Hakan Demirtas University of Illinois at Chicago hedeker@uic.edu
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 informationAnalysis of Repeated Measures and Longitudinal Data in Health Services Research
Analysis of Repeated Measures and Longitudinal Data in Health Services Research Juned Siddique, Donald Hedeker, and Robert D. Gibbons Abstract This chapter reviews statistical methods for the analysis
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 informationWhy analyze as ordinal? Mixed Models for Longitudinal Ordinal Data Don Hedeker University of Illinois at Chicago
Why analyze as ordinal? Mixed Models for Longitudinal Ordinal Data Don Hedeker University of Illinois at Chicago hedeker@uic.edu www.uic.edu/ hedeker/long.html Efficiency: Armstrong & Sloan (1989, Amer
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 informationReview. Timothy Hanson. Department of Statistics, University of South Carolina. Stat 770: Categorical Data Analysis
Review Timothy Hanson Department of Statistics, University of South Carolina Stat 770: Categorical Data Analysis 1 / 22 Chapter 1: background Nominal, ordinal, interval data. Distributions: Poisson, binomial,
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 informationLongitudinal Modeling with Logistic Regression
Newsom 1 Longitudinal Modeling with Logistic Regression Longitudinal designs involve repeated measurements of the same individuals over time There are two general classes of analyses that correspond to
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 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 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 informationInternational Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: , ISSN(Online): Vol.9, No.9, pp , 2016
International Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: 0974-4304, ISSN(Online): 2455-9563 Vol.9, No.9, pp 477-487, 2016 Modeling of Multi-Response Longitudinal DataUsing Linear Mixed Modelin
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 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 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 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 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 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 informationMultilevel Modeling of Non-Normal Data. Don Hedeker Department of Public Health Sciences University of Chicago.
Multilevel Modeling of Non-Normal Data Don Hedeker Department of Public Health Sciences University of Chicago email: hedeker@uchicago.edu https://hedeker-sites.uchicago.edu/ Hedeker, D. (2005). Generalized
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 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 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 informationGrowth Mixture Model
Growth Mixture Model Latent Variable Modeling and Measurement Biostatistics Program Harvard Catalyst The Harvard Clinical & Translational Science Center Short course, October 28, 2016 Slides contributed
More informationMixed Models for Longitudinal Ordinal and Nominal Outcomes. Don Hedeker Department of Public Health Sciences University of Chicago
Mixed Models for Longitudinal Ordinal and Nominal Outcomes Don Hedeker Department of Public Health Sciences University of Chicago hedeker@uchicago.edu https://hedeker-sites.uchicago.edu/ Hedeker, D. (2008).
More informationCategorical Predictor Variables
Categorical Predictor Variables We often wish to use categorical (or qualitative) variables as covariates in a regression model. For binary variables (taking on only 2 values, e.g. sex), it is relatively
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 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 informationModeling Longitudinal Count Data with Excess Zeros and Time-Dependent Covariates: Application to Drug Use
Modeling Longitudinal Count Data with Excess Zeros and : Application to Drug Use University of Northern Colorado November 17, 2014 Presentation Outline I and Data Issues II Correlated Count Regression
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 12 1 / 34 Correlated data multivariate observations clustered data repeated measurement
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 informationThe lmm Package. May 9, Description Some improved procedures for linear mixed models
The lmm Package May 9, 2005 Version 0.3-4 Date 2005-5-9 Title Linear mixed models Author Original by Joseph L. Schafer . Maintainer Jing hua Zhao Description Some improved
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 informationMIXED-EFFECTS POLYNOMIAL REGRESSION MODELS
CHAPTER 5 MIXED-EFFECTS POLYNOMIAL REGRESSION MODELS 5.1 INTRODUCTION In many situations, it is too simplistic to assume that the change across time is linear. For example, it may be that the outcome changes
More informationBIOS 2083: Linear Models
BIOS 2083: Linear Models Abdus S Wahed September 2, 2009 Chapter 0 2 Chapter 1 Introduction to linear models 1.1 Linear Models: Definition and Examples Example 1.1.1. Estimating the mean of a N(μ, σ 2
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 informationContrasting Marginal and Mixed Effects Models Recall: two approaches to handling dependence in Generalized Linear Models:
Contrasting Marginal and Mixed Effects Models Recall: two approaches to handling dependence in Generalized Linear Models: Marginal models: based on the consequences of dependence on estimating model parameters.
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 informationPackage mlmmm. February 20, 2015
Package mlmmm February 20, 2015 Version 0.3-1.2 Date 2010-07-07 Title ML estimation under multivariate linear mixed models with missing values Author Recai Yucel . Maintainer Recai Yucel
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 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 informationPackage lmm. R topics documented: March 19, Version 0.4. Date Title Linear mixed models. Author Joseph L. Schafer
Package lmm March 19, 2012 Version 0.4 Date 2012-3-19 Title Linear mixed models Author Joseph L. Schafer Maintainer Jing hua Zhao Depends R (>= 2.0.0) Description Some
More informationML estimation: Random-intercepts logistic model. and z
ML estimation: Random-intercepts logistic model log p ij 1 p = x ijβ + υ i with υ i N(0, συ) 2 ij Standardizing the random effect, θ i = υ i /σ υ, yields log p ij 1 p = x ij β + σ υθ i with θ i N(0, 1)
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 informationModel Assumptions; Predicting Heterogeneity of Variance
Model Assumptions; Predicting Heterogeneity of Variance Today s topics: Model assumptions Normality Constant variance Predicting heterogeneity of variance CLP 945: Lecture 6 1 Checking for Violations of
More informationWeighted Least Squares
Weighted Least Squares The standard linear model assumes that Var(ε i ) = σ 2 for i = 1,..., n. As we have seen, however, there are instances where Var(Y X = x i ) = Var(ε i ) = σ2 w i. Here w 1,..., w
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 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 informationModule 6 Case Studies in Longitudinal Data Analysis
Module 6 Case Studies in Longitudinal Data Analysis Benjamin French, PhD Radiation Effects Research Foundation SISCR 2018 July 24, 2018 Learning objectives This module will focus on the design of longitudinal
More informationResearch Design: Topic 18 Hierarchical Linear Modeling (Measures within Persons) 2010 R.C. Gardner, Ph.d.
Research Design: Topic 8 Hierarchical Linear Modeling (Measures within Persons) R.C. Gardner, Ph.d. General Rationale, Purpose, and Applications Linear Growth Models HLM can also be used with repeated
More informationOutline. Linear OLS Models vs: Linear Marginal Models Linear Conditional Models. Random Intercepts Random Intercepts & Slopes
Lecture 2.1 Basic Linear LDA 1 Outline Linear OLS Models vs: Linear Marginal Models Linear Conditional Models Random Intercepts Random Intercepts & Slopes Cond l & Marginal Connections Empirical Bayes
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 informationFigure 36: Respiratory infection versus time for the first 49 children.
y BINARY DATA MODELS We devote an entire chapter to binary data since such data are challenging, both in terms of modeling the dependence, and parameter interpretation. We again consider mixed effects
More informationRandom Coefficients Model Examples
Random Coefficients Model Examples STAT:5201 Week 15 - Lecture 2 1 / 26 Each subject (or experimental unit) has multiple measurements (this could be over time, or it could be multiple measurements on a
More informationIndependence (Null) Baseline Model: Item means and variances, but NO covariances
CFA Example Using Forgiveness of Situations (N = 1103) using SAS MIXED (See Example 4 for corresponding Mplus syntax and output) SAS Code to Read in Mplus Data: * Import data from Mplus, becomes var1-var23
More informationWeighted Least Squares
Weighted Least Squares The standard linear model assumes that Var(ε i ) = σ 2 for i = 1,..., n. As we have seen, however, there are instances where Var(Y X = x i ) = Var(ε i ) = σ2 w i. Here w 1,..., w
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 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 informationDesigning Multilevel Models Using SPSS 11.5 Mixed Model. John Painter, Ph.D.
Designing Multilevel Models Using SPSS 11.5 Mixed Model John Painter, Ph.D. Jordan Institute for Families School of Social Work University of North Carolina at Chapel Hill 1 Creating Multilevel Models
More informationEconometrics of Panel Data
Econometrics of Panel Data Jakub Mućk Meeting # 2 Jakub Mućk Econometrics of Panel Data Meeting # 2 1 / 26 Outline 1 Fixed effects model The Least Squares Dummy Variable Estimator The Fixed Effect (Within
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 informationApplication of Item Response Theory Models for Intensive Longitudinal Data
Application of Item Response Theory Models for Intensive Longitudinal Data Don Hedeker, Robin Mermelstein, & Brian Flay University of Illinois at Chicago hedeker@uic.edu Models for Intensive Longitudinal
More informationAnalysis of Longitudinal Data. Patrick J. Heagerty PhD Department of Biostatistics University of Washington
Analsis of Longitudinal Data Patrick J. Heagert PhD Department of Biostatistics Universit of Washington 1 Auckland 2008 Session Three Outline Role of correlation Impact proper standard errors Used to weight
More informationClassification. Chapter Introduction. 6.2 The Bayes classifier
Chapter 6 Classification 6.1 Introduction Often encountered in applications is the situation where the response variable Y takes values in a finite set of labels. For example, the response Y could encode
More informationSample Size and Power Considerations for Longitudinal Studies
Sample Size and Power Considerations for Longitudinal Studies Outline Quantities required to determine the sample size in longitudinal studies Review of type I error, type II error, and power For continuous
More informationProblems. Suppose both models are fitted to the same data. Show that SS Res, A SS Res, B
Simple Linear Regression 35 Problems 1 Consider a set of data (x i, y i ), i =1, 2,,n, and the following two regression models: y i = β 0 + β 1 x i + ε, (i =1, 2,,n), Model A y i = γ 0 + γ 1 x i + γ 2
More informationFor more information about how to cite these materials visit
Author(s): Kerby Shedden, Ph.D., 2010 License: Unless otherwise noted, this material is made available under the terms of the Creative Commons Attribution Share Alike 3.0 License: http://creativecommons.org/licenses/by-sa/3.0/
More informationWell-developed and understood properties
1 INTRODUCTION TO LINEAR MODELS 1 THE CLASSICAL LINEAR MODEL Most commonly used statistical models Flexible models Well-developed and understood properties Ease of interpretation Building block for more
More informationLDA Midterm Due: 02/21/2005
LDA.665 Midterm Due: //5 Question : The randomized intervention trial is designed to answer the scientific questions: whether social network method is effective in retaining drug users in treatment programs,
More information,..., θ(2),..., θ(n)
Likelihoods for Multivariate Binary Data Log-Linear Model We have 2 n 1 distinct probabilities, but we wish to consider formulations that allow more parsimonious descriptions as a function of covariates.
More informationModelling the Covariance
Modelling the Covariance Jamie Monogan Washington University in St Louis February 9, 2010 Jamie Monogan (WUStL) Modelling the Covariance February 9, 2010 1 / 13 Objectives By the end of this meeting, participants
More informationModelling heterogeneous variance-covariance components in two-level multilevel models with application to school effects educational research
Modelling heterogeneous variance-covariance components in two-level multilevel models with application to school effects educational research Research Methods Festival Oxford 9 th July 014 George Leckie
More informationmultilevel modeling: concepts, applications and interpretations
multilevel modeling: concepts, applications and interpretations lynne c. messer 27 october 2010 warning social and reproductive / perinatal epidemiologist concepts why context matters multilevel models
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 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 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 informationMonday 7 th Febraury 2005
Monday 7 th Febraury 2 Analysis of Pigs data Data: Body weights of 48 pigs at 9 successive follow-up visits. This is an equally spaced data. It is always a good habit to reshape the data, so we can easily
More informationProduct Held at Accelerated Stability Conditions. José G. Ramírez, PhD Amgen Global Quality Engineering 6/6/2013
Modeling Sub-Visible Particle Data Product Held at Accelerated Stability Conditions José G. Ramírez, PhD Amgen Global Quality Engineering 6/6/2013 Outline Sub-Visible Particle (SbVP) Poisson Negative Binomial
More informationLongitudinal Data Analysis of Health Outcomes
Longitudinal Data Analysis of Health Outcomes Longitudinal Data Analysis Workshop Running Example: Days 2 and 3 University of Georgia: Institute for Interdisciplinary Research in Education and Human Development
More informationLinear Regression. In this problem sheet, we consider the problem of linear regression with p predictors and one intercept,
Linear Regression In this problem sheet, we consider the problem of linear regression with p predictors and one intercept, y = Xβ + ɛ, where y t = (y 1,..., y n ) is the column vector of target values,
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 informationChapter 2 The Simple Linear Regression Model: Specification and Estimation
Chapter The Simple Linear Regression Model: Specification and Estimation Page 1 Chapter Contents.1 An Economic Model. An Econometric Model.3 Estimating the Regression Parameters.4 Assessing the Least Squares
More informationChapter 1: Linear Regression with One Predictor Variable also known as: Simple Linear Regression Bivariate Linear Regression
BSTT523: Kutner et al., Chapter 1 1 Chapter 1: Linear Regression with One Predictor Variable also known as: Simple Linear Regression Bivariate Linear Regression Introduction: Functional relation between
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 informationStatistical Methods III Statistics 212. Problem Set 2 - Answer Key
Statistical Methods III Statistics 212 Problem Set 2 - Answer Key 1. (Analysis to be turned in and discussed on Tuesday, April 24th) The data for this problem are taken from long-term followup of 1423
More informationJoint Modeling of Longitudinal Item Response Data and Survival
Joint Modeling of Longitudinal Item Response Data and Survival Jean-Paul Fox University of Twente Department of Research Methodology, Measurement and Data Analysis Faculty of Behavioural Sciences Enschede,
More informationMaster s Written Examination - Solution
Master s Written Examination - Solution Spring 204 Problem Stat 40 Suppose X and X 2 have the joint pdf f X,X 2 (x, x 2 ) = 2e (x +x 2 ), 0 < x < x 2
More informationSpecifying Latent Curve and Other Growth Models Using Mplus. (Revised )
Ronald H. Heck 1 University of Hawai i at Mānoa Handout #20 Specifying Latent Curve and Other Growth Models Using Mplus (Revised 12-1-2014) The SEM approach offers a contrasting framework for use in analyzing
More information