Estimating a Piecewise Growth Model with Longitudinal Data that Contains Individual Mobility across Clusters
|
|
- Oswin Walker
- 5 years ago
- Views:
Transcription
1 Estimating a Piecewise Growth Model with Longitudinal Data that Contains Individual Mobility across Clusters Audrey J. Leroux Georgia State University
2 Piecewise Growth Model (PGM) PGMs are beneficial for potentially nonlinear data, because they break up curvilinear growth trajectories into separate linear components This modeling approach is useful when wanting to compare growth rates during two or more different time periods. For example: Compare growth rates during two or more different time periods Longitudinal data before treatment as well as during treatment Longitudinal data during treatment as well as follow-up data after treatment Etc.
3 Three-Level PGM Three-level PGMs are used to model the clustering of individuals, such as when students are nested within schools, classrooms, districts, etc. in educational research. 3
4 Baseline Three-Level PGM Level 1: Y tij = π 0ij + π 1ij Time 1tij + π ij Time tij + e tij, e tij ~ N 0, σ e where Time 1tij and Time tij are coded to represent piecewise growth Of course, can include predictors Level : Level 3: π 0ij = β 00j + r 0ij π 1ij = β 10j + r 1ij π ij = β 0j + r ij, β 00j = γ u 00j β 10j = γ u 10j, β 0j = γ 00 + u 0j u 00j u 10j u 0j r 0ij r 1ij r ij ~ MVN ~ MVN 0 0 0, 0 0 0, τ u00 τ u00,u10 τ r0 τ r0,r1 τ r1 τ r0,r τ r1,r τ r τ u10 τ u00,u0 τ u10,u0 τ u0 4
5 Three-Level PGM To estimate this particular model presented with initial status and two slopes varying across both individuals and clusters, a minimum of four time-points must be included for identification of the model. One more time-point than the number of growth parameters (three). Otherwise, modeling the slopes as fixed or constraints on the level-1 residual variance can be placed if utilizing fewer time-points (see McCoach, O'Connell, Reis, & Levitt, 006; Palardy, 010). 5
6 Three-Level PGM Three-level PGMs utilized in previous research have assumed pure clustering of individuals across time or removed individuals from the analysis who changed clusters. BUT in reality, individual mobility across clusters is frequently encountered in longitudinal studies. Incorrect model specification in the presence of cluster mobility negatively impacts parameter estimates (Chung & Beretvas, 01; Grady, 010; Grady & Beretvas, 010; Leroux, 014; Leroux & Beretvas, 017a, Leroux & Beretvas, 017b; Luo & Kwok, 009; Luo & Kwok, 01; Meyers & Beretvas, 006). Generally, leads to inaccurate estimates of between-clusters variance components and standard errors of the fixed effects. 6
7 Longitudinal Data with Mobile Students Student Fall K Spring K Spring 1 st Spring 3 rd Spring 5 th Sch. 1 Sch. 1 Sch. Sch. 1 Sch. Sch. 3 Sch. 1 Sch. Sch. 3 Sch. 4 Sch. 1 Sch. Sch. 3 Sch. 4 Sch. 5 A B C D E 7
8 Multiple Membership Data Elementary ES1 ES ES3 ES4 ES5 ES6 ES7 ES8 School Student A B C D E F G H I J K L M N O P Q R S T U V Some units of a lower-level classification are members of more than one higherlevel classification. 8
9 Multiple Membership Random Effects Model (MMREM) Models the contribution to the outcome, Y, of each level- unit of which the level- 1 unit is a member. E.g., For student i who is a member of a set of ESs {j}, the unconditional model s L1 equation is: Y i j = β 0 j + r i j At L: β 0 j = γ 00 + h j w ih u 0h 9
10 MMREM Single equation: Y i j = γ 00 + h j w ih u 0h + r i j r i j ~ N 0, σ and u 0h ~ N 0, τ u00 where the user specifies the weights to represent the hypothesized contribution of each L unit (here, elementary school) For each student i: h j w ih = 1 10
11 MMREM For non-mobile student B (attending ES1): Y B ES1 = γ 00 + u 0 ES1 + r B ES1 For mobile student A, attending ES1 and ES: Y A ES1,ES = γ u 0 ES u 0 ES + r A ES1,ES For mobile student Q, attending ES6, ES7, and ES8: Y Q ES6,ES7,ES8 = γ u 0 ES u 0 ES u 0 ES8 + r Q ES6,ES7,ES8 11
12 Conditional MMREM Can include L1 and L predictors. At L1: And at L: Y i j = β 0 j + β 1 j X i j + r i j β 0 j = γ 00 + β 1 j = γ 10 + h j h j w ih γ 01 Z h + u 0h w ih γ 11 Z h + u 1h 1
13 Conditional MMREM The following multivariate normal distribution is assumed for the level- residuals: u 0 j u 1 j ~ N 0 0, τ u00 τ u10 τ u11 Note: Contribution of each ES s Z (for mobile students) is modeled and weighted in the same way as are schools effects (the u s). 13
14 Conditional MMREM Private (coded 1) Public (coded 0) For mobile student A, attending private ES1 and public ES: Y A ES1,ES = γ 00 + γ u 0 ES u 0 ES + r A ES1,ES 14
15 Cross-Classified Multiple Membership Longitudinal Data School in 1ES1 1ES 1ES3 1ES4 Time 1 Student A B C D E F G H I J K L M N O P Q R S T U V School after +ES1 +ES +ES3 +ES4 Time 1 Grady & Beretvas (010) 15
16 Cross-Classified Multiple Membership Longitudinal Data School in 1ES1 1ES 1ES3 1ES4 Time 1 Student A B C D E F G H I J K L M N O P Q R S T U V School after +ES1 +ES +ES3 +ES4 Time 1 Grady & Beretvas (010) 16
17 Purpose of Current Study The current study proposes a three-level PGM to handle mobile students who change schools (clusters) during the period of data collection. The proposed cross-classified multiple membership PGM (CCMM-PGM) will be derived, justified, and explained using a real dataset. 17
18 Purpose of Current Study This extension is of particular importance when repeated measures over time are captured for students within schools because there is a high probability that at least some substantial proportion of students change schools during a study s time period. 38.5% of people aged 5-17 years moved within 005 to 010 (Ihrke & Faber, 01) 5% of those between the ages 5-17 relocated within the same county From 01 to 013, 1% of people between the ages 5-17 years old moved 69% of those moves occurred within the same county (U.S. Census Bureau, 013) 13% of students changed schools 4 or more times between kindergarten and 8th grade (U.S. Government accounting office, 010) 18
19 Baseline CCMM-PGM Level 1: Y ti j1, j = π 0i j1, j + π 1i j1, j TIME 1ti j 1, j + π i j1, j TIME ti j 1, j + e ti j1, j Level : Intercept (initial status) π 0i j1, j = β 00 j1, j + r 0i j1, j π 1i j1, j = β 10 j1, j + r 1i j1, j π i j1, j = β 0 j1, j + r i j1, j β 00 j1, j = γ u 00j1 0 Subscripts j 1 and {j } index the first and set of subsequent schools attended by a student. No subsequent school residual for initial status Level 3: 1 st Slope β 10 j1, j = γ u 10j1 0 + w 1tih u 100h h j β 0 j1, j = γ u 0j1 0 + w tih u 00h h j Note two different weights because each slope is associated with different time-points nd Slope Cross-classification of first and subsequent schools 19
20 Data ECLS-K data were used with time nested within students nested within schools. Multiple membership structure due to some students switching elementary schools across the course of data collection Time-points: Fall of kindergarten and springs of kindergarten, 1 st, 3 rd, and 5 th grade Outcome: Math IRT-scaled scores Growth rates from K 1 st grade appeared faster than those from 1 st 5 th grade Gender (1 = female; 0 = male) and school type (1 = private; 0 = public) 10,906 students (9.8% mobile) from 970 schools 0
21 Descriptive Statistics Variable Name M SD N Outcome Math achievement in Fall Kindergarten Y 1ij ,74 Math achievement in Spring Kindergarten Y ij ,664 Math achievement in Spring 1 st Grade Y 3ij ,803 Math achievement in Spring 3 rd Grade Y 4ij ,764 Math achievement in Spring 5 th Grade Y 5ij ,801 Variable Name Percentage N Level- variable Female student 49.78% 5,49 FEMALE Male student ij 50.% 5,477 Level-3 variable Private school 3.51% 8 PRIVATE Public school j 76.49% 74 1
22 Coding of Time Variables Grades Fall K Spring K Spring 1 st Spring 3 rd Spring 5 th Time 1tij Interpretation of πs π 0ij status in Fall K π 1ij growth rate period 1 Time tij π ij growth rate period Exploratory analyses suggested a two-piece growth model because growth rates from kindergarten through 1 st grade appeared faster than those from 1 st through 5 th grade
23 Analyses Baseline and conditional versions of the following models were estimated: CCMM-PGM: appropriately took into account student mobility First school-pgm: ignored mobility by only modeling effect of the first school attended Delete-PGM: ignored mobility by deleting students who changed schools Weights are based on the proportion of time-points a student was associated with a school. 3
24 Coding Schemes for Weights Student Fall K Spring K Schools Weights (K 1 st ) Weights (1 st 5 th ) Spring Spring Spring 1 st 1 st 3 rd 5 th School nd School 1 st School nd School 3 rd School A S1 S1 S1 S1 S B S1 S1 S1 S1 S 1 0 3/4 1/4 0 0 C S1 S1 S1 S S 1 0 1/ 1/ 0 0 D S1 S1 S S S 1/ 1/ 1/4 3/4 0 0 E S1 S S S S F S1 S S S S /4 1/4 0 0 G S1 S S S3 S / 1/ 0 0 H S1 S S3 S3 S3 1/ 1/ 1/4 3/4 0 0 I S1 S S3 S3 S4 1/ 1/ 1/4 1/ 1/4 0 J S1 S S3 S4 S4 1/ 1/ 1/4 1/4 1/ 0 4 th School K S1 S S3 S4 S5 1/ 1/ 1/4 1/4 1/4 1/4 4
25 Estimation Models were fit using R with MCMC estimation using Rjags to interface with Just Another Gibbs Sampler (JAGS). Non-informative normal priors were used for fixed effects parameters and inverse-wishart distributions for the covariance matrices. Burn-in period of 5,000 iterations and an additional 50,000 iterations Parameter and SE estimates were compared, as well as model fit using the deviance information criterion value (DIC; Spiegelhalter, Best, Carlin, & van der Linde, 00). 5
26 Baseline Fixed Effects Estimating Model CCMM-PGM 1 First School-PGM 1 Delete-PGM Parameter Coeff. Est. (SE) Coeff. Est. (SE) Coeff. Est. (SE) Model for initial status Intercept γ (0.171) γ (0.171) γ (0.194) Model for 1 st slope Intercept γ (0.150) γ (0.146) γ (0.171) Model for nd slope Intercept γ (0.066) γ (0.066) γ (0.078) DIC 451, , ,
27 Baseline Random Effects Estimating Model CCMM-PGM First School-PGM Delete-PGM Parameter Coeff. Est. (SE) Coeff. Est. (SE) Coeff. Est. (SE) Level-1 variance between Measures σ (0.631) σ (0.66) σ (0.706) Initial status variance between Students 1 st schools τ r0 τ u0j (1.057) (1.69) τ r0 τ u (1.00) (1.88) τ r0 τ u (1.73) (1.495) 1 st slope variance between Students 1 st schools Subsequent schools τ r1 τ u1j1 τ u1 j (1.13) (1.139) (0.990) τ r1 τ u (1.113) (0.97) τ r1 τ u (1.79) (1.13) nd slope variance between Students 1 st schools Subsequent schools τ r τ uj1 τ u j (0.14) (0.349) (0.414) τ r τ u (0.4) (0.184) τ r τ u (0.54) (0.49) 7
28 Conditional Fixed Effects Estimating Model CCMM-PGM 1 First School-PGM 1 Delete-PGM Parameter Coeff. Est. (SE) Coeff. Est. (SE) Coeff. Est. (SE) Model for initial status Intercept FEMALE Sch1_PRIVATE Model for 1 st slope Intercept FEMALE Sch1_PRIVATE SubSch_PRIVATE Model for nd slope Intercept FEMALE Sch1_PRIVATE SubSch_PRIVATE γ 0000 γ 0100 γ 0010 γ 1000 γ 1100 γ 1010 γ 1001 γ 000 γ 100 γ 010 γ (0.18) (0.176) (0.385) (0.163) (0.186) (1.59) (1.569) (0.07) (0.078) (0.407) (0.45) γ 000 γ 010 γ 001 γ 100 γ 110 γ 101 γ 00 γ 10 γ (0.168) (0.187) (0.374) (0.156) (0.185) (0.340) (0.073) (0.077) (0.168) γ 000 γ 010 γ 001 γ 100 γ 110 γ 101 γ 00 γ 10 γ DIC 454, , ,754.3 (0.04) (0.16) (0.43) (0.197) (0.7) (0.403) (0.086) (0.096) (0.180) 8
29 Conditional Random Effects Estimating Model CCMM-PGM First School-PGM Delete-PGM Parameter Coeff. Est. (SE) Coeff. Est. (SE) Coeff. Est. (SE) Level-1 variance between Measures σ (0.657) σ (0.657) σ (0.711) Initial status variance between Students 1 st schools τ r0 τ u0j (0.97) (1.041) τ r0 τ u (1.043) (1.08) τ r0 τ u (1.95) (1.15) 1 st slope variance between Students 1 st schools Subsequent schools τ r1 τ u1j1 τ u1 j (1.107) (1.089) (0.93) τ r1 τ u (1.08) (0.908) τ r1 τ u (1.337) (1.111) nd slope variance between Students 1 st schools Subsequent schools τ r τ uj1 τ u j (0.17) (0.353) (0.395) τ r τ u (0.0) (0.190) τ r τ u (0.55) (0.40) 9
30 Implications Ignoring mobility could lead to inaccurate conclusions about: The intercept and slope estimates in a three-level PGM if one were to delete mobile cases The impact of cluster-level predictors (regardless if you delete or ignore mobile individuals) The impact of both level- and level-3 predictors if one were to delete mobile cases. 30
31 Implications Researchers using a PGM ignoring multiple membership data should be careful when making inferences about the nature of variability in growth rates. For the delete-pgm, the SE estimates of the other variances were larger, which could then lead to erroneous conclusions about random effects if mobile individuals were removed from analysis. The CCMM-PGM fit to the data better than the first school-pgm. Because of these findings, a simulation study will be conducted this summer, so stay tuned 31
32 Thank you! 3
Ron Heck, Fall Week 3: Notes Building a Two-Level Model
Ron Heck, Fall 2011 1 EDEP 768E: Seminar on Multilevel Modeling rev. 9/6/2011@11:27pm Week 3: Notes Building a Two-Level Model We will build a model to explain student math achievement using student-level
More informationIntroducing Generalized Linear Models: Logistic Regression
Ron Heck, Summer 2012 Seminars 1 Multilevel Regression Models and Their Applications Seminar Introducing Generalized Linear Models: Logistic Regression The generalized linear model (GLM) represents and
More 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 informationRon Heck, Fall Week 8: Introducing Generalized Linear Models: Logistic Regression 1 (Replaces prior revision dated October 20, 2011)
Ron Heck, Fall 2011 1 EDEP 768E: Seminar in Multilevel Modeling rev. January 3, 2012 (see footnote) Week 8: Introducing Generalized Linear Models: Logistic Regression 1 (Replaces prior revision dated October
More informationAn Introduction to Mplus and Path Analysis
An Introduction to Mplus and Path Analysis PSYC 943: Fundamentals of Multivariate Modeling Lecture 10: October 30, 2013 PSYC 943: Lecture 10 Today s Lecture Path analysis starting with multivariate regression
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 informationAn Introduction to Path Analysis
An Introduction to Path Analysis PRE 905: Multivariate Analysis Lecture 10: April 15, 2014 PRE 905: Lecture 10 Path Analysis Today s Lecture Path analysis starting with multivariate regression then arriving
More informationPart 8: GLMs and Hierarchical LMs and GLMs
Part 8: GLMs and Hierarchical LMs and GLMs 1 Example: Song sparrow reproductive success Arcese et al., (1992) provide data on a sample from a population of 52 female song sparrows studied over the course
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 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 informationPath Analysis. PRE 906: Structural Equation Modeling Lecture #5 February 18, PRE 906, SEM: Lecture 5 - Path Analysis
Path Analysis PRE 906: Structural Equation Modeling Lecture #5 February 18, 2015 PRE 906, SEM: Lecture 5 - Path Analysis Key Questions for Today s Lecture What distinguishes path models from multivariate
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 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 informationSociology 593 Exam 2 Answer Key March 28, 2002
Sociology 59 Exam Answer Key March 8, 00 I. True-False. (0 points) Indicate whether the following statements are true or false. If false, briefly explain why.. A variable is called CATHOLIC. This probably
More informationA multivariate multilevel model for the analysis of TIMMS & PIRLS data
A multivariate multilevel model for the analysis of TIMMS & PIRLS data European Congress of Methodology July 23-25, 2014 - Utrecht Leonardo Grilli 1, Fulvia Pennoni 2, Carla Rampichini 1, Isabella Romeo
More informationMeasurement Invariance (MI) in CFA and Differential Item Functioning (DIF) in IRT/IFA
Topics: Measurement Invariance (MI) in CFA and Differential Item Functioning (DIF) in IRT/IFA What are MI and DIF? Testing measurement invariance in CFA Testing differential item functioning in IRT/IFA
More informationBayesian non-parametric model to longitudinally predict churn
Bayesian non-parametric model to longitudinally predict churn Bruno Scarpa Università di Padova Conference of European Statistics Stakeholders Methodologists, Producers and Users of European Statistics
More informationThe Application and Promise of Hierarchical Linear Modeling (HLM) in Studying First-Year Student Programs
The Application and Promise of Hierarchical Linear Modeling (HLM) in Studying First-Year Student Programs Chad S. Briggs, Kathie Lorentz & Eric Davis Education & Outreach University Housing Southern Illinois
More informationMultilevel Statistical Models: 3 rd edition, 2003 Contents
Multilevel Statistical Models: 3 rd edition, 2003 Contents Preface Acknowledgements Notation Two and three level models. A general classification notation and diagram Glossary Chapter 1 An introduction
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 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 informationSummary of Talk Background to Multilevel modelling project. What is complex level 1 variation? Tutorial dataset. Method 1 : Inverse Wishart proposals.
Modelling the Variance : MCMC methods for tting multilevel models with complex level 1 variation and extensions to constrained variance matrices By Dr William Browne Centre for Multilevel Modelling Institute
More informationExploiting TIMSS and PIRLS combined data: multivariate multilevel modelling of student achievement
Exploiting TIMSS and PIRLS combined data: multivariate multilevel modelling of student achievement Second meeting of the FIRB 2012 project Mixture and latent variable models for causal-inference and analysis
More informationTwo-Level Models for Clustered* Data
Two-Level Models for Clustered* Data Today s Class: Fixed vs. Random Effects for Modeling Clustered Data ICC and Design Effects in Clustered Data Group-Mean-Centering vs. Grand-Mean Centering Model Extensions
More informationHierarchical Linear Modeling. Lesson Two
Hierarchical Linear Modeling Lesson Two Lesson Two Plan Multivariate Multilevel Model I. The Two-Level Multivariate Model II. Examining Residuals III. More Practice in Running HLM I. The Two-Level Multivariate
More informationMachine Learning Linear Regression. Prof. Matteo Matteucci
Machine Learning Linear Regression Prof. Matteo Matteucci Outline 2 o Simple Linear Regression Model Least Squares Fit Measures of Fit Inference in Regression o Multi Variate Regession Model Least Squares
More informationOutline. Clustering. Capturing Unobserved Heterogeneity in the Austrian Labor Market Using Finite Mixtures of Markov Chain Models
Capturing Unobserved Heterogeneity in the Austrian Labor Market Using Finite Mixtures of Markov Chain Models Collaboration with Rudolf Winter-Ebmer, Department of Economics, Johannes Kepler University
More informationA Study of Statistical Power and Type I Errors in Testing a Factor Analytic. Model for Group Differences in Regression Intercepts
A Study of Statistical Power and Type I Errors in Testing a Factor Analytic Model for Group Differences in Regression Intercepts by Margarita Olivera Aguilar A Thesis Presented in Partial Fulfillment of
More informationCourse Introduction and Overview Descriptive Statistics Conceptualizations of Variance Review of the General Linear Model
Course Introduction and Overview Descriptive Statistics Conceptualizations of Variance Review of the General Linear Model EPSY 905: Multivariate Analysis Lecture 1 20 January 2016 EPSY 905: Lecture 1 -
More informationTime-Invariant Predictors in Longitudinal Models
Time-Invariant Predictors in Longitudinal Models Today s Topics: What happens to missing predictors Effects of time-invariant predictors Fixed vs. systematically varying vs. random effects Model building
More informationTABLE OF CONTENTS INTRODUCTION TO MIXED-EFFECTS MODELS...3
Table of contents TABLE OF CONTENTS...1 1 INTRODUCTION TO MIXED-EFFECTS MODELS...3 Fixed-effects regression ignoring data clustering...5 Fixed-effects regression including data clustering...1 Fixed-effects
More informationTime-Invariant Predictors in Longitudinal Models
Time-Invariant Predictors in Longitudinal Models Topics: Summary of building unconditional models for time Missing predictors in MLM Effects of time-invariant predictors Fixed, systematically varying,
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 informationReview of Multiple Regression
Ronald H. Heck 1 Let s begin with a little review of multiple regression this week. Linear models [e.g., correlation, t-tests, analysis of variance (ANOVA), multiple regression, path analysis, multivariate
More 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 informationAdditional Notes: Investigating a Random Slope. When we have fixed level-1 predictors at level 2 we show them like this:
Ron Heck, Summer 01 Seminars 1 Multilevel Regression Models and Their Applications Seminar Additional Notes: Investigating a Random Slope We can begin with Model 3 and add a Random slope parameter. If
More informationTime-Invariant Predictors in Longitudinal Models
Time-Invariant Predictors in Longitudinal Models Today s Class (or 3): Summary of steps in building unconditional models for time What happens to missing predictors Effects of time-invariant predictors
More informationMultilevel/Mixed Models and Longitudinal Analysis Using Stata
Multilevel/Mixed Models and Longitudinal Analysis Using Stata Isaac J. Washburn PhD Research Associate Oregon Social Learning Center Summer Workshop Series July 2010 Longitudinal Analysis 1 Longitudinal
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 informationPart 7: Hierarchical Modeling
Part 7: Hierarchical Modeling!1 Nested data It is common for data to be nested: i.e., observations on subjects are organized by a hierarchy Such data are often called hierarchical or multilevel For example,
More informationAP Statistics Bivariate Data Analysis Test Review. Multiple-Choice
Name Period AP Statistics Bivariate Data Analysis Test Review Multiple-Choice 1. The correlation coefficient measures: (a) Whether there is a relationship between two variables (b) The strength of the
More informationMixed- Model Analysis of Variance. Sohad Murrar & Markus Brauer. University of Wisconsin- Madison. Target Word Count: Actual Word Count: 2755
Mixed- Model Analysis of Variance Sohad Murrar & Markus Brauer University of Wisconsin- Madison The SAGE Encyclopedia of Educational Research, Measurement and Evaluation Target Word Count: 3000 - Actual
More informationClass Notes: Week 8. Probit versus Logit Link Functions and Count Data
Ronald Heck Class Notes: Week 8 1 Class Notes: Week 8 Probit versus Logit Link Functions and Count Data This week we ll take up a couple of issues. The first is working with a probit link function. While
More informationGeneralized Linear Models for Non-Normal Data
Generalized Linear Models for Non-Normal Data Today s Class: 3 parts of a generalized model Models for binary outcomes Complications for generalized multivariate or multilevel models SPLH 861: Lecture
More informationTime-Invariant Predictors in Longitudinal Models
Time-Invariant Predictors in Longitudinal Models Topics: What happens to missing predictors Effects of time-invariant predictors Fixed vs. systematically varying vs. random effects Model building strategies
More informationReview of the General Linear Model
Review of the General Linear Model EPSY 905: Multivariate Analysis Online Lecture #2 Learning Objectives Types of distributions: Ø Conditional distributions The General Linear Model Ø Regression Ø Analysis
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 informationCourse Introduction and Overview Descriptive Statistics Conceptualizations of Variance Review of the General Linear Model
Course Introduction and Overview Descriptive Statistics Conceptualizations of Variance Review of the General Linear Model PSYC 943 (930): Fundamentals of Multivariate Modeling Lecture 1: August 22, 2012
More informationSparse Matrix Methods and Mixed-effects Models
Sparse Matrix Methods and Mixed-effects Models Douglas Bates University of Wisconsin Madison 2010-01-27 Outline Theory and Practice of Statistics Challenges from Data Characteristics Linear Mixed Model
More informationREVIEW 8/2/2017 陈芳华东师大英语系
REVIEW Hypothesis testing starts with a null hypothesis and a null distribution. We compare what we have to the null distribution, if the result is too extreme to belong to the null distribution (p
More informationSociology 593 Exam 2 March 28, 2002
Sociology 59 Exam March 8, 00 I. True-False. (0 points) Indicate whether the following statements are true or false. If false, briefly explain why.. A variable is called CATHOLIC. This probably means that
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 informationUsing Decision Trees to Evaluate the Impact of Title 1 Programs in the Windham School District
Using Decision Trees to Evaluate the Impact of Programs in the Windham School District Box 41071 Lubbock, Texas 79409-1071 Phone: (806) 742-1958 www.depts.ttu.edu/immap/ Stats Camp offers week long sessions
More informationPackage effectfusion
Package November 29, 2016 Title Bayesian Effect Fusion for Categorical Predictors Version 1.0 Date 2016-11-21 Author Daniela Pauger [aut, cre], Helga Wagner [aut], Gertraud Malsiner-Walli [aut] Maintainer
More informationMULTILEVEL IMPUTATION 1
MULTILEVEL IMPUTATION 1 Supplement B: MCMC Sampling Steps and Distributions for Two-Level Imputation This document gives technical details of the full conditional distributions used to draw regression
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 informationAddressing Analysis Issues REGRESSION-DISCONTINUITY (RD) DESIGN
Addressing Analysis Issues REGRESSION-DISCONTINUITY (RD) DESIGN Overview Assumptions of RD Causal estimand of interest Discuss common analysis issues In the afternoon, you will have the opportunity to
More informationUnivariate Normal Distribution; GLM with the Univariate Normal; Least Squares Estimation
Univariate Normal Distribution; GLM with the Univariate Normal; Least Squares Estimation PRE 905: Multivariate Analysis Spring 2014 Lecture 4 Today s Class The building blocks: The basics of mathematical
More informationHierarchical Generalized Linear Models. ERSH 8990 REMS Seminar on HLM Last Lecture!
Hierarchical Generalized Linear Models ERSH 8990 REMS Seminar on HLM Last Lecture! Hierarchical Generalized Linear Models Introduction to generalized models Models for binary outcomes Interpreting parameter
More informationAssessing the relation between language comprehension and performance in general chemistry. Appendices
Assessing the relation between language comprehension and performance in general chemistry Daniel T. Pyburn a, Samuel Pazicni* a, Victor A. Benassi b, and Elizabeth E. Tappin c a Department of Chemistry,
More informationTime Invariant Predictors in Longitudinal Models
Time Invariant Predictors in Longitudinal Models Longitudinal Data Analysis Workshop Section 9 University of Georgia: Institute for Interdisciplinary Research in Education and Human Development Section
More 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 informationDuration of Unemployment - Analysis of Deviance Table for Nested Models
Duration of Unemployment - Analysis of Deviance Table for Nested Models February 8, 2012 The data unemployment is included as a contingency table. The response is the duration of unemployment, gender and
More informationFinal Exam - Solutions
Ecn 102 - Analysis of Economic Data University of California - Davis March 19, 2010 Instructor: John Parman Final Exam - Solutions You have until 5:30pm to complete this exam. Please remember to put your
More informationappstats8.notebook October 11, 2016
Chapter 8 Linear Regression Objective: Students will construct and analyze a linear model for a given set of data. Fat Versus Protein: An Example pg 168 The following is a scatterplot of total fat versus
More informationAdvanced Quantitative Data Analysis
Chapter 24 Advanced Quantitative Data Analysis Daniel Muijs Doing Regression Analysis in SPSS When we want to do regression analysis in SPSS, we have to go through the following steps: 1 As usual, we choose
More informationIntroduction to Bayesian Statistics and Markov Chain Monte Carlo Estimation. EPSY 905: Multivariate Analysis Spring 2016 Lecture #10: April 6, 2016
Introduction to Bayesian Statistics and Markov Chain Monte Carlo Estimation EPSY 905: Multivariate Analysis Spring 2016 Lecture #10: April 6, 2016 EPSY 905: Intro to Bayesian and MCMC Today s Class An
More informationLogistic Regression: Regression with a Binary Dependent Variable
Logistic Regression: Regression with a Binary Dependent Variable LEARNING OBJECTIVES Upon completing this chapter, you should be able to do the following: State the circumstances under which logistic regression
More informationFixed Effects, Invariance, and Spatial Variation in Intergenerational Mobility
American Economic Review: Papers & Proceedings 2016, 106(5): 400 404 http://dx.doi.org/10.1257/aer.p20161082 Fixed Effects, Invariance, and Spatial Variation in Intergenerational Mobility By Gary Chamberlain*
More informationSTA 303 H1S / 1002 HS Winter 2011 Test March 7, ab 1cde 2abcde 2fghij 3
STA 303 H1S / 1002 HS Winter 2011 Test March 7, 2011 LAST NAME: FIRST NAME: STUDENT NUMBER: ENROLLED IN: (circle one) STA 303 STA 1002 INSTRUCTIONS: Time: 90 minutes Aids allowed: calculator. Some formulae
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 informationDescribing Nonlinear Change Over Time
Describing Nonlinear Change Over Time Longitudinal Data Analysis Workshop Section 8 University of Georgia: Institute for Interdisciplinary Research in Education and Human Development Section 8: Describing
More informationStat 135 Fall 2013 FINAL EXAM December 18, 2013
Stat 135 Fall 2013 FINAL EXAM December 18, 2013 Name: Person on right SID: Person on left There will be one, double sided, handwritten, 8.5in x 11in page of notes allowed during the exam. The exam is closed
More informationThe Wishart distribution Scaled Wishart. Wishart Priors. Patrick Breheny. March 28. Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/11
Wishart Priors Patrick Breheny March 28 Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/11 Introduction When more than two coefficients vary, it becomes difficult to directly model each element
More informationLinear Models 1. Isfahan University of Technology Fall Semester, 2014
Linear Models 1 Isfahan University of Technology Fall Semester, 2014 References: [1] G. A. F., Seber and A. J. Lee (2003). Linear Regression Analysis (2nd ed.). Hoboken, NJ: Wiley. [2] A. C. Rencher and
More informationBayesian Causal Mediation Analysis for Group Randomized Designs with Homogeneous and Heterogeneous Effects: Simulation and Case Study
Multivariate Behavioral Research, 50:316 333, 2015 Copyright C Taylor & Francis Group, LLC ISSN: 0027-3171 print / 1532-7906 online DOI: 10.1080/00273171.2014.1003770 Bayesian Causal Mediation Analysis
More informationCIVL 7012/8012. Collection and Analysis of Information
CIVL 7012/8012 Collection and Analysis of Information Uncertainty in Engineering Statistics deals with the collection and analysis of data to solve real-world problems. Uncertainty is inherent in all real
More informationSUPPLEMENTARY APPENDIX MATERIALS [not intended for publication]
SUPPLEMENTARY APPENDIX MATERIALS [not intended for publication] Appendix Figure A1: Schooling and democracy trends in Kenya Panel A: Proportion of population completing primary schooling in Kenya, by birth
More informationChapter 8: Sampling Distributions. A survey conducted by the U.S. Census Bureau on a continual basis. Sample
Chapter 8: Sampling Distributions Section 8.1 Distribution of the Sample Mean Frequently, samples are taken from a large population. Example: American Community Survey (ACS) A survey conducted by the U.S.
More informationBayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features. Yangxin Huang
Bayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features Yangxin Huang Department of Epidemiology and Biostatistics, COPH, USF, Tampa, FL yhuang@health.usf.edu January
More informationChapter 4 Regression with Categorical Predictor Variables Page 1. Overview of regression with categorical predictors
Chapter 4 Regression with Categorical Predictor Variables Page. Overview of regression with categorical predictors 4-. Dummy coding 4-3 4-5 A. Karpinski Regression with Categorical Predictor Variables.
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 informationPublished online: 19 Jun 2015.
This article was downloaded by: [University of Wisconsin - Madison] On: 23 June 2015, At: 12:17 Publisher: Routledge Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office:
More informationUtilizing Hierarchical Linear Modeling in Evaluation: Concepts and Applications
Utilizing Hierarchical Linear Modeling in Evaluation: Concepts and Applications C.J. McKinney, M.A. Antonio Olmos, Ph.D. Kate DeRoche, M.A. Mental Health Center of Denver Evaluation 2007: Evaluation and
More informationEPSY 905: Fundamentals of Multivariate Modeling Online Lecture #7
Introduction to Generalized Univariate Models: Models for Binary Outcomes EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #7 EPSY 905: Intro to Generalized In This Lecture A short review
More informationA Bayesian multi-dimensional couple-based latent risk model for infertility
A Bayesian multi-dimensional couple-based latent risk model for infertility Zhen Chen, Ph.D. Eunice Kennedy Shriver National Institute of Child Health and Human Development National Institutes of Health
More informationTechnical and Practical Considerations in applying Value Added Models to estimate teacher effects
Technical and Practical Considerations in applying Value Added Models to estimate teacher effects Pete Goldschmidt, Ph.D. Educational Psychology and Counseling, Development Learning and Instruction Purpose
More informationBayesian Networks in Educational Assessment
Bayesian Networks in Educational Assessment Estimating Parameters with MCMC Bayesian Inference: Expanding Our Context Roy Levy Arizona State University Roy.Levy@asu.edu 2017 Roy Levy MCMC 1 MCMC 2 Posterior
More informationOne-Way ANOVA. Some examples of when ANOVA would be appropriate include:
One-Way ANOVA 1. Purpose Analysis of variance (ANOVA) is used when one wishes to determine whether two or more groups (e.g., classes A, B, and C) differ on some outcome of interest (e.g., an achievement
More informationSimple Linear Regression: One Qualitative IV
Simple Linear Regression: One Qualitative IV 1. Purpose As noted before regression is used both to explain and predict variation in DVs, and adding to the equation categorical variables extends regression
More informationIntroduction to Matrix Algebra and the Multivariate Normal Distribution
Introduction to Matrix Algebra and the Multivariate Normal Distribution Introduction to Structural Equation Modeling Lecture #2 January 18, 2012 ERSH 8750: Lecture 2 Motivation for Learning the Multivariate
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 informationMultilevel Analysis, with Extensions
May 26, 2010 We start by reviewing the research on multilevel analysis that has been done in psychometrics and educational statistics, roughly since 1985. The canonical reference (at least I hope so) is
More informationMultilevel Modeling Day 2 Intermediate and Advanced Issues: Multilevel Models as Mixed Models. Jian Wang September 18, 2012
Multilevel Modeling Day 2 Intermediate and Advanced Issues: Multilevel Models as Mixed Models Jian Wang September 18, 2012 What are mixed models The simplest multilevel models are in fact mixed models:
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 informationSpatio-Temporal Threshold Models for Relating UV Exposures and Skin Cancer in the Central United States
Spatio-Temporal Threshold Models for Relating UV Exposures and Skin Cancer in the Central United States Laura A. Hatfield and Bradley P. Carlin Division of Biostatistics School of Public Health University
More informationTutorial 6: Tutorial on Translating between GLIMMPSE Power Analysis and Data Analysis. Acknowledgements:
Tutorial 6: Tutorial on Translating between GLIMMPSE Power Analysis and Data Analysis Anna E. Barón, Keith E. Muller, Sarah M. Kreidler, and Deborah H. Glueck Acknowledgements: The project was supported
More informationThree-Level Modeling for Factorial Experiments With Experimentally Induced Clustering
Three-Level Modeling for Factorial Experiments With Experimentally Induced Clustering John J. Dziak The Pennsylvania State University Inbal Nahum-Shani The University of Michigan Copyright 016, Penn State.
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 information1 Introduction. 2 Example
Statistics: Multilevel modelling Richard Buxton. 2008. Introduction Multilevel modelling is an approach that can be used to handle clustered or grouped data. Suppose we are trying to discover some of the
More information