library(modsem)
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)
'
m1 <-# Outer Model
X =~ x1 + x2 + x3
Y =~ y1 + y2 + y3
Z =~ z1 + z2 + z3
# Inner model
Y ~ X + Z + X:Z
'
modsem(m1, oneInt)
est <- extract_lavaan(est)
lav_object <- bootstrapLavaan(lav_object, R = 500) bootstrap <-