Dear Dr mike :

I have a question about MASEM. My example is to study the effect of social support on depression.I found four variables : subject support (zz), objective support (kz) , the using of support (zl), depression (yy). Among these , zz, kz, zl were loaded under a higher-order factor ,called "Alpha", that represents social support. The following is my R code :

a<-c(1,0.428,0.41,-0.327,0.428,1,0.292,-0.199,0.41,0.292,1,-0.198,-0.327,-0.199,-0.198,1)

m1<-matrix(a,nrow=4,ncol=4,dimnames=list(rowname,colname))

rowname=c("zz","kz","zl","yy")

rowname

colname=c("zz","kz","zl","yy")

colname

m1

b<-c(1,0.482,0.239,-0.183,0.482,1,0.305,-0.229,0.239,0.305,1,-0.256,-0.183,-0.229,-0.256,1)

m2<-matrix(a,nrow=4,ncol=4,dimnames=list(rowname,colname))

rowname=c("zz","kz","zl","yy")

rowname

colname=c("zz","kz","zl","yy")

colname

m2

c<-c(1,0.474,0.360,-0.044,0.474,1,0.42,-0.113,0.360,0.42,1,-0.185,-0.044,-0.113,-0.185,1)

m3<-matrix(c,nrow=4,ncol=4,dimnames=list(rowname,colname))

rowname=c("zz","kz","zl","yy")

rowname

colname=c("zz","kz","zl","yy")

colname

m3

d<-c(1,0.330,0.294,-0.246,0.330,1,0.261,-0.264,0.294,0.261,1,-0.217,-0.246,-0.264,-0.217,1)

m4<-matrix(d,nrow=4,ncol=4,dimnames=list(rowname,colname))

rowname=c("zz","kz","zl","yy")

rowname

colname=c("zz","kz","zl","yy")

colname

m4

e<-c(1,0.25,0.24,-0.21,0.25,1,0.15,-0.19,0.24,0.15,1,-0.23,-0.21,-0.19,-0.23,1)

m5<-matrix(e,nrow=4,ncol=4,dimnames=list(rowname,colname))

rowname=c("zz","kz","zl","yy")

rowname

colname=c("zz","kz","zl","yy")

colname

m5

l1<-list(m1,m2,m3,m4,m5)

l2<-list(A=m1,B=m2,C=m3,D=m4,E=m5)

names(l2)

l2

n<-c(211,123,459,698,310)

n

#### Load the metaSEM library for TSSEM

library(metaSEM)

#### Inspect the data for inspection

head(l2)

head(n)

#### Fixed-effects model: Stage 1 analysis

fixed1 <- tssem1(my.df=l2, n=n, method="FEM")

summary(fixed1)

coef(fixed1)

I didn't know how to specify, so I could't go on . I look forward to reading your responses. Thank you!

regards,

xiarizhiyue0727