Hi,
I have the model below. It is running fine but I am getting pvalue=0.00 for b1, b2, b3, and b4. On the other hand a1, a2, a3, a4 has no such problem.
Any Help?
mediation_model_bip <- '
#direct effects
PhysAct ~ a1 * BIP_PRS + Sex + Age + PC1 + PC2 + PC3 + PC4 + PC5 + PC6 + PC7 + PC8 + PC9 + PC10
Screen ~ a2 * BIP_PRS + Sex + Age + PC1 + PC2 + PC3 + PC4 + PC5 + PC6 + PC7 + PC8 + PC9 + PC10
Healthy_Food ~ a3 * BIP_PRS + Sex + Age + PC1 + PC2 + PC3 + PC4 + PC5 + PC6 + PC7 + PC8 + PC9 + PC10
Meat ~ a4 * BIP_PRS + Sex + Age + PC1 + PC2 + PC3 + PC4 + PC5 + PC6 + PC7 + PC8 + PC9 + PC10
BMI ~ c * BIP_PRS + b1 * PhysAct + b2 * Screen + b3 * Healthy_Food + b4 * Meat
#indirect effect (a * b)
indirect1 := a1 * b1
indirect2 := a2 * b2
indirect3 := a3 * b3
indirect4 := a4 * b4
#total effect (c + indirect)
total := c + indirect1 + indirect2 + indirect3 + indirect4
PhysAct ~~ Screen
PhysAct ~~ Healthy_Food
PhysAct ~~ Meat
Screen ~~ Healthy_Food
Screen ~~ Meat
Healthy_Food ~~ Meat
'
mediation_results_bip <- sem(mediation_model_bip, data = df2, verbose = TRUE, ordered = 'PhysAct', se = 'bootstrap',estimator = 'DWLS')
summary(mediation_results_bip, standardized = TRUE, fit.measures = TRUE)