Phi1<-matrix(c(1,"0.3*cor1","0.3*cor2","0.3*cor1",1,"0.3*cor3","0.3*cor2","0.3*cor3",1),ncol=3, nrow=3) Psi1<-Diag(c("0.2*e1","0.2*e2","0.2*e3","0.2*e4","0.2*e5","0.2*e6","0.2*e7","0.2*e8")) S2<-bdiagMat(list(Psi1, Phi1)) LambdaX<-matrix(c(".3*X9_X1",".3*X9_X2",".3*x9*X3",rep(0,3),".3*X9_X7",rep(0,5),".3*x10_x5",0,".3*x10_X7",".3*x10_X8",0, ".3*x11_X2",0,".3*x11_X4",0,".3*X11_X6",0,0),ncol=3,nrow=8) AX<-rbind(cbind(matrix(0,ncol=8,nrow=8),LambdaX),matrix(0,ncol=11,nrow=3)) FX<-create.Fmatrix(c(1,1,1,1,1,1,1,1,0,0,0), as.mxMatrix=FALSE) randomX<-tssem2(randomB5Full,Amatrix=AX, Smatrix=S2, Fmatrix=FX, diag.constraints=FALSE, intervals="LB") summary(randomX) myplotX<-meta2semPlot(randomX) semPaths(myplotX, whatLabels="est", nCharEdges=10, nCharNodes=10, color="blue")