Added methods sample_dense_matrix, construct_bn, compute_cond_probs_for_chrons_2009 and hc_for_optimal_treatment