library("metaSEM") my.n <- c(3865,129,321,172,234,408,249,301,1659,294,16337,1845,4240,3802,99,834,281,1105,2042,10380,1270,16230,9602,271,1344,511,112,5862,994,548,432,1314,10518,896,1088,621,3887,2353,391,253,386,215,693,1775,345,191,564,196,2002,1038,227,4482,16543,15059,108,1675,940,2055,2157,38991,1536,9055,246,8933,876,162,2823,2236,1160,4280,21433,20683,783,39701,1939,16797,11011,5431,727,8239,11017,1859,2620,339,402,1930,177,2199,6973,87,519,547,764,229,1610,33997,15865,400,64436,39,9193,4465,8415,6208,210,297,3782,1324,218,183,183,183,183,1190,300,13617) my.vec <- readStackVec("data.dat") my.vec <- lapply(my.vec,function(x) {dimnames(x) <- list(c("OB3","OB1","WE","IE","FK","EE","NVI"), c("OB3","OB1","WE","IE","FK","EE","NVI")) x}) #Specify the Model’s A and S Matrices f<-c("0.2*index2OB3",0, "0.2*i_WE", "0.2*i_IE", "0.2*i_FK", "0.2*i_EE","0.2*i_NVI") F<-matrix(f,ncol = 1,nrow = 7) F A1 <- matrix(c(0, 0, 0, 0, 0, 0, 0, "0.2*OB32OB1", 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), nrow=7, ncol=7, byrow=TRUE) A1 A<-rbind(cbind(A1, F), matrix (c(0,0,0 , 0,0,0, 0,0), ncol=8, nrow=1)) A dimnames(A)<-list(c("OB3","OB1","WE","IE","FK","EE","NVI","index"), c("OB3","OB1","WE","IE","FK","EE","NVI","index")) A phi<-matrix(c(1),nrow = 1,ncol = 1) psi<-matrix(c("0.3*e1",0,0,0,0,0,0, 0, "0.3*e2",0,0,0,0,0, 0,0, "0.3*e3",0,0,0,0, 0,0,0, "0.3*e4",0,0,0, 0,0,0,0, "0.3*e5",0,0, 0,0,0,0,0,"0.3*e6",0, 0,0,0,0,0,0,"0.3*e7"), nrow=7, ncol=7, byrow=TRUE) S<-bdiagMat(list(psi,phi)) dimnames(S)<- list(c("OB3","OB1","WE","IE","FK","EE","NVI","index"), c("OB3","OB1","WE","IE","FK","EE","NVI","index")) S FM<-create.Fmatrix(c(1,1,1,1,1,1,1,0), as.mxMatrix = FALSE) dimnames(FM)<- list(c("OB3","OB1","WE","IE","FK","EE","NVI"), c("OB3","OB1","WE","IE","FK","EE","NVI","index")) FM #Fit the Model Random1 <- tssem1(my.vec, my.n, method = "REM") summary(Random1) FullMediation<-tssem2(Random1, Amatrix = A,Smatrix = S,Fmatrix = FM, diag.constraints = FALSE) summary(FullMediation)