The adoptability of the heart to external and internal stimuli is reflected by heart rate variability (HRV). Reduced HRV can be a predictor of post-infarction mortality. In this study, we propose an automated system to predict and diagnose congestive heart failure using short-term heart rate variability analysis. Based on the nonlinear, nonstationary, and highly complex dynamics of congestive heart failure, we extracted multimodal features to capture the temporal, spectral, and complex dynamics. Recently, the Bayesian inference approach has been recognized as an attractive option for the deeper analysis of static features, in order to perform a comprehensive analysis of extracted nodes (features). We computed the gray level co-occurrence (GLCM) features from congestive heart failure signals and then ranked them based on ROC methods. This study focused on utilizing the dissimilarity feature, which is ranked as highly important, as a target node for the empirical analysis of dynamic profiling and optimization, in order to explain the nonlinear dynamics of GLCM features extracted from heart failure signals, and distinguishing CHF from NSR. We applied Bayesian inference and Pearson’s correlation (PC). The association, in terms of node force and mapping, was computed. The higher-ranking target node was used to compute the posterior probability, total effect, arc contribution, network profile, and compression. The highest value of ROC was obtained for dissimilarity, at 0.3589. Based on the information-gain algorithm, the highest strength of the relationship was obtained between nodes “dissimilarity” and “cluster performance” (1.0146), relative to mutual information (81.33%). Moreover, the highest relative binary significance was yielded for dissimilarity for 1/3rd (80.19%), 2/3rd (74.95%) and 3/3rd (100%). The results revealed that the proposed methodology can provide further in-depth insights for the early diagnosis and prognosis of congestive heart failure.