Attachment | Size |
---|---|

path diagram | 26.42 KB |

Dear Mike,

A colleague and I studied the factor structure of a scale using the metaSEM package. The reviewers asked us for residual covariance matrice of structural model we construct. I couldn't see a command for this in the metaSEM package. There is vcov but this command gives sampling covariance matrices. Is there a way to see residual covariance matrices?

Sincerely

The script that we used to construct our model is ad follow:

library(metaSEM)

library(lavaan)

sdeğ <- 6

değad <- c("m1","m2","m3","m4","m5","m6")

etiket <- list(değad,değad)

korver <- list()

for (i in 1:nrow(veri)){

korver[[i]] <- vec2symMat(as.matrix(veri[i,3:17]),diag = FALSE)

dimnames(korver[[i]]) <- etiket

}

korver

REM1 <- tssem1(Cov=korver, n=veri$N, method="REM", RE.type="Diag")

summary(REM1)

model<-"

Faktor1=~M1+M2+M3

Faktor2=~M4+M5+M6

Faktor1~~Faktor2"

RAM <- lavaan2RAM(model, obs.variables=c("M1","M2","M3","M4","M5", "M6"),

A.notation="on", S.notation="with")

REM2 <- tssem2(REM1, RAM=RAM, intervals="LB")

summary(REM2)

Call:

wls(Cov = pooledS, aCov = aCov, n = tssem1.obj$total.n, RAM = RAM,

Amatrix = Amatrix, Smatrix = Smatrix, Fmatrix = Fmatrix,

diag.constraints = diag.constraints, cor.analysis = cor.analysis,

intervals.type = intervals.type, mx.algebras = mx.algebras,

model.name = model.name, suppressWarnings = suppressWarnings,

silent = silent, run = run)

95% confidence intervals: Likelihood-based statistic

Coefficients:

Estimate Std.Error lbound ubound z value Pr(>|z|)

M1onFaktor1 0.75023 NA 0.69404 0.80902 NA NA

M2onFaktor1 0.81600 NA 0.75594 0.87878 NA NA

M3onFaktor1 0.72925 NA 0.67467 0.78664 NA NA

M4onFaktor2 0.79731 NA 0.72519 0.87357 NA NA

M5onFaktor2 0.80009 NA 0.72521 0.87772 NA NA

M6onFaktor2 0.79901 NA 0.72891 0.87430 NA NA

Faktor1withFaktor2 0.17314 NA 0.13774 0.20912 NA NA

Goodness-of-fit indices:

Value

Sample size 5769.0000

Chi-square of target model 11.0801

DF of target model 8.0000

p value of target model 0.1972

Number of constraints imposed on "Smatrix" 0.0000

DF manually adjusted 0.0000

Chi-square of independence model 2034.0182

DF of independence model 15.0000

RMSEA 0.0082

RMSEA lower 95% CI 0.0000

RMSEA upper 95% CI 0.0186

SRMR 0.0193

TLI 0.9971

CFI 0.9985

AIC -4.9199

BIC -58.2020

OpenMx status1: 0 ("0" or "1": The optimization is considered fine.

Other values indicate problems.)

plot(REM2)

The attached approach should work.

Thank you very much. it worked!

Dear Cheung,

We want to investigate whether and how the regression coefficients of the structural models differ across studies with sample of children or adult. In other words, is it possible to do Random effect model subgroup analysis?

Yes, but not straightforward.

Instead of running the stage 2 model with tssem2(...), you may export the model and data to a mxmodel with tssem2(..., run=FALSE).

You may then use mxFitFunctionMultigroup() in OpenMx to conduct a multiple-group analysis directly.

Dear Mike,

Thank you very much for your replay. I'm pretty new to openmx. How to export my model and data to mxModel(). If we have stage 1 result (REM), RAM model and data file (ls), how can we write the mxModel and conduct tssem2.

if you don't mind can you write me the script or send me an example?

Best regards.

It can be complicated.

It is easier to apply the OSMASEM approach by using a dummy code of children or adults as a predictor. There are a couple of examples, e.g.,

https://github.com/mikewlcheung/code-in-articles/tree/master/Schutte%20Keng%20and%20Cheung%202021

https://github.com/mikewlcheung/code-in-articles/tree/master/Jak%20and%20Cheung%202020

Thanks, I'm reviewing both the codes and your article "Emotional Intelligence Mediates the Connection Between Mindfulness and Gratitude".