-
Notifications
You must be signed in to change notification settings - Fork 21
Description
Hi :)
Thank you for creating this great statistical test and package, really helpful!
I have a conceptual questions about evaluating the NCT results. I'm estimating two networks of depression symptoms post-intervention, one for a control group and one for the intervention group.
`network_model_T1_SCCM <- estimateNetwork(ana_data_T1_SCCM,
default = "EBICglasso",
tuning = 0.5)
network_model_T1_TAU <- estimateNetwork(ana_data_T1_TAU,
default = "EBICglasso",
tuning = 0.5)
NCT_SCCMTAUT1 <- NCT(network_model_T1_SCCM, network_model_T1_TAU, it=1000, test.edges = TRUE, abs=TRUE)
summary(NCT_SCCMTAUT1)
`
I get the following, non-significant results for network invariance and global strength invariance:
NETWORK INVARIANCE TEST
Test statistic M: 0.2550473
p-value 0.3326673
GLOBAL STRENGTH INVARIANCE TEST
Global strength per group: 3.663337 3.700252
Test statistic S: 0.036915
p-value 0.8591409
However, when I check out individual edges, I do actually find individual significant edge differences between the networks, e.g....
37 phq9_interest phq9_appetite 0.018981019 0.19195629
48 phq9_insomnia phq9_guilt 0.048951049 0.15244224
67 phq9_noenergy phq9_psychomotor 0.003996004 0.25504734
It was my understanding that a significant networ invariance test, being an omnibus test, indicates that there is at least one significant edge difference, i.e. I shouldn't be finding significant results in the edge invariance test.
What am I missing? Is this still a reasonable result or does this indicate mistakes in my analysis? Or did I misunderstand the network invariance test?
Thank you in advance!
Best,
Max