Generate lavaan syntax for latent change score models

In lcsm: Univariate and Bivariate Latent Change Score Modelling

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(lcsm)

Univariate LCSM syntax

The example below shows how to specify a generic univariate latent change score model using the function specify_uni_lcsm().

uni_lcsm_syntax <- specify_uni_lcsm(timepoints = 5,
                                    var = "x",  
                                    change_letter = "g",
                                    model = list(alpha_constant = TRUE, 
                                                 beta = TRUE, 
                                                 phi = TRUE))

# To get a readable output use cat() function
cat(uni_lcsm_syntax)

Add lavaan syntax

uni_lcsm_syntax_add <- specify_uni_lcsm(timepoints = 3,
                                    var = "x",  
                                    change_letter = "g",
                                    add = "# JUST ANOTHER COMMENT",
                                    model = list(alpha_constant = TRUE, 
                                                 beta = TRUE, 
                                                 phi = TRUE))

# To get a readable output use cat() function
cat(uni_lcsm_syntax_add)

Bivariate LCSM syntax

True Score to Change Score

# Specify bivariate LCSM
lavaan_bi_lcsm_delta_01 <- specify_bi_lcsm(timepoints = 5, 
                                           var_x = "x",
                                           model_x = list(alpha_constant = TRUE, 
                                                          beta = TRUE, 
                                                          phi = TRUE),
                                           var_y = "y",  
                                           model_y = list(alpha_constant = TRUE, 
                                                          beta = TRUE, 
                                                          phi = TRUE),  
                                           coupling = list(delta_lag_xy = TRUE, 
                                                           delta_lag_yx = TRUE),
                                           change_letter_x = "g",
                                           change_letter_y = "j")

# To get a readable output use cat() function
cat(lavaan_bi_lcsm_delta_01)

Change Score to Change Score

# Specify bivariate LCSM
lavaan_bi_lcsm_xi_01 <- specify_bi_lcsm(timepoints = 5, 
                                        var_x = "x",
                                        model_x = list(alpha_constant = TRUE, 
                                                       beta = TRUE, 
                                                       phi = TRUE),
                                        var_y = "y",  
                                        model_y = list(alpha_constant = TRUE, 
                                                       beta = TRUE, 
                                                       phi = TRUE),  
                                        coupling = list(xi_lag_xy = TRUE, 
                                                        xi_lag_yx = TRUE),
                                        change_letter_x = "g",
                                        change_letter_y = "j")

# To get a readable output use cat() function
cat(lavaan_bi_lcsm_xi_01)

Add model specifications to lavaan code

# Specify bivariate LCSM
lavaan_bi_lcsm_xi_add <- specify_bi_lcsm(timepoints = 3, 
                                        var_x = "x",
                                        model_x = list(alpha_constant = TRUE, 
                                                       beta = TRUE, 
                                                       phi = TRUE),
                                        var_y = "y",  
                                        model_y = list(alpha_constant = TRUE, 
                                                       beta = TRUE, 
                                                       phi = TRUE),  
                                        coupling = list(xi_lag_xy = TRUE, 
                                                        xi_lag_yx = TRUE),
                                        add = "age ~~ age_j2 * j2",
                                        change_letter_x = "g",
                                        change_letter_y = "j")

# To get a readable output use cat() function
cat(lavaan_bi_lcsm_xi_add)

TODOs

• show different examples
• highlight the comments that explain the different sections of the code
• show how this syntax can be adapted and used outside of this package in lavaan

lcsm documentation built on March 7, 2023, 6:45 p.m.