fitAbn function

Fit an additive Bayesian network model