Discovering high-probable structural features has not been studied so extensively. Madigan
and York (1995) propose a Markov chain Monte Carlo (MCMC) method in the space of network
structures. Friedman and Koller (2003) design a more efficient MCMC procedure in the space of
variable orders. These algorithms output approximate posterior probabilities of structural features.
The approximation quality is not guaranteed in finite runs.