tests/testthat/test-standardizedSolution_boot_ci_eq.R

In semhelpinghands: Helper Functions for Structural Equation Modeling

# skip_on_cran()
# skip_if(!interactive(),
#         message = "standardizedSolution_boot_ci not tested if not interactive")

library(testthat)
library(semhelpinghands)

# Adapted from https://lavaan.ugent.be/tutorial/syntax2.html

library(lavaan)
data(HolzingerSwineford1939)
model <-
"
visual  =~ x1 + v2*x2 + v2*x3
textual =~ x4 + x5 + x6
speed   =~ x7 + x8 + x9
"

set.seed(1234)
system.time(fit <- cfa(model = model,
                       data = HolzingerSwineford1939,
                       se = "boot",
                       bootstrap = 50))
ci_boot <- standardizedSolution_boot_ci(fit,
                                        save_boot_est_std = TRUE,
                                        force_run = TRUE)
head(attr(ci_boot, "boot_est_std"))
get_std <- function(object) {
    lavaan::standardizedSolution(object)$est.std
  }

fit2 <- cfa(model = model,
            data = HolzingerSwineford1939,
            se = "none")
set.seed(1234)
boot_ci_test <- bootstrapLavaan(fit2, R = 50,
                                FUN = get_std)
set.seed(1234)
bl_est <- bootstrapLavaan(fit2, R = 50,
                                FUN = coef)
fit_bl_est <- fit
fit_bl_est@boot$coef <- bl_est
ci_boot_bl_est <- standardizedSolution_boot_ci(fit_bl_est,
                                        save_boot_est_std = TRUE,
                                        force_run = TRUE)
head(attr(ci_boot_bl_est, "boot_est_std"))

# Cannot compare with get_std results because, even with same seed,
# bootstrapLavaan and se="boot" does not result in the same set
# of bootstrap samples.

test_that("Compare boot estimates directly", {
    expect_equal(
        attr(ci_boot_bl_est, "boot_est_std"),
        boot_ci_test,
        ignore_attr = TRUE
      )
  })