Bayesian learning of multiple directed networks from observational data