Now it integrates over effect parameters, based on a pmf that is fed into the BayesNet object. The default is a beta distribution with alpha=beta=2.