A simulated dataset.
# \donttest{
syntax <- "
X =~ x1 + x2 + x3
Z =~ z1 + z2 + z3
Y =~ y1 + y2 + y3
W =~ w1 + w2 + w3
Y ~ X + Z + (1 + X + Z | cluster)
W ~ X + Z + (1 + X + Z | cluster)
"
fit <- pls(syntax, data = randomSlopes)
fit
#> plssem (0.1.4) ended normally after 77 iterations
#> lhs op rhs label est se z pvalue ci.lower ci.upper
#> 1 X =~ x1 0.913 NA NA NA NA NA
#> 2 X =~ x2 0.843 NA NA NA NA NA
#> 3 X =~ x3 0.877 NA NA NA NA NA
#> 4 Z =~ z1 0.915 NA NA NA NA NA
#> 5 Z =~ z2 0.828 NA NA NA NA NA
#> 6 Z =~ z3 0.877 NA NA NA NA NA
#> 7 Y =~ y1 0.910 NA NA NA NA NA
#> 8 Y =~ y2 0.850 NA NA NA NA NA
#> 9 Y =~ y3 0.874 NA NA NA NA NA
#> 10 W =~ w1 0.908 NA NA NA NA NA
#> 11 W =~ w2 0.838 NA NA NA NA NA
#> 12 W =~ w3 0.882 NA NA NA NA NA
#> 13 Y ~ X 0.291 NA NA NA NA NA
#> 14 Y ~ Z 0.456 NA NA NA NA NA
#> 15 W ~ X 0.382 NA NA NA NA NA
#> 16 W ~ Z 0.250 NA NA NA NA NA
#> 17 X ~~ X 1.000 NA NA NA NA NA
#> 18 X ~~ Z 0.174 NA NA NA NA NA
#> 19 Z ~~ Z 1.000 NA NA NA NA NA
#> 20 Y ~~ Y 0.662 NA NA NA NA NA
#> 21 W ~~ W 0.759 NA NA NA NA NA
#> 22 x1 ~~ x1 0.167 NA NA NA NA NA
#> 23 x2 ~~ x2 0.290 NA NA NA NA NA
#> 24 x3 ~~ x3 0.231 NA NA NA NA NA
#> 25 z1 ~~ z1 0.163 NA NA NA NA NA
#> 26 z2 ~~ z2 0.314 NA NA NA NA NA
#> 27 z3 ~~ z3 0.230 NA NA NA NA NA
#> 28 y1 ~~ y1 0.172 NA NA NA NA NA
#> 29 y2 ~~ y2 0.278 NA NA NA NA NA
#> 30 y3 ~~ y3 0.236 NA NA NA NA NA
#> 31 w1 ~~ w1 0.175 NA NA NA NA NA
#> 32 w2 ~~ w2 0.297 NA NA NA NA NA
#> 33 w3 ~~ w3 0.222 NA NA NA NA NA
#> 34 Y~1 ~~ Y~1 0.085 NA NA NA NA NA
#> 35 Y~X ~~ Y~1 -0.004 NA NA NA NA NA
#> 36 Y~Z ~~ Y~1 -0.027 NA NA NA NA NA
#> 37 Y~X ~~ Y~X 0.018 NA NA NA NA NA
#> 38 Y~Z ~~ Y~X 0.011 NA NA NA NA NA
#> 39 Y~Z ~~ Y~Z 0.107 NA NA NA NA NA
#> 40 W~1 ~~ W~1 0.058 NA NA NA NA NA
#> 41 W~X ~~ W~1 0.003 NA NA NA NA NA
#> 42 W~Z ~~ W~1 0.011 NA NA NA NA NA
#> 43 W~X ~~ W~X 0.096 NA NA NA NA NA
#> 44 W~Z ~~ W~X 0.010 NA NA NA NA NA
#> 45 W~Z ~~ W~Z 0.151 NA NA NA NA NA
# }