Maximum likelihood fitting of acyclic directed mixed graphs to binary   data
  Acyclic directed mixed graphs, also known as semi-Markov models represent the conditional independence structure induced on an observed margin by a DAG model with latent variables. In this paper we present the first method for fitting these models to binary data using maximum likelihood estimation. 
Robin J. Evans, Thomas S. Richardson
stat.ME cs.AI