Getting Pvalue=0.00 for sem model

15 views
Skip to first unread message

zillur rahman

unread,
Sep 26, 2024, 5:49:53 PMSep 26
to lavaan
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)

Yago Luksevicius de Moraes

unread,
Sep 27, 2024, 9:55:03 AMSep 27
to lavaan
Hi,

What's the problem?
It just means that b1-b4 are statistically significant with p < 0.01, but a1-a4 are not.

Best,
Yago

Reply all
Reply to author
Forward
0 new messages