SEM-mediation analysis: different results in JASP and R
Currently, I am trying to perform a mediation analysis using structural equation modeling in JASP. I wanted to check these results in R by copying the model from the Lavaan syntax option and testing this using the 'sem' function from Lavaan. However, I do not get the same results as in JASP when I conduct this mediation analysis in R this way. Because my p values differ in such a way that some are insignificant in R while they are significant in JASP, I am desperately trying to find out where this discrepancy comes from. Specifically, I used the bootstrap method with 1000 replications and the bias-corrected percentile type. Everything else was kept at the default.
What I entered in R was the following:
sem <- lavaan::sem(fullmodel, data = dataset, se = "bootstrap", bootstrap = 1000)
se = TRUE, zstat = TRUE, pvalue = TRUE, ci = TRUE,
standardized = TRUE,
fmi = FALSE, level = 0.95, boot.ci.type = "bca.simple")
As said before "fullmodel" was simply copy-pasted from what was printed as the Lavaan syntax in JASP.
Both direct, indirect and total effects give me different results.
Is there anything I am missing in these functions?
Thanks in advance for the help!