# 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)
## Display the model
plot(random_2, what="path")
## Display the model with the parameter estimates
plot(random_2, color="yellow")
## Load the library
library("semPlot")
# }
Run the code above in your browser using DataLab