data(ECSImobi)
ecsi <- sempls(model=ECSImobi, data=mobi)
### Bootstrapping
set.seed(123)
ecsiBoot <- bootsempls(ecsi, nboot=200, start="ones", verbose=TRUE)
summary(ecsiBoot, type="perc", level=0.95)
## inspectation of bootstrap samples
parallelplot(ecsiBoot, subset=1:ncol(ecsiBoot$t), reflinesAt=0)
# only inspecting the path coefficients
parallelplot(ecsiBoot, pattern="beta", reflinesAt=c(0,1))
densityplot(ecsiBoot, pattern="beta")
# only inspecting the outer loadings
parallelplot(ecsiBoot, pattern="lam")
# only inspecting the outer loadings for Loyalty
parallelplot(ecsiBoot, pattern="lam7", type="perc", level=0.90,
main="Loyalty
200 bootstrapped outer loadings")
Run the code above in your browser using DataLab