Applied
This presentation is the property of its rightful owner.
Sponsored Links
1 / 116

Joseph J. Sudano, Jr., PhD Center for Health Care Research and Policy PowerPoint PPT Presentation


  • 79 Views
  • Uploaded on
  • Presentation posted in: General

Applied Structural Equation Modeling for Dummies, by Dummies February 22, 2013 Indiana University, Bloomington. Joseph J. Sudano, Jr., PhD Center for Health Care Research and Policy Case Western Reserve University at The MetroHealth System Adam T. Perzynski , PhD

Download Presentation

Joseph J. Sudano, Jr., PhD Center for Health Care Research and Policy

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Joseph j sudano jr phd center for health care research and policy

Applied Structural Equation Modeling for Dummies, by DummiesFebruary 22, 2013Indiana University, Bloomington

Joseph J. Sudano, Jr., PhD

Center for Health Care Research and Policy

Case Western Reserve University at The MetroHealth System

Adam T. Perzynski, PhD

Center for Health Care Research and Policy

Case Western Reserve University at The MetroHealth System


Acknowledgements

Acknowledgements

  • Thanks Joe.

  • Thanks to Bill Pridemore and all of you here at IU.

  • Thanks to Doug Gunzler.

  • Thanks to Kyle Kercher.


Rejected titles for this talk february 22 2013 indiana university bloomington

Rejected Titles for this TalkFebruary 22, 2013Indiana University, Bloomington

Joseph J. Sudano, Jr., PhD

Center for Health Care Research and Policy

Case Western Reserve University at The MetroHealth System

Adam T. Perzynski, PhD

Center for Health Care Research and Policy

Case Western Reserve University at The MetroHealth System


Structural equation modeling for fashion week

Structural Equation Modeling for Fashion Week


We have lots of models

We have lots of Models!


Structural equation modelin fer pirates

Structural Equation Modelin’ fer Pirates


Structural equation modelin fer pirates1

Structural Equation Modelin’ fer Pirates

SEM be a statistical technique for testin' and estimatin' causal relations usin' a combination o' statistical data and qualitative causal assumptions

*From Wikipedia


Assumptions

Assumptions

  • I do not actually assume you are dummies

  • Feel free to assume what you want about me

  • I do not assume you will be experts in SEM after this presentation

  • I assume you know something about means and regression (hopefully)


Outline

Outline

  • Important SEM Resources

  • Measurement (and measurement error)

  • Examples

    • Measurement Invariance

    • Latent Class Analysis

    • Latent Growth Mixture Modeling

  • Model Specification


Outline1

Outline

  • Important SEM Resources

  • Measurement (and measurement error)

  • Examples

    • Measurement Invariance

    • Latent Class Analysis

    • Latent Growth Mixture Modeling

  • Model Specification


Sem resources

SEM Resources


Sem resources1

SEM Resources


Sem resources2

SEM Resources


Sem resources statmodel com

SEM Resources: Statmodel.com


Sem resources3

SEM Resources


Sem resources4

SEM Resources


Sem resources5

SEM Resources


Sem resources6

SEM Resources


Sem resources7

SEM Resources


Sem resources8

SEM Resources


Sem resources9

SEM Resources


Outline2

Outline

  • Important SEM Resources

  • Measurement (and measurement error)

  • Examples

    • Measurement Invariance

    • Latent Class Analysis

    • Latent Growth Mixture Modeling

  • Model Specification


Measurement models

Measurement Models

  • A special type of causal models

  • Survey items are assumed to have measurement error

    • Each question has its own amount of error

  • Your answer to a survey question is causally related to a latent, unobserved variable.


Perfect measurement

health

Self-rated health

Perfect Measurement

1.0?


Causality and the latent concept of health

Causality and the Latent Concept of Health

  • In general, how would you describe your health?

  • We assume that every individual varies along an infinite continuum from best possible health to worst possible health.

  • When any given individual answers this question, they are approximating their position on this latent continuum.


Imperfect measurement

health

Self-rated health

e4

Imperfect Measurement

1.0

Variance > 0

< 1.0


Measurement models using multiple indicators

Measurement Models using Multiple Indicators

  • Single items are unreliable

  • Single cases prevent generalizability

  • Use multiple indicators and large samples to estimate the values of the latent, unobservered variables or factors

  • The SF36 uses multiple indicators describing multiple factors in order to measure health more reliably.


Outline3

Outline

  • Important SEM Resources

  • Measurement (and measurement error)

  • Examples

    • Measurement Invariance

    • Latent Class Analysis

    • Latent Growth Mixture Modeling

  • Model Specification


Joseph j sudano jr phd center for health care research and policy

Acknowledgement: This study was funded by Grant number R01-AG022459 from the NIH National Institute on Aging.

Measuring Disparities: Bias in Self-reported Health Among Spanish-speaking PatientsJ.J. Sudano1,2, A.T. Perzynski1,2, T.E. Love2, S.A. Lewis1,B. Ruo3, D.W. Baker31 The MetroHealth System, Cleveland, OH; 2 Case Western Reserve University School of Medicine, Cleveland, OH; 3 Northwestern University Feinberg School of Medicine


Measurement model of the sf36

Measurement Model of the SF36


Objective significance

Objective & Significance

  • Do observed differences in SRH reflect true differences in health?

    • Cultural and language differences may create measurement bias

    • If outcomes aren’t measuring the same thing in different groups, we have a problem


Measurement equivalence factorial invariance

Measurement Equivalence &Factorial Invariance

  • It is only possible to properly interpret group differences after measurement equivalence has been established (Horn & McArdle, 1992; Steenkamp & Baumgartner, 1998).

  • “It may be the case that the groups differ … but it also may be the case that extraneous influences are giving rise to the observed difference.” Meredith & Teresi (2006 p. S69)

  • The external validity of any conclusion regarding group differences rests securely on whether the measurement equivalence of the scale has been established (Borsboom, 2006).


Cross sectional study

Cross-sectional Study

  • N= 1281

  • Medical patients categorized into four groups:White, Black, English-speaking Hispanic and Spanish-speaking Hispanic.

  • Multigroup Confirmatory Factor Analysis (MGCFA)


Two types of invariance

Two Types of Invariance

  • Metric (Weak) Invariance

    • Are the item factor loadings equivalent across groups?

    • Is a one unit change in the item equal to a one unit change in the factor score for all groups?

  • Scalar (Strong) Invariance

    • Are the item intercepts equivalent across groups?

    • Unequal intercepts results in unequal scaling of factor scores


Weak invariance

health

Self-rated health

e4

Weak invariance

What happens to the model fit when we constrain all of these paths (loadings) to be equal across groups?


Joseph j sudano jr phd center for health care research and policy

The Unconstrained Model Fits the Data Well


Joseph j sudano jr phd center for health care research and policy

The Unconstrained Model fits the data well

The model with factor loadings constrained still fits the data well.


Metric weak invariance was confirmed

Metric (Weak) Invariance was Confirmed


I forget what an intercept is

I forget what an intercept is

  • Scalar (Strong) Invariance

    • Are the item intercepts equivalent across groups?

  • Intercept: the intercept in a multiple regression model is the mean for the response when all of the explanatory variables take on the value 0.

  • Could be called the “starting point”


Joseph j sudano jr phd center for health care research and policy

The Unconstrained Model fits the data well

The model with factor loadings constrained still fits the data well.

Constraining the intercepts results in a worsening of model fit


Joseph j sudano jr phd center for health care research and policy

The model with factor loadings constrained still fits the data well.

Constraining the intercepts results in a worsening of model fit

The fit is still poor if you allow intercepts for English-speaking Hispanics to vary


Joseph j sudano jr phd center for health care research and policy

The model with factor loadings constrained still fits the data well.

The fit is acceptable if you allow intercepts for Spanish speaking Hispanics to vary


Scalar strong invariance is not confirmed

Scalar (Strong) Invariance is NOT Confirmed

Measurement equivalence of the SF36 does not exist for Spanish speaking Hispanics


Intercepts are lower for spanish speaking hispanics on nearly all items

Intercepts are lower for Spanish-speaking Hispanics on nearly all items

Measurement equivalence of the SF36 does not exist for Spanish speaking Hispanics


Joseph j sudano jr phd center for health care research and policy

Use of English Rating Categories on Twiter

Using of Spanish Rating Categories on Twitter


Outline4

Outline

  • Important SEM Resources

  • Measurement (and measurement error)

  • Examples

    • Measurement Invariance

    • Latent Class Analysis

    • Latent Growth Mixture Modeling

  • Model Specification


Everywhere and nowhere latent class analysis of knowledge of the spread of hepatitis c

Everywhere and Nowhere: Latent Class Analysis of Knowledge of the Spread of Hepatitis C

Adam T. Perzynski, PhD

E-mail: [email protected]


Introduction

Introduction

  • Hepatitis C is a widespread and serious disease that affects the liver.

  • 170 million people worldwide are infected.

  • 3.9 million Americans infected with HCV. (AHRQ 2003)

  • More Americans die every year from Chronic HCV infection than from HIV


Hcv transmission

HCV Transmission

  • Blood

    • Injection Drug Use

    • Blood Transfusions

    • Needle Sticks

    • Shared Household Items (Razor or Toothbrush)

  • Sexual transmission of HCV is recognized but is infrequent.

  • HCV is not transmitted by Coughing, Kissing, Sneezing, Touching, Bathrooms, Fecal Matter, or Contaminated Food


Sample

Sample

  • Behavior Risk Factor Surveillance System (BRFSS), 2001, Arizona

  • Conducted by the Centers for Disease Control (CDC)

  • The world’s largest telephone survey

  • Nearly 200,000 people participated in 2001


Measure

Measure

  • Do you think hepatitis C can be spread thru?

    • Sneezing or Coughing

    • Kissing

    • Unprotected Sex

    • Food or Water

    • Sharing Needles to Inject Street Drugs

    • Using the Same Bathroom

    • Contact with the Blood of an Infected Person


Methods of analysis

Methods of Analysis

  • Analyzed with Mplus

  • Analysis proceeded in several stages

    • Exploratory Factor Analysis

    • Confirmatory Factor Analysis

    • Cluster Analysis (Not reported)

    • Latent Class Analysis

    • Mixture Modeling

  • Robust estimation for binary indicators

  • Missing Values Imputation using Full Information Maximum Likelihood Estimation (FIML)


Distribution of outcome variables

Distribution of Outcome Variables


Means and standard deviations

Means and Standard Deviations


Exploratory factor analysis

Exploratory Factor Analysis

  • Scree plot, Eigenvalues, and Root Mean Square Residuals more or less supported a two factor solution


What is different about lca

What is different about LCA?

  • Instead of assuming that the latent variable is continuous (infinitely poor to infinitely good)

  • We assume the latent variable is categorical.

  • Membership in “hidden” empirical forms determines answers rather than a single latent continuum.


Three latent classes

Three Latent Classes

  • The Two Category and Four Category models do not fit the data as well as as the Three Category model.

  • HCV is Nowhere

    • N = 1683 (The largest class!)

  • Full Awareness of how HCV is Spread

    • N = 930

  • HCV is Everywhere

    • N = 479


Joseph j sudano jr phd center for health care research and policy

Figure 3: Estimated Probabilities of Knowing How HCV is Spread by Class Membership


Additional analyses

Additional Analyses

  • What predicts membership in each latent class?

  • Do the relationships between variables vary inside of a particular class?

  • Mixture Modeling

    • Simultaneously test continuous and categorical predictors of class membership.


Outline5

Outline

  • Important SEM Resources

  • Measurement (and measurement error)

  • Examples

    • Measurement Invariance

    • Latent Class Analysis

    • Latent Growth Mixture Modeling

  • Model Specification


Longitudinal patterns of depressive symptoms in the health and retirement study

Longitudinal Patterns of Depressive Symptoms in the Health and Retirement Study

Adam T. Perzynski, PhD & Joseph S. Sudano, Jr., PhD

Center for Health Care Research and Policy

Case Western Reserve University and MetroHealth

Presentation at the Annual Meeting of the Gerontological Society of America on November 22, 2010


Introduction1

Introduction

  • This is another measurement study

  • Explore the use of Latent Class Growth Analysis to model changes in depressive symptoms over time in the Health and Retirement Study.

  • Most studies compare the change in means scores between two waves.

  • A small number of studies have modeled change as a single growth trajectory


Change in means between waves

Change in Means Between Waves

  • Often we simply calculate the mean depressive symptoms at Wave 1 (baseline).

  • Subract it from the mean at Wave 2 (followup).


What is a trajectory

What is a trajectory?

  • Regrettably, the term “trajectory” has taken on multiple meanings across disciplines and research studies.

  • A broad, inclusive definition of trajectory modeling is the analysis of patterns of change or stability.

  • Confusion is possible between aggregate trajectories which summarize an overall average pattern of change for a population and disaggregated trajectories which examine multiple potential trajectories of different shapes (George 2006).


Example of a single growth trajectory

Example of a single growth trajectory


Continuous latent growth curve analysis

Continuous Latent Growth Curve Analysis

  • LGA / LGCA

  • Studies in older adults (ie George and Lynch 2003) typically find that the slope of the latent growth curve for depressive symptoms is small and positive, and that the slope of the curve is steepest in the oldest cohorts.


Example from george and lynch 2003

Example from George and Lynch (2003)


Example of an lga finding

Example of an LGA finding


Lga estimates a single aggregate trajectory

LGA estimates a single Aggregate trajectory

  • Assumes that the average population starting point (intercept for the growth curve) and average amount of change (slope) are a sufficient depiction of variation over time in depressive symptoms.

  • If discrete subtypes of depressive symptom trajectories exist, but are ignored (as in single latent growth curve and autoregressive models) the magnitude of associations could be grossly misestimated.


What is latent class growth analysis

What is Latent Class Growth Analysis?

  • Latent Class Growth Analysis (LCGA), also referred to as growth mixture modeling, belongs to a family of statistical techniques referred to as general latent variable modeling or GLVM.


Why would we ever think we should use lcga

Why would we ever think we should use LCGA?

  • Studying the mean change or using a single trajectory for everyone assumes uniform heterogeneity in the population.

  • Researchers use familiar methods and typically assume that the underlying (latent or real) distribution of variables is continuous.

  • We have theoretical reasons to suspect that underlying distributions could be categorical.

  • Life course theorists (Dannefer) specifically caution that intracohort differentiation is unlikely to be homogeneous.


Why would we use lcga

Why would we use LCGA?

  • We think individuals and cohorts diverge over time

  • Cumulative change differentiates individuals and cohorts.


Prior lcga models of depression or depressive symptoms

Prior LCGA Models of Depression or Depressive Symptoms

  • LCGA models and closely related Longitudinal Latent Class Analysis (LLCA) have been used to estimate models of depressive symptoms in prior studies of

    • maternity (Campbell et al 2009; Mora et al 2009)

    • childhood and adolescence (Meadows et al 2006)

    • adolescence through young adulthood (Olino et al 2009)

    • response to antidepressants among adults (Muthen et al, 2007; Hunter et al 2009)

    • patients who have had a cardiovascular event (Kaptein et al 2006).


Methods

Methods

  • 5,195 age-eligible respondents from the 1992 Health and Retirement Study cohort, who completed interviews in all seven waves through 2004.

  • Depressive symptoms in HRS are measured using a dichotomous, 8-item version of the CES-D. Analysis begins with Wave 2 data due to a change in response categories from Wave 1.

  • Using MPlus, we compared the fit of LCGA models of two to eight classes while also accounting for the HRS complex sampling design.

  • We then tested the effect of a small number of covariates. This is very similar to a multinomial logistic regression.


Demographic characteristics

Demographic characteristics

  • Gender

    • 60.3% female

  • Race/ethnicity

    • 76.4% non-Hispanic White

    • 14.4% Black

    • 7.4% Hispanic

    • 1.8% other racial/ethnic groups

  • Age

    • Median=55

  • Education

    • Mean=12.4 years (SD=3.0).


Rule for determining the number of latent classes

Rule for Determining the number of Latent Classes

  • “How many trajectories are there?”

  • Measures of model fit including:

    • Lo-Mendell-Rubin Test (LMR test)

    • log-likelihood (LL)

    • Bayesian Information Criteria (BIC) (Vuong, 1989; Muthen, 2004; Muthen, & Muthen, 2005; Nylund et al, 2007).

  • Here we will use the LMR Test

  • Where k is the number of latent classes, this test gives a p-value for the k-1 versus the k-class model when running the k-class model (Vuong, 1989; Muthen, B. 2005).

  • The first time p > .05, k-1 is the preferred number of classes.


Results

Results

  • How many classes are there?

  • What do the classes look like?

  • How is this different from looking at means or single trajectory?

  • Are any demographic variables associated with being in a particular class?


How many classes are there

How many Classes are there?


How many classes are there1

How many Classes are there?


What do the classes look like

What do the classes look like?


How is this different from looking at means or a single trajectory

How is this different from looking at Means or a Single Trajectory?

Online at: http://spreadsheets.google.com/pub?key=0ApRkae54BRnudEYyUGdXZWlES3Z4VzZ6akNaOFFiekE&gid=5


Does anything influence the chances of being in a particular class

Does anything influence the chances of being in a particular class?


Does anything influence the chances of being in a particular class1

Does anything influence the chances of being in a particular class?

  • Females, African Americans and those with fewer years of education have a higher probability of being in the Many Symptoms trajectory.


Outline6

Outline

  • Important SEM Resources

  • Measurement (and measurement error)

  • Examples

    • Measurement Invariance

    • Latent Class Analysis

    • Latent Growth Mixture Modeling

  • Model Specification


Model specification

Model Specification

  • Choosing the model that best represents the data structure and addresses the research questions of interest can be a daunting task.

  • Brief overview of model specification tests and procedures.


Model specification1

Model Specification

  • “First, your return to shore was not part of our negotiations nor our agreement so I must do nothing. And secondly, you must be a pirate for the pirate's code to apply and you're not. And thirdly, the code is more what you'd call ‘guidelines’ than actual rules.”

  • Captain Barbossa from Pirates of the Caribbean: The Curse of the Black Pearl (2003)


Model specification2

Model Specification

  • In model specification a researcher can use:

    • logic, theory and prior empirical evidence to choose the initial model

    • model comparison testing to compare the initial model to competing models

    • a combination of theory, prior evidence, and the results of the model comparison testing to decide upon which model or models are appropriate for a given study


Nested or not nested

Nested or Not Nested?


Chi square test

Chi Square Test

  • The Chi-square statistic is computed and used to test whether the model does fit the data well.

  • It is the basis for most other fit tests.

  • Along with other fit tests we use it to evaluate whether to include or exclude model paths relating measures to each other for a given study.


Chi square test1

Chi Square Test

  • Also called the discrepancy function

  • If not significant, the model is regarded as acceptable.*


Chi square test2

Chi Square Test*

  • Some limitations are:

  • Complex models with many parameters

  • With large samples, models will most often be rejected, sometimes unfairly

  • Where multivariate non-normality is present, the chi-square fit index is inaccurate.

  • Modified tests (The Satorra-Bentler scaled chi-square) are available.


Modification indices

Modification Indices

  • Modification indices can be calculated individually for every path that is fixed to zero, by estimating a chi-square test statistic with one df.

  • The higher the value of the modification index for a causal path, the better the predicted improvement in overall model fit if that path were added to the model.

  • Jöreskogsuggested that a modification index should be at least five before the researcher considers adding the causal path and modifying the hypothesized model.


R squared

R-squared

  • In linear regression analysis, we interpret the r2 value as the amount of variation in the response that can be explained by the regressors in the model.

  • In SEM, it is pretty much the same*

  • *Not exactly, but that is beyond the “for dummies” version of this talk


Aic bic and bcc

AIC, BIC (and BCC)

  • Bayesian Information Criterion (BIC)

  • AkaikeInformation Criterion (AIC)

  • Based on the chi-squared test statistic

  • While the models under comparison can be nested or non-nested, in both these tests, as with all tests in this section, for a truly direct comparison, we prefer that the same observed measures are used in the models we are comparing.

  • Both BIC and AIC feature the goodness-of-fit term for our model , derived directly from the discrepancy function when applicable, along with a penalty term.


Aic bic and bcc1

AIC, BIC (and BCC)

  • Cannot identify if a model has good fit. Only if one model fits better than another.

  • The lower the value of BIC, AIC and BCC, the better the fit.

  • BCC penalizes for model complexity more than AIC and BIC.

  • BIC penalizes for model complexity more than AIC.


Specification search

Specification Search

  • Allows researchers to choose a model from among a number of candidates.

  • Exploratory

  • Should be guided by theory


Specification search1

Specification Search

  • Given the model in Figure 1, with 7 unknown paths, the number of models is equivalent to 27=128 possible specifications of the model.

  • 128 different possible models!


Specification search2

Specification Search

  • The unconstrained model (The one with all seven ambiguous paths in the model) demonstrates satisfactory overall model fit

    • CFI=.95

    • TLI=.92

    • RMSEA =.07

    • chisq= 476.33, DF=63


Table results of specification search in amos

Table. Results of Specification Search in AMOS

Notes: Reported R2 values are for the equation in each model with the endogenous PF latent variable with the interpretation of total explained variance in physical functioning given all other paths in the model. C is the chi-squared test statistic and df are the associated degrees of freedom.


Specification search3

Specification Search

  • When a number of models are plausible, specification tests can be used as evidence for verification of or improvement over an initial model.


Guidelines

‘Guidelines’

  • a researcher is ultimately left to decide if the results of the specification tests are unjustly in favor of a certain model due to complexity or sample size, rather than the meaning behind the causal paths.

  • Thus the specification tests act more like guidelines, rather than strict codes dictating the “best” fitting model.


Selected strengths limitations of sem

Selected Strengths & Limitations of SEM

  • Strengths

    • Very flexible

    • Estimate and correct for measurement error

  • Limitations

    • Large sample sizes

    • Challenging to learn

    • Need lots of hands-on experience to learn

    • Need a strong theoretical basis

      • It’s easy to mis-specify a model if you have no idea what you are doing.


Applied structural equation modeling for everyone february 22 2013 indiana university bloomington

Applied Structural Equation Modeling for Everyone!February 22, 2013Indiana University, Bloomington

Joseph J. Sudano, Jr., PhD

Center for Health Care Research and Policy

Case Western Reserve University at The MetroHealth System

Adam T. Perzynski, PhD

Center for Health Care Research and Policy

Case Western Reserve University at The MetroHealth System


  • Login