In milanwiedemann/minvariance: Testing Measurement Invariance

• Lavaan tutorial Link (by Julie Wood)

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

library(minvariance)
library(tidyverse)
library(lavaan)
library(semPlot)

Load data

data_url <- "https://quantdev.ssri.psu.edu/sites/qdev/files/ECLS_Science.csv"

df <- read_csv(file = url(data_url)) %>% 
  select("id", 
         "s_g3", "r_g3", "m_g3", 
         "s_g5", "r_g5", "m_g5", 
         "s_g8", "r_g8", "m_g8")

# Rename so functions work
timepoints <- list(t3 = c("s_g3", "r_g3", "m_g3"), 
                   t5 = c("s_g5", "r_g5", "m_g5"), 
                   t8 = c("s_g8", "r_g8", "m_g8"))

lavaan objects

minvariance(data = df, 
            var_list = timepoints, 
            return = "lavaan_objects")

Fit statistics

minvariance(data = df, 
            var_list = timepoints, 
            return = "fit_statistics") %>% 
  knitr::kable(digits = 3)

Model tests

minvariance(data = df, 
                 var_list = timepoints, 
                 return = "model_tests") %>% 
  knitr::kable(digits = 3)

lavaan syntax

lminvar_syntax <- minvariance(data = df, 
                              var_list = timepoints, 
                              return = "lavaan_syntax") 

Configural

lminvar_syntax$configural %>% 
  cat()

Weak

lminvar_syntax$weak %>% 
  cat()

Strong

lminvar_syntax$strong %>% 
  cat()

Strict

lminvar_syntax$strict %>% 
  cat()

milanwiedemann/minvariance documentation built on May 9, 2021, 9:05 p.m.