To reduce outliers and simplify the calculation burden, we first retained 4677 prognostic-associated mRNAs selected by the univariate Cox proportional hazards regression (P < 0.05). A weighted correlation network analysis (WGCNA) was further performed on prognostic-associated mRNAs to construct scale-free gene co-expression networks, by using the R package “WGCNA”, with min-ModuleSize of 20, mergeCutHeight of 0.25, and unsigned topological overlap measure (TOM) [25]. An appropriate soft-threshold power was selected according to standard scale-free distribution (scale independence and mean connectivity). The module eigengenes (MEs) based on the first principal component were calculated for each module to measure the correlation between modules and clinical information including immunotype. The experimental protein-protein interactions (PPIs) data for genes within the immunotype-related module were retrieved from the STRING database, by specifying the interaction score of more than 0.4 ( The network was visualized within Cytoscape software (version 3.5.1.), and the hub genes were selected based on the Matthews correlation coefficient (MCC) score by using “cytoHubba” plug.

Note: The content above has been extracted from a research article, so it may not display correctly.

Please log in to submit your questions online.
Your question will be posted on the Bio-101 website. We will send your questions to the authors of this protocol and Bio-protocol community members who are experienced with this method. you will be informed using the email address associated with your Bio-protocol account.

We use cookies on this site to enhance your user experience. By using our website, you are agreeing to allow the storage of cookies on your computer.