Joint estimation of linear non-Gaussian acyclic models