The weighted gene co-expression network analysis (WGCNA) method constructs co-expressed gene modules based on the linear similarity between paired gene expressions. Linear correlations are the main form of similarity between genes, however, nonlinear correlations still existed and had always been ignored. We proposed a modified network analysis method, WGCNA-P + M, which combines Pearson's correlation coefficient and the maximum information coefficient (MIC) as the similarity measures to assess the linear and nonlinear correlations between genes, respectively. Taking two real datasets, GSE44861 and liver hepatocellular carcinoma (TCGA-LIHC), as examples, we compared the gene modules constructed by WGCNA-P + M and WGCNA from four perspectives: the "Usefulness" score, GO enrichment analysis on genes in the gray module, prediction performance of the top hub gene, survival analysis and literature reports on different hub genes. The results showed that the modules obtained by WGCNA-P + M are more biological meaningful, the hub genes obtained from WGCNA-P + M have more potential cancer genes.