model <- 'Y~X'
stability <- data.frame(X = c(.3, .4, .5), Y = c(.3, .5, .6))
dat <- data.frame(Y = rnorm(500, 0, 1), X = rnorm(500, 0, 1), Z = rnorm(500, 0, 1))
output <- stim(data = dat, model = model, stability = stability)
lavaanSummary(output, subset = c(1,2))
Run the code above in your browser using DataLab