I have a philosophical question about causal relationships inference. Suppose we have UCBAdmission data, where different gender and different department have different admission rates.
> d=as.data.frame(UCBAdmissions)
> d
Admit Gender Dept Freq
1 Admitted Male A 512
2 Rejected Male A 313
3 Admitted Female A 89
4 Rejected Female A 19
...
...
21 Admitted Male F 22
22 Rejected Male F 351
23 Admitted Female F 24
24 Rejected Female F 317
My question is can we do any causal relationships inference on this data? I would say both gender and department would case difference on admission, but how can I know it is gender cause department or the other way around?
In other words, how can I do Bayesian network structure learning on this data? Is it even possible?