The construction of the p1-model for a directed network is described in an Appendix Additional file 1: Appendix I. For the gene-gene relationship network with undirected edges, the description of the p1-model can be simplified by using only two Bernoulli variables Yij0 and Yij1 instead of four as follows: Yijk=1ifuij=k,0otherwise