Bayesian learning of multiple directed networks from observational data. (23rd September 2020)