A proper understanding of complex biological networks facilitates a better perception of those diseases that plague systems and efficient production of drug targets, which is one of the major research questions under the personalized medicine. However, the description of these complexities is challenging due to the associated continuous, high-dimensional, correlated and very sparse data. The Copula Gaussian Graphical Model (CGGM), which is based on the representation of the multivariate normal distribution via marginal and copula terms, is one of the successful modeling approaches to presenting such types of problematic datasets. This study shows its novelty by using CGGM in modeling the steady-state activation of biological networks and making inference of the model parameters under the Bayesian setting. In this regard, the Reversible Jump Markov Chain Monte Carlo (RJMCMC) algorithm is suggested in order to estimate the plausible interactions (conditional dependence) between the systems' elements, which are proteins or genes. Furthermore, the open-source R codes of RJMCMC are generated for CGGM in different dimensional networks. In this regard, real datasets are applied, and the accuracy of estimates via F-measure is evaluated. From the results, it is observed that CGGM with RJMCMC is successful in presenting real and complex systems with higher accuracy. (C) 2019 Sharif University of Technology. All rights reserved.