tests/testthat/test-lav_betaselect_mod_boot_3.R

In betaselectr: Betas-Select in Structural Equation Models and Linear Models

skip_on_cran()
# Parallel processing
# Long test

# testthat::test_file("./tests/testthat/test_std_selected_lavaan_mod_boot_3.R")

library(manymome)

#Load a test data of 500 cases
data(data_test_medmod)
library(lavaan)
dat <- data_test_medmod
dat$group <- rep(c("gp1", "gp2"), nrow(dat))[seq_len(nrow(dat))]
dat$iv_mod <- dat$iv * dat$mod
head(dat)

mod <-
"
dv ~ med + cov2
"
fit <- sem(mod,
           data = dat)

test_that("Get do_boot results", {
  fit_boot <- sem(mod,
                  data = dat,
                  se = "boot",
                  bootstrap = 50,
                  iseed = 234567)

  out_boot <- lav_betaselect(fit_boot,
                             standardized = TRUE,
                             not_to_standardize = c("dv", "med", "cov2"),
                             progress = FALSE,
                             std_se = "bootstrap")

  out_boot_lav <- lav_betaselect(fit,
                                 standardized = TRUE,
                                 not_to_standardize = c("dv", "med", "cov2"),
                                 progress = FALSE,
                                 std_se = "bootstrap",
                                 bootstrap = 50,
                                 iseed = 234567)
  a <- round(out_boot_lav[1:4, "std.p.se"], 2)
  b <- round(out_boot[1:4, "se"], 2)
  i <- !is.na(a)
  expect_equal(a[i],
               b[i],
               tolerance = 1e-3,
               ignore_attr = TRUE)

})