using lavaan functions

library(modsem)
#> This is modsem (1.0.11). Please report any bugs!

If you’re using one of the product indicator approaches, you might want to use some lavaan functions on top of the estimated lavaan-object. To do so you can extract the lavaan-object using the extract_lavaan() function.

library(lavaan)
#> This is lavaan 0.6-19
#> lavaan is FREE software! Please report any bugs.

m1 <- '
# Outer Model
X =~ x1 + x2 + x3
Y =~ y1 + y2 + y3
Z =~ z1 + z2 + z3

# Inner model
Y ~ X + Z + X:Z
'

est <- modsem(m1, oneInt)
lav_object <- extract_lavaan(est)
bootstrap <- bootstrapLavaan(lav_object, R = 500)