Covariance Models (*) X i : (n i p) design matrix for fixed effects β : (p 1) regression coefficient for fixed effects
|
|
- Myrtle Gordon
- 5 years ago
- Views:
Transcription
1 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 effects Note: see pg 60 for specific examples Note: FLW Appendix A = Gentle Intro to Matrices 139 Heagerty, 2006
2 Covariance Models (*) Mixed Models Laird & Ware (1982) Z i : (n i q) design matrix for random effects b i : (q 1) vector of random effects e i : (n i 1) vector of errors For the random components of the model we typically assume: b i N (0, D) e i N (0, R i ) 140 Heagerty, 2006
3 140-1 Heagerty, 2006 Laird & Ware Chair, Dept Biostatistics HSPH Associate Dean HSPH
4 LMM and components of variation (*) This yields a covariance structure: cov(y i ) = Z i DZ T i }{{} + R i }{{} between-cluster var + within-cluster var We assume that observations on different subjects are independent Note: This is a matrix (compact) way of writing the covariance for any possibe pair Y ij, Y ik, and represents the variance and covariance details that we presented on pp 60-1 and Heagerty, 2006
5 LMM and components of variation Within-Subject: Independence Model : R i = σ 2 I or general diagonal matrix Then, assuming normal errors we have that Y i = (Y i1, Y i2,, Y i,ni ) are conditionally independent given b i This model assumes that the within-subject errors do not have any serial correlation 142 Heagerty, 2006
6 More on Covariance Models Within-Subject: Serial Models Linear mixed models assume that each subject follows his/her own line In some situations the dependence is more local meaning that observations close in time are more similar than those far apart in time One model that we introduced is called the autoregressive model where: cov(e ij, e ik ) = σ 2 ρ j k 143 Heagerty, 2006
7 ID=1, rho=09 ID=2, rho=09 ID=3, rho=09 ID=4, rho=09 y y y y months months months months ID=5, rho=09 ID=6, rho=09 ID=7, rho=09 ID=8, rho=09 y y y y months months months months ID=9, rho=09 ID=10, rho=09 ID=11, rho=09 ID=12, rho=09 y y y y months months months months ID=13, rho=09 ID=14, rho=09 ID=15, rho=09 ID=16, rho=09 y y y y months months months months ID=17, rho=09 ID=18, rho=09 ID=19, rho=09 ID=20, rho=09 y y y y months months months months Heagerty, 2006
8 More on Covariance Models Autoregressive Correlation Assume t ij = j, n i 4: 1 ρ ρ 2 ρ 3 ρ 1 ρ ρ 2 corr(e i ) = ρ 2 ρ 1 ρ ρ 3 ρ 2 ρ Heagerty, 2006
9 More on Covariance Models Autoregressive Correlation Assume t ij = j: corr(e i ) = 1 ρ ρ 2 ρ (n 1) ρ 1 ρ ρ (n 2) ρ 2 ρ 1 ρ (n 3) ρ (n 1) ρ (n 2) ρ (n 3) Heagerty, 2006
10 More on Covariance Models Autoregressive Correlation corr(e i ) = Assume t ij unique: 1 ρ t i1 t i2 ρ t i1 t i3 ρ t i1 t in ρ t i2 t i1 1 ρ t i2 t i3 ρ t i2 t in ρ t i3 t i1 ρ t i3 t i2 1 ρ t i3 t in ρ t in t i1 ρ t in t i2 ρ t in t i Heagerty, 2006
11 More on Covariance Models Mixed + Serial Diggle (1988) proposed the following model Y ij = X ij β + b i,0 + W i (t ij ) + ɛ ij 147 Heagerty, 2006
12 Covariance Models Mixed + Serial The most general type of covariance model will combine some random effects with some additional aspects that characterize within-subject serial correlation One such model contains three sources of random variation: random intercept b i,0 serial process W i (t ij ) measurement error ɛ ij 148 Heagerty, 2006
13 We assume: var(b i,0 ) = ν 2 cov[w (s), W (t)] = σ 2 ρ s t var(ɛ ij ) = τ 2 Then: Total Variance = ν 2 + σ 2 + τ 2 Covariance(Y ij, Y ik ) = ν 2 + σ 2 ρ t ij t ik 149 Heagerty, 2006
14 Covariance Models Mixed + Serial of variation? Q: How to biologically interpret these three sources random intercept: This represents a trait of the subject FEV1 child size not captured by age and height CD4 subject s normal steady-state level serial variation: This represents a state for the subject FEV1 child current health status (infected with PseudoA) CD4 subject s current immune status (diet? treatment?) measurement error: This represents the instrumentation or process used to generate the final quantitative measurement FEV1 result of only one trial with expiration CD4 blood sample, lab processing 150 Heagerty, 2006
15 150-1 Heagerty, 2006
16 EDA for Covariance Structure Numerical Summaries Empirical covariance & correlation Variogram Define: R ij = Y ij X ij β = b i,0 + W (t ij ) + ɛ ij 151 Heagerty, 2006
17 Note: var(r ij ) = ν 2 + σ 2 + τ 2 E [ ] 1 2 (R ij R ik ) 2 = σ 2 (1 ρ t ij t ik ) + τ 2 R ij R ik = (b i,0 + W i (t ij ) + ɛ ij ) (b i,0 + W i (t ik ) + ɛ ik ) = [W i (t ij ) W i (t ik )] + [ɛ ij ɛ ik ] Plot: 1 2 ( R ij R ik ) 2 versus t ij t ik 152 Heagerty, 2006
18 152-1 Heagerty, 2006 Variogram E( dr^2 ) Total Variance deltatime
19 152-2 Heagerty, 2006 Variogram: Key Features When t ij = t ik : E [ ] 1 2 (R ij R ik ) 2 = σ 2 (1 ρ t ij t ik ) + τ 2 = σ 2 (1 ρ 0 ) + τ 2 = τ 2 = measurement error variance
20 152-3 Heagerty, 2006 Variogram: Key Features When t ij >> t ik : E [ ] 1 2 (R ij R ik ) 2 (large time separation) = σ 2 (1 ρ t ij t ik ) + τ 2 = σ 2 (1 ρ ) + τ 2 = σ 2 + τ 2 = serial and measurement error variances
21 dt dr FEV1 residual variogram Heagerty, 2006
22 Recall Simple Linear Regression (**) In simple linear regression we fit the model E(Y i X i ) = β 0 + β 1 X i We can write the estimate of the slope, β 1 as follows: β 1 = 1 i (X (X i X) 2 i X) (Y i Y ) This method is sometimes called ordinary least squares, or OLS i 153 Heagerty, 2006
23 Recall Simple Linear Regression (**) In some applications we still want to fit the regression model: E(Y i X i ) = β 0 + β 1 X i But now we want to assign weights, w i, to each observation Using the weights leads to weighted least squares (WLS) We can write the estimate of the slope, β 1 as follows: β 1 (w) = 1 i w (X i (X i X) 2 i X) w i (Y i Y ) With longitudinal data we have a method of estimation that generalizes this to allow covariance weights i 154 Heagerty, 2006
24 Estimation of β with known Σ i (**) Weighted least squares: In univariate regression, WLS yields estimates of β that minimize the objective function Q(β) = N w i (Y i X i β) 2 where W i is an (n i n i ) positive definite symmetric matrix i=1 Analogously, the multivariate version of WLS finds the value of the parameter β(w ) that minimizes Q W (β) = N (Y i X i β) T W i (Y i X i β) i=1 155 Heagerty, 2006
25 Estimation of β with known Σ i (**) It s straight forward to see that U(β) = β Q W (β) = 2 N X T i W i (Y i X i β) i=1 This is a general way to statistically define the regression estimator a solution to equations In general W i is chosen as the inverse of Σ i 156 Heagerty, 2006
26 The solution to the minimization solves U(β) = 0 and yields ( N ) 1 ( N ) β(w ) = X T i W i X i X T i W i Y i i=1 i=1 157 Heagerty, 2006
27 Properties of β(w ) (**) Given X 1, X 2, X N and W 1, W 2, W N ( ] N ) 1 ( N ) E [ β(w ) = X T i W i X i X T i W i E[Y i ] i=1 i=1 ( N ) 1 ( N ) = X T i W i X i X T i W i X i β = β i=1 i=1 Notice that the estimate β(w ) is unbiased no matter what weighting scheme is used 158 Heagerty, 2006
28 Properties of β(w ) (**) 1 When W i is correctly specified as the inverse of the variance of Y i then: W i = Σ 1 i ] var [ β(σ 1 ) = ( i ) 1 X T i Σ 1 i X i When we use gllamm, SAS PROC MIXED, or S+ lme this is what is returned to provide standard errors for the estimated regression coefficients 159 Heagerty, 2006
29 Properties of β(w ) (**) 2 When W i is not the inverse of the variance of Y i then: W i Σ 1 i ] var [ β(w ) = bread ( ) {}}{ A 1 X T i W i Σ i W i X i i }{{} cheese bread {}}{ A 1 A = i X T i W i X i More on this sandwich later 160 Heagerty, 2006
30 Likelihood Estimation for Linear Mixed Models (**) Parameters: β : regression parameter, fixed effects coefficient α : variance components α D(α) and R(α) where cov(y i ) = Z i DZ T i + R i 161 Heagerty, 2006
31 Normality: E(Y i ) = X i β cov(y i ) = Σ(α) f(y i ; β, α) = Σ 1/2 (2π) ni/2 [ exp 1 ] 2 (Y i X i β) T Σ 1 (Y i X i β) Maximum Likelihood: Find the values for the regression coefficients, β, and the variance components that maximizes the likelihood eg put the highest available probability on the observed data 162 Heagerty, 2006
32 RA Fisher Heagerty, 2006
33 ML versus REML There is a variant of ML estimation known as REML Residual ML Restricted ML REML is used to provide slightly less biased estimates of variance components However, be careful using REML when you change the covariates in your model since one can not use changes in REML log likelihoods to test for fixed effects Useful for a single fitted model, or to compare covariance models with a fixed regression model 163 Heagerty, 2006
34 Inference in the Linear Mixed Model In practice: (1) Saturated mean model & explore the covariance (2) Fix the covariance & explore the mean 164 Heagerty, 2006
35 Likelihood Ratio Tests Fixed Effects Standard likelihood theory can be applied to test H 0 : β 2 = 0 where E[Y ] = [X 1, X 2 ] β 1 β 2 = X 1 β 1 + X 2 β 2 [1]Full Model: E[Y ] = X 1 β 1 + X 2 β 2 [0]Reduced Model: E[Y ] = X 1 β Heagerty, 2006
36 Likelihood Ratio Tests Fixed Effects In this case we have (when null hypothesis is true): Likelihood Ratio = L ML( β 1, β 2, α; ML using model 1) L ML ( β 1, 0, α; ML using model 0) LRstatistic = 2 log Likelihood Ratio = 2 log L ML,1 2 log L ML,0 χ 2 (q) Where q is the number of coefficients that are set to zero in the reduced model 166 Heagerty, 2006
37 Other Tests Fixed Effects (*) We also have for a general linear contrast A and a hypothesis H 0 : Aβ = 0 Wald Test: ( ) (A β) T Avar( β)a T 1 (A β) χ 2 (q) F Test: F = (A β) T (A var( β)a T ) 1 (A β) rank(a) F (ndf = rank(a), ddf) 167 Heagerty, 2006
38 LMM: Selection of the Covariance Matrix A Model that fits the data Compare the fitted covariance to the empirical assessment of it: Σ i = Z i DZ T i + R i ( α) versus cov(y i µ i ) γ( ) = τ 2 + σ 2 [1 ρ( )] versus empirical variogram var(y ij ) = τ 2 + σ 2 + ν 2 versus empirical variance 168 Heagerty, 2006
39 Look at the maximized likelihood: Compare 2 log L AIC, BIC Don t lose sight of the goals of analysis If covariance selection is to obtain valid model based standard errors then we can assess the impact on β and se s We can also calculate an empirical (sandwich) variance estimate 169 Heagerty, 2006
40 Inference in the Linear Mixed Model (*) Likelihood Ratio Tests Variance Components We may want to test whether we have random intercepts and slopes, or just random intercepts H 0 : D = D versus H 1 : D = D 11 D 12 D 21 D 22 Q: What is the distribution of the likelihood ratio statistic LRstat = 2 log L ML( θ ML,model 1 ) L ML ( θ ML,model 0 ) 170 Heagerty, 2006
41 LR Testing for Variance Components (*) D 22 = 0 is on the boundary of the parameter space!!! This violates the standard assumption that we use to justify the χ 2 (p 1 p 0 ) distribution of the LR statistic We appeal to results in Stram and Lee (1994) that build upon results in Self & Liang (1987) showing that LR stat is a mixture of χ 2 Note: For a fixed mean strucure we can use the LR based on either ML or REML (Why?) See: Verbeke and Molenberghs (1997) pages Heagerty, 2006
42 171-1 Heagerty, 2006
43 171-2 Heagerty, 2006
44 171-3 Heagerty, 2006 S+ LMM Program: # # cfkids-cda-newlmmq # # # # PURPOSE: Use linear mixed models to characterize longitudinal # change by gender and genotype # # AUTHOR: P Heagerty # # DATE: 00/07/10 Revised 14Feb2002 # # # # ##### ##### Read data ##### # source("cfkids-readq") # # ##### ##### Trellis plots of individuals and groups #####
45 # # Create Grouped Data Set # ntotal <- cumsum( unlist( lapply( split( cfkids$id, cfkids$id), length ) ) ) cfsubset <- groupeddata( fev1 ~ age id, outer = ~ factor(f508)*female, data = cfkids[ 1:ntotal[(8*4*1)], ] ) # cfkids <- groupeddata( fev1 ~ age id, outer = ~ factor(f508)*female, data = cfkids ) # # trellis plot, by id, first 1 pages, 8x4 # postscript( file="cfkids-trellis1ps", horiz=f ) plot( cfsubset, layout = c(4,8) ) graphicsoff() postscript( file="cfkids-trellis2ps", horiz=t ) par( pch="" ) plot( cfkids, outer = ~ factor(f508)*factor(female), layout=c(3,2), aspect=1 ) graphicsoff() # ##### ##### Linear Mixed Models ##### # options( contrasts=c("contrtreatment","contrhelmert") ) Heagerty, 2006
46 171-5 Heagerty, 2006 # ### Intercept only fit0 <- lme( fev1 ~ age0 + agel + female*agel + factor(f508)*agel, method = "ML", random = restruct( ~ 1 id, pdclass="pdsymm", REML=F), data = cfkids ) summary( fit0 ) ### Intercept plus Slope fit1 <- lme( fev1 ~ age0 + agel + female*agel + factor(f508)*agel, method = "ML", random = restruct( ~ 1 + agel id, pdclass="pdsymm", REML=F), data = cfkids ) summary( fit1 ) ### EDA for serial correlation postscript( file="cfkids-variogramps", horiz=t ) plot( Variogram( fit0, form = ~ age id, restype="response" ) ) graphicsoff() ### Intercept plus AR(1) fit2a <- lme( fev1 ~ age0 + agel + female*agel + factor(f508)*agel,
47 171-6 Heagerty, 2006 method = "ML", random = restruct( ~ 1 id, pdclass="pdsymm", REML=F), correlation = corar1( form = ~ 1 id ), data = cfkids ) summary( fit2a ) ### another way fit2b <- lme( fev1 ~ age0 + agel + female*agel + factor(f508)*agel, method = "ML", random = restruct( ~ 1 id, pdclass="pdsymm", REML=F), correlation = corexp( form = ~ agel id, nugget=f), data = cfkids ) summary( fit2b ) ### another way fit2c <- lme( fev1 ~ age0 + agel + female*agel + factor(f508)*agel, method = "ML", random = restruct( ~ 1 id, pdclass="pdsymm", REML=F), correlation = corcar1( form = ~ agel id ), data = cfkids ) summary( fit2c ) fit2 <- fit2b ### Intercept plus AR(1) plus measurement error
48 171-7 Heagerty, 2006 fit3 <- lme( fev1 ~ age0 + agel + female*agel + factor(f508)*agel, method = "ML", random = restruct( ~ 1 id, pdclass="pdsymm", REML=F), correlation = corexp( form = ~ agel id, nugget=t), data = cfkids ) summary( fit3 ) # ##### compare these models # anova( fit0, fit1, fit2, fit3 ) # ##### ##### Residual Analysis -- using fit3 ##### # popres <- resid( fit3, level=0 ) clusterres <- resid( fit3, level=1 ) print( var( popres ) ) print( var( clusterres ) ) # postscript( file="cfkids-newresidualsps", horiz=f ) par( mfrow=c(2,1) ) plot( cfkids$age0, popres, pch="" ) lines( smoothspline( cfkids$age0, popres, df=5 ) ) title("residuals (pop) vs Age0") abline( h=0, lty=2 )
49 171-8 Heagerty, 2006 plot( cfkids$agel, popres, pch="" ) lines( smoothspline( cfkids$agel, popres, df=5 ) ) title("residuals (pop) vs AgeL") abline( h=0, lty=2 ) graphicsoff() # postscript( file="cfkids-newresiduals2ps", horiz=f ) par( mfrow=c(2,1) ) plot( cfkids$age0, clusterres, pch="" ) lines( smoothspline( cfkids$age0, clusterres, df=5 ) ) abline( h=0, lty=2 ) title("residuals (cluster) vs Age0") b0 <- unlist( fit2$coefficients$random ) age0 <- unlist( lapply( split( cfkids$age0, cfkids$id ), min ) ) plot( age0, b0 ) lines( smoothspline( age0, b0, df=5 ) ) abline( h=0, lty=2 ) title("eb b0 versus Age0") graphicsoff() # ##### Do we need a quadratic age0??? # fit4 <- lme( fev1 ~ age0 + age0^2 + agel + female*agel + factor(f508)*agel, method = "ML", random = restruct( ~ 1 id, pdclass="pdsymm", REML=F), correlation = corexp( form = ~ agel id, nugget=t), data = cfkids )
50 summary( fit4 ) # anova( fit3, fit4 ) # # end-of-file Heagerty, 2006
51 Heagerty, fev age
52 fev age Heagerty, 2006
53 Heagerty, 2006 Fit 0 Random Intercepts Linear mixed-effects model fit by maximum likelihood Data: cfkids AIC BIC loglik Random effects: Formula: ~ 1 id (Intercept) Residual StdDev: Fixed effects: fev1 ~ age0 + agel + female * agel + factor(f508) * agel Value StdError DF t-value p-value (Intercept) <0001 age <0001 agel female factor(f508) factor(f508) female:agel agelfactor(f508) agelfactor(f508) Number of Observations: 1513 Number of Groups: 200
54 Heagerty, Semivariogram Distance
55 Heagerty, 2006 Fit 1 Random Intercepts and Slopes Linear mixed-effects model fit by maximum likelihood Data: cfkids AIC BIC loglik Random effects: StdDev Corr Formula: ~ 1 + agel id (Intercept) (Inter Structure: General positive-definite agel Residual Fixed effects: fev1 ~ age0 + agel + female * agel + factor(f508) * agel Value StdError DF t-value p-value (Intercept) <0001 age <0001 agel female factor(f508) factor(f508) female:agel agelfactor(f508) agelfactor(f508) Number of Observations: 1513 Number of Groups: 200
56 Heagerty, 2006 Fit 2a Random Intercepts + AR(1) errors Linear mixed-effects model fit by maximum likelihood Data: cfkids AIC BIC loglik Random effects: Correlation Structure: AR(1) Formula: ~ 1 id Formula: ~ 1 id (Intercept) Residual Parameter estimate(s): StdDev: Phi Fixed effects: fev1 ~ age0 + agel + female * agel + factor(f508) * agel Value StdError DF t-value p-value (Intercept) <0001 age <0001 agel female factor(f508) factor(f508) female:agel agelfactor(f508) agelfactor(f508) Number of Observations: 1513 Number of Groups: 200
57 Heagerty, 2006 Fit 2b Random Intercepts + corexp errors Linear mixed-effects model fit by maximum likelihood Data: cfkids AIC BIC loglik Random effects: Correlation Structure: Exponential spatial corr Formula: ~ 1 id Formula: ~ agel id (Intercept) Residual Parameter estimate(s): StdDev: range Fixed effects: fev1 ~ age0 + agel + female * agel + factor(f508) * agel Value StdError DF t-value p-value (Intercept) <0001 age <0001 agel female factor(f508) factor(f508) female:agel agelfactor(f508) agelfactor(f508) Number of Observations: 1513 Number of Groups: 200
58 Heagerty, 2006 Fit 2c Random Intercepts + CAR(1) errors Linear mixed-effects model fit by maximum likelihood Data: cfkids AIC BIC loglik Random effects: Correlation Structure: Continuous AR(1) Formula: ~ 1 id Formula: ~ agel id (Intercept) Residual Parameter estimate(s): StdDev: Phi Fixed effects: fev1 ~ age0 + agel + female * agel + factor(f508) * agel Value StdError DF t-value p-value (Intercept) <0001 age <0001 agel female factor(f508) factor(f508) female:agel agelfactor(f508) agelfactor(f508) Number of Observations: 1513 Number of Groups: 200
59 Heagerty, 2006 Fit 3 Random Intercepts + corexp + meas error Linear mixed-effects model fit by maximum likelihood Data: cfkids AIC BIC loglik Random effects: Correlation Structure: Exponential spatial corr Formula: ~ 1 id Formula: ~ agel id (Intercept) Residual Parameter estimate(s): StdDev: range nugget Fixed effects: fev1 ~ age0 + agel + female * agel + factor(f508) * agel Value StdError DF t-value p-value (Intercept) <0001 age <0001 agel female factor(f508) factor(f508) female:agel agelfactor(f508) agelfactor(f508) Number of Observations: 1513 Number of Groups: 200
60 Heagerty, 2006 ANOVA Model df AIC BIC loglik Test LRatio p-value fit fit vs <0001 fit vs e-04 fit vs <0001
61 LMM Summary Observe Y i, i = 1, 2,, m independent clusters Model: (Laird & Ware, 1982) Y i = X i β + Z i b i + e i β is the coefficient that is common to all clusters (fixed across clusters) b i is the deviation of the coefficient that varies from cluster to cluster (random across clusters) (β j + b j,i ) is the coefficient of X i,j for cluster i 172 Heagerty, 2006
62 b i N (0, D) between-cluster e i N (0, R i ) within-cluster Estimation/Inference: WLS, ML Covariance model choice leads to WLS but estimated regression coefficient is unbiased for any choice of weight (covariance) Covariance model choice determines the standard error estimates for the regression coefficients correct covariance model is needed for correct standard errors 173 Heagerty, 2006
General Linear Model for Correlated Data
General Linear Model for Correlated Data Objectives: Weighted least squares methods Moment estimation Sandwich variance Linear mixed models. Models for covariance Maximum likelihood and REML Empirical
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 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 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 informationPart III: Linear mixed effects models
Part III: Linear mixed effects models 159 BIO 245, Spring 2018 Dental growth data Females Males Length, mm 16 18 20 22 24 26 28 30 32 11 11 4 3 2 11 4 7 11 4 3 3 2 8 4 8 5 8 58 7 1 7 5 6 9 57 2 1 9 12
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 informationRepeated measures, part 2, advanced methods
enote 12 1 enote 12 Repeated measures, part 2, advanced methods enote 12 INDHOLD 2 Indhold 12 Repeated measures, part 2, advanced methods 1 12.1 Intro......................................... 3 12.2 A
More informationRestricted Maximum Likelihood in Linear Regression and Linear Mixed-Effects Model
Restricted Maximum Likelihood in Linear Regression and Linear Mixed-Effects Model Xiuming Zhang zhangxiuming@u.nus.edu A*STAR-NUS Clinical Imaging Research Center October, 015 Summary This report derives
More informationIntroduction to General and Generalized Linear Models
Introduction to General and Generalized Linear Models Mixed effects models - II Henrik Madsen, Jan Kloppenborg Møller, Anders Nielsen April 16, 2012 H. Madsen, JK. Møller, A. Nielsen () Chapman & Hall
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 informationStat 579: Generalized Linear Models and Extensions
Stat 579: Generalized Linear Models and Extensions Linear Mixed Models for Longitudinal Data Yan Lu April, 2018, week 14 1 / 64 Data structure and Model t1 t2 tn i 1st subject y 11 y 12 y 1n1 2nd subject
More information13. October p. 1
Lecture 8 STK3100/4100 Linear mixed models 13. October 2014 Plan for lecture: 1. The lme function in the nlme library 2. Induced correlation structure 3. Marginal models 4. Estimation - ML and REML 5.
More 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 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 informationThese slides illustrate a few example R commands that can be useful for the analysis of repeated measures data.
These slides illustrate a few example R commands that can be useful for the analysis of repeated measures data. We focus on the experiment designed to compare the effectiveness of three strength training
More informationIntroduction and Background to Multilevel Analysis
Introduction and Background to Multilevel Analysis Dr. J. Kyle Roberts Southern Methodist University Simmons School of Education and Human Development Department of Teaching and Learning Background and
More informationHierarchical Linear Models (HLM) Using R Package nlme. Interpretation. 2 = ( x 2) u 0j. e ij
Hierarchical Linear Models (HLM) Using R Package nlme Interpretation I. The Null Model Level 1 (student level) model is mathach ij = β 0j + e ij Level 2 (school level) model is β 0j = γ 00 + u 0j Combined
More informationIntruction to General and Generalized Linear Models
Intruction to General and Generalized Linear Models Mixed Effects Models IV Henrik Madsen Anna Helga Jónsdóttir hm@imm.dtu.dk April 30, 2012 Henrik Madsen Anna Helga Jónsdóttir (hm@imm.dtu.dk) Intruction
More informationIntroduction to the Analysis of Hierarchical and Longitudinal Data
Introduction to the Analysis of Hierarchical and Longitudinal Data Georges Monette, York University with Ye Sun SPIDA June 7, 2004 1 Graphical overview of selected concepts Nature of hierarchical models
More 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 informationMixed models with correlated measurement errors
Mixed models with correlated measurement errors Rasmus Waagepetersen October 9, 2018 Example from Department of Health Technology 25 subjects where exposed to electric pulses of 11 different durations
More 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 informationModeling the Covariance
Modeling the Covariance Jamie Monogan University of Georgia February 3, 2016 Jamie Monogan (UGA) Modeling the Covariance February 3, 2016 1 / 16 Objectives By the end of this meeting, participants should
More informationLongitudinal Data Analysis
Longitudinal Data Analysis Mike Allerhand This document has been produced for the CCACE short course: Longitudinal Data Analysis. No part of this document may be reproduced, in any form or by any means,
More informationThe consequences of misspecifying the random effects distribution when fitting generalized linear mixed models
The consequences of misspecifying the random effects distribution when fitting generalized linear mixed models John M. Neuhaus Charles E. McCulloch Division of Biostatistics University of California, San
More informationAnalysis of Longitudinal Data. Patrick J. Heagerty PhD Department of Biostatistics University of Washington
Analysis of Longitudinal Data Patrick J Heagerty PhD Department of Biostatistics University of Washington Auckland 8 Session One Outline Examples of longitudinal data Scientific motivation Opportunities
More informationIntroduction to Linear Mixed Models: Modeling continuous longitudinal outcomes
1/64 to : Modeling continuous longitudinal outcomes Dr Cameron Hurst cphurst@gmail.com CEU, ACRO and DAMASAC, Khon Kaen University 4 th Febuary, 2557 2/64 Some motivational datasets Before we start, I
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 informationLecture 2: Linear Models. Bruce Walsh lecture notes Seattle SISG -Mixed Model Course version 23 June 2011
Lecture 2: Linear Models Bruce Walsh lecture notes Seattle SISG -Mixed Model Course version 23 June 2011 1 Quick Review of the Major Points The general linear model can be written as y = X! + e y = vector
More 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 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 informationLISA Short Course Series Generalized Linear Models (GLMs) & Categorical Data Analysis (CDA) in R. Liang (Sally) Shan Nov. 4, 2014
LISA Short Course Series Generalized Linear Models (GLMs) & Categorical Data Analysis (CDA) in R Liang (Sally) Shan Nov. 4, 2014 L Laboratory for Interdisciplinary Statistical Analysis LISA helps VT researchers
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 informationNon-independence due to Time Correlation (Chapter 14)
Non-independence due to Time Correlation (Chapter 14) When we model the mean structure with ordinary least squares, the mean structure explains the general trends in the data with respect to our dependent
More informationThe linear mixed model: introduction and the basic model
The linear mixed model: introduction and the basic model Analysis of Experimental Data AED The linear mixed model: introduction and the basic model of 39 Contents Data with non-independent observations
More informationLecture 3: Linear Models. Bruce Walsh lecture notes Uppsala EQG course version 28 Jan 2012
Lecture 3: Linear Models Bruce Walsh lecture notes Uppsala EQG course version 28 Jan 2012 1 Quick Review of the Major Points The general linear model can be written as y = X! + e y = vector of observed
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 informationCorrelation and the Analysis of Variance Approach to Simple Linear Regression
Correlation and the Analysis of Variance Approach to Simple Linear Regression Biometry 755 Spring 2009 Correlation and the Analysis of Variance Approach to Simple Linear Regression p. 1/35 Correlation
More informationMath 423/533: The Main Theoretical Topics
Math 423/533: The Main Theoretical Topics Notation sample size n, data index i number of predictors, p (p = 2 for simple linear regression) y i : response for individual i x i = (x i1,..., x ip ) (1 p)
More information22s:152 Applied Linear Regression. Returning to a continuous response variable Y...
22s:152 Applied Linear Regression Generalized Least Squares Returning to a continuous response variable Y... Ordinary Least Squares Estimation The classical models we have fit so far with a continuous
More information22s:152 Applied Linear Regression. In matrix notation, we can write this model: Generalized Least Squares. Y = Xβ + ɛ with ɛ N n (0, Σ)
22s:152 Applied Linear Regression Generalized Least Squares Returning to a continuous response variable Y Ordinary Least Squares Estimation The classical models we have fit so far with a continuous response
More informationFinal Review. Yang Feng. Yang Feng (Columbia University) Final Review 1 / 58
Final Review Yang Feng http://www.stat.columbia.edu/~yangfeng Yang Feng (Columbia University) Final Review 1 / 58 Outline 1 Multiple Linear Regression (Estimation, Inference) 2 Special Topics for Multiple
More informationMatrix Approach to Simple Linear Regression: An Overview
Matrix Approach to Simple Linear Regression: An Overview Aspects of matrices that you should know: Definition of a matrix Addition/subtraction/multiplication of matrices Symmetric/diagonal/identity matrix
More informationAn R # Statistic for Fixed Effects in the Linear Mixed Model and Extension to the GLMM
An R Statistic for Fixed Effects in the Linear Mixed Model and Extension to the GLMM Lloyd J. Edwards, Ph.D. UNC-CH Department of Biostatistics email: Lloyd_Edwards@unc.edu Presented to the Department
More informationMaximum Likelihood (ML) Estimation
Econometrics 2 Fall 2004 Maximum Likelihood (ML) Estimation Heino Bohn Nielsen 1of32 Outline of the Lecture (1) Introduction. (2) ML estimation defined. (3) ExampleI:Binomialtrials. (4) Example II: Linear
More informationSTAT 526 Advanced Statistical Methodology
STAT 526 Advanced Statistical Methodology Fall 2017 Lecture Note 10 Analyzing Clustered/Repeated Categorical Data 0-0 Outline Clustered/Repeated Categorical Data Generalized Linear Mixed Models Generalized
More informationA brief introduction to mixed models
A brief introduction to mixed models University of Gothenburg Gothenburg April 6, 2017 Outline An introduction to mixed models based on a few examples: Definition of standard mixed models. Parameter estimation.
More informationIntroduction to 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 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 informationFirst Year Examination Department of Statistics, University of Florida
First Year Examination Department of Statistics, University of Florida August 20, 2009, 8:00 am - 2:00 noon Instructions:. You have four hours to answer questions in this examination. 2. You must show
More informationStat/F&W Ecol/Hort 572 Review Points Ané, Spring 2010
1 Linear models Y = Xβ + ɛ with ɛ N (0, σ 2 e) or Y N (Xβ, σ 2 e) where the model matrix X contains the information on predictors and β includes all coefficients (intercept, slope(s) etc.). 1. Number of
More informationStat 209 Lab: Linear Mixed Models in R This lab covers the Linear Mixed Models tutorial by John Fox. Lab prepared by Karen Kapur. ɛ i Normal(0, σ 2 )
Lab 2 STAT209 1/31/13 A complication in doing all this is that the package nlme (lme) is supplanted by the new and improved lme4 (lmer); both are widely used so I try to do both tracks in separate Rogosa
More informationOutline. Mixed models in R using the lme4 package Part 3: Longitudinal data. Sleep deprivation data. Simple longitudinal data
Outline Mixed models in R using the lme4 package Part 3: Longitudinal data Douglas Bates Longitudinal data: sleepstudy A model with random effects for intercept and slope University of Wisconsin - Madison
More informationAdvantages of Mixed-effects Regression Models (MRM; aka multilevel, hierarchical linear, linear mixed models) 1. MRM explicitly models individual
Advantages of Mixed-effects Regression Models (MRM; aka multilevel, hierarchical linear, linear mixed models) 1. MRM explicitly models individual change across time 2. MRM more flexible in terms of repeated
More 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 informationLecture 3 Linear random intercept models
Lecture 3 Linear random intercept models Example: Weight of Guinea Pigs Body weights of 48 pigs in 9 successive weeks of follow-up (Table 3.1 DLZ) The response is measures at n different times, or under
More informationUsing Estimating Equations for Spatially Correlated A
Using Estimating Equations for Spatially Correlated Areal Data December 8, 2009 Introduction GEEs Spatial Estimating Equations Implementation Simulation Conclusion Typical Problem Assess the relationship
More information20. REML Estimation of Variance Components. Copyright c 2018 (Iowa State University) 20. Statistics / 36
20. REML Estimation of Variance Components Copyright c 2018 (Iowa State University) 20. Statistics 510 1 / 36 Consider the General Linear Model y = Xβ + ɛ, where ɛ N(0, Σ) and Σ is an n n positive definite
More information36-463/663: Hierarchical Linear Models
36-463/663: Hierarchical Linear Models Lmer model selection and residuals Brian Junker 132E Baker Hall brian@stat.cmu.edu 1 Outline The London Schools Data (again!) A nice random-intercepts, random-slopes
More informationMixed effects models
Mixed effects models The basic theory and application in R Mitchel van Loon Research Paper Business Analytics Mixed effects models The basic theory and application in R Author: Mitchel van Loon Research
More 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 informationSTAT 100C: Linear models
STAT 100C: Linear models Arash A. Amini June 9, 2018 1 / 56 Table of Contents Multiple linear regression Linear model setup Estimation of β Geometric interpretation Estimation of σ 2 Hat matrix Gram matrix
More information[y i α βx i ] 2 (2) Q = i=1
Least squares fits This section has no probability in it. There are no random variables. We are given n points (x i, y i ) and want to find the equation of the line that best fits them. We take the equation
More informationLecture 9 STK3100/4100
Lecture 9 STK3100/4100 27. October 2014 Plan for lecture: 1. Linear mixed models cont. Models accounting for time dependencies (Ch. 6.1) 2. Generalized linear mixed models (GLMM, Ch. 13.1-13.3) Examples
More informationModeling the scale parameter ϕ A note on modeling correlation of binary responses Using marginal odds ratios to model association for binary responses
Outline Marginal model Examples of marginal model GEE1 Augmented GEE GEE1.5 GEE2 Modeling the scale parameter ϕ A note on modeling correlation of binary responses Using marginal odds ratios to model association
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 informationModel Estimation Example
Ronald H. Heck 1 EDEP 606: Multivariate Methods (S2013) April 7, 2013 Model Estimation Example As we have moved through the course this semester, we have encountered the concept of model estimation. Discussions
More informationUNIVERSITY OF TORONTO Faculty of Arts and Science
UNIVERSITY OF TORONTO Faculty of Arts and Science December 2013 Final Examination STA442H1F/2101HF Methods of Applied Statistics Jerry Brunner Duration - 3 hours Aids: Calculator Model(s): Any calculator
More informationStatistical Distribution Assumptions of General Linear Models
Statistical Distribution Assumptions of General Linear Models Applied Multilevel Models for Cross Sectional Data Lecture 4 ICPSR Summer Workshop University of Colorado Boulder Lecture 4: Statistical Distributions
More informationMixed models in R using the lme4 package Part 2: Longitudinal data, modeling interactions
Mixed models in R using the lme4 package Part 2: Longitudinal data, modeling interactions Douglas Bates Department of Statistics University of Wisconsin - Madison Madison January 11, 2011
More informationMixed models in R using the lme4 package Part 5: Generalized linear mixed models
Mixed models in R using the lme4 package Part 5: Generalized linear mixed models Douglas Bates 2011-03-16 Contents 1 Generalized Linear Mixed Models Generalized Linear Mixed Models When using linear mixed
More informationLecture 3.1 Basic Logistic LDA
y Lecture.1 Basic Logistic LDA 0.2.4.6.8 1 Outline Quick Refresher on Ordinary Logistic Regression and Stata Women s employment example Cross-Over Trial LDA Example -100-50 0 50 100 -- Longitudinal Data
More informationLab 3: Two levels Poisson models (taken from Multilevel and Longitudinal Modeling Using Stata, p )
Lab 3: Two levels Poisson models (taken from Multilevel and Longitudinal Modeling Using Stata, p. 376-390) BIO656 2009 Goal: To see if a major health-care reform which took place in 1997 in Germany was
More informationMixed models in R using the lme4 package Part 5: Generalized linear mixed models
Mixed models in R using the lme4 package Part 5: Generalized linear mixed models Douglas Bates Madison January 11, 2011 Contents 1 Definition 1 2 Links 2 3 Example 7 4 Model building 9 5 Conclusions 14
More informationDe-mystifying random effects models
De-mystifying random effects models Peter J Diggle Lecture 4, Leahurst, October 2012 Linear regression input variable x factor, covariate, explanatory variable,... output variable y response, end-point,
More informationIntroduction to Hierarchical Data Theory Real Example. NLME package in R. Jiang Qi. Department of Statistics Renmin University of China.
Department of Statistics Renmin University of China June 7, 2010 The problem Grouped data, or Hierarchical data: correlations between subunits within subjects. The problem Grouped data, or Hierarchical
More informationTime-Series Regression and Generalized Least Squares in R*
Time-Series Regression and Generalized Least Squares in R* An Appendix to An R Companion to Applied Regression, third edition John Fox & Sanford Weisberg last revision: 2018-09-26 Abstract Generalized
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 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 equivalence of the Maximum Likelihood and a modified Least Squares for a case of Generalized Linear Model
Applied and Computational Mathematics 2014; 3(5): 268-272 Published online November 10, 2014 (http://www.sciencepublishinggroup.com/j/acm) doi: 10.11648/j.acm.20140305.22 ISSN: 2328-5605 (Print); ISSN:
More informationLatent Variable Models for Binary Data. Suppose that for a given vector of explanatory variables x, the latent
Latent Variable Models for Binary Data Suppose that for a given vector of explanatory variables x, the latent variable, U, has a continuous cumulative distribution function F (u; x) and that the binary
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 informationRegression Models - Introduction
Regression Models - Introduction In regression models there are two types of variables that are studied: A dependent variable, Y, also called response variable. It is modeled as random. An independent
More informationSTAT Financial Time Series
STAT 6104 - Financial Time Series Chapter 4 - Estimation in the time Domain Chun Yip Yau (CUHK) STAT 6104:Financial Time Series 1 / 46 Agenda 1 Introduction 2 Moment Estimates 3 Autoregressive Models (AR
More informationGeneralized Linear. Mixed Models. Methods and Applications. Modern Concepts, Walter W. Stroup. Texts in Statistical Science.
Texts in Statistical Science Generalized Linear Mixed Models Modern Concepts, Methods and Applications Walter W. Stroup CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint
More informationVector autoregressions, VAR
1 / 45 Vector autoregressions, VAR Chapter 2 Financial Econometrics Michael Hauser WS17/18 2 / 45 Content Cross-correlations VAR model in standard/reduced form Properties of VAR(1), VAR(p) Structural VAR,
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More informationOutline. Statistical inference for linear mixed models. One-way ANOVA in matrix-vector form
Outline Statistical inference for linear mixed models Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark general form of linear mixed models examples of analyses using linear mixed
More informationModeling the Mean: Response Profiles v. Parametric Curves
Modeling the Mean: Response Profiles v. Parametric Curves Jamie Monogan University of Georgia Escuela de Invierno en Métodos y Análisis de Datos Universidad Católica del Uruguay Jamie Monogan (UGA) Modeling
More informationMultivariate Regression
Multivariate Regression The so-called supervised learning problem is the following: we want to approximate the random variable Y with an appropriate function of the random variables X 1,..., X p with the
More informationPAPER 206 APPLIED STATISTICS
MATHEMATICAL TRIPOS Part III Thursday, 1 June, 2017 9:00 am to 12:00 pm PAPER 206 APPLIED STATISTICS Attempt no more than FOUR questions. There are SIX questions in total. The questions carry equal weight.
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More informationRestricted Likelihood Ratio Tests in Nonparametric Longitudinal Models
Restricted Likelihood Ratio Tests in Nonparametric Longitudinal Models Short title: Restricted LR Tests in Longitudinal Models Ciprian M. Crainiceanu David Ruppert May 5, 2004 Abstract We assume that repeated
More informationOrdinary Least Squares Regression
Ordinary Least Squares Regression Goals for this unit More on notation and terminology OLS scalar versus matrix derivation Some Preliminaries In this class we will be learning to analyze Cross Section
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 informationA Multivariate Two-Sample Mean Test for Small Sample Size and Missing Data
A Multivariate Two-Sample Mean Test for Small Sample Size and Missing Data Yujun Wu, Marc G. Genton, 1 and Leonard A. Stefanski 2 Department of Biostatistics, School of Public Health, University of Medicine
More informationMixed models in R using the lme4 package Part 7: Generalized linear mixed models
Mixed models in R using the lme4 package Part 7: Generalized linear mixed models Douglas Bates University of Wisconsin - Madison and R Development Core Team University of
More informationModels for Clustered Data
Models for Clustered Data Edps/Psych/Soc 589 Carolyn J Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Spring 2019 Outline Notation NELS88 data Fixed Effects ANOVA
More informationSTATISTICS 174: APPLIED STATISTICS FINAL EXAM DECEMBER 10, 2002
Time allowed: 3 HOURS. STATISTICS 174: APPLIED STATISTICS FINAL EXAM DECEMBER 10, 2002 This is an open book exam: all course notes and the text are allowed, and you are expected to use your own calculator.
More informationModels for Clustered Data
Models for Clustered Data Edps/Psych/Stat 587 Carolyn J Anderson Department of Educational Psychology c Board of Trustees, University of Illinois Fall 2017 Outline Notation NELS88 data Fixed Effects ANOVA
More information