## 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")
# 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")
# ## End(Not run)
Run the code above in your browser using DataLab