# NOT RUN {
data(Cheung00)
## Variable labels
labels <- colnames(Cheung00$data[[1]])
## Full mediation model
S <- create.mxMatrix(c("1",
".2*cov_att_sn", "1",
0, 0, ".2*e_bi",
0, 0, 0, ".2*e_beh"),
type="Symm", as.mxMatrix=FALSE, byrow=TRUE)
dimnames(S) <- list(labels, labels)
S
A <- matrix(c("0","0","0","0",
"0","0","0","0",
".2*att2bi", ".2*sn2bi", "0", "0",
"0", "0", ".2*bi2beh", "0"),
byrow=TRUE, 4, 4)
dimnames(A) <- list(labels, labels)
A
#### Random-effects model
## Stage 1 analysis
random_1 <- tssem1(Cheung00$data, Cheung00$n, method="REM", RE.type="Symm", acov="individual")
summary(random_1)
## Stage 2 analysis
random_2 <- tssem2(random_1, Amatrix=A, Smatrix=S, intervals.type="LB", diag.constraints=TRUE)
summary(random_2)
## Load the library
library("semPlot")
my.plot1 <- meta2semPlot(random_2)
semPaths(my.plot1, whatLabels="est", col="yellow")
# }
Run the code above in your browser using DataLab