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?