I have a dataset of 150 samples of 60 binary variables. The pcalg binCItest results are not impressive and i can not find a straightforward way to do this using bnlearn or catnet. I appreciate if you let me know your suggestion, especially if accompanied by a code snippet. Thanks