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("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf"),
                     c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf"))
x})


#Specify the Model’s A and S Matrices


A <- matrix(c(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.2*TTOres2LicRev", "0.2*UNIres2LicRev", "0.2*REGres2LicRev", "0.2*uSize2LicRev", 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=8, ncol=8, byrow=TRUE)
A
dimnames(A)<-list(c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf"), c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf"))
A
S<-matrix(c(1,0,0,0,0,0,0,0,
            0, 1,0,0,0,0,0,0,
            0,0, 1,0,0,0,0,0,
            0,0,0, 1,0,0,0,0,
            0,0,0,0, "0.3*e5",0,0,0,
            0,0,0,0,0, "0.3*e6",0,0,
            0,0,0,0,0,0,"0.3*e7",0,
            0,0,0,0,0,0,0,"0.3*e8"),
          nrow=8, ncol=8, byrow=TRUE)
dimnames(S)<- list(c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf"), c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf"))
S
FM<-create.Fmatrix(c(1,1,1,1,1,1,1,1), as.mxMatrix = FALSE)
dimnames(FM)<- list(c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf"),  c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf"))
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 = TRUE)

summary(FullMediation)


#Specify the Model’s A and S Matrices
f<-c("0.2*TTOonRES","0.2*UNIonRES", 0, "0", "0", "0","0","0")
F<-matrix(f,ncol = 1,nrow = 8)
F

A1 <- matrix(c(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.2*REGres2LicRev", "0.2*uSize2LicRev", 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=8, ncol=8, byrow=TRUE)
A1
A<-rbind(cbind(A1, F), matrix (c(0,0,0 , 0,"0.2*Resources2LicRev",0, 0,0,0), ncol=9, nrow=1))
A
dimnames(A)<-list(c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf","Resources"), c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf","Resources"))
A
phi<-matrix(c(1),nrow = 1,ncol = 1)
psi<-matrix(c("0.3*e1",0,0,0,0,0,0,0,
              0, "0.3*e2",0,0,0,0,0,0,
              0,0, 1,0,0,0,0,0,
              0,0,0, 1,0,0,0,0,
              0,0,0,0, "0.3*e5",0,0,0,
              0,0,0,0,0, "0.3*e6",0,0,
              0,0,0,0,0,0,"0.3*e7",0,
              0,0,0,0,0,0,0,"0.3*e8"),
            nrow=8, ncol=8, byrow=TRUE)
S<-bdiagMat(list(psi,phi))
dimnames(S)<- list(c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf","Resources"), c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf","Resources"))
S
FM<-create.Fmatrix(c(1,1,1,1,1,1,1,1,0), as.mxMatrix = FALSE)
dimnames(FM)<- list(c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf"),  c("TTO","UNI","REG","uSize","LicRev","LicNum","Spin","Perf","Resources"))
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 = TRUE)

summary(FullMediation)