-
Notifications
You must be signed in to change notification settings - Fork 0
/
Scripts_Details.rtf
27 lines (26 loc) · 2.51 KB
/
Scripts_Details.rtf
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
{\rtf1\ansi\ansicpg1252\cocoartf1265\cocoasubrtf210
{\fonttbl\f0\fswiss\fcharset0 Helvetica;}
{\colortbl;\red255\green255\blue255;}
\paperw11900\paperh16840\margl1440\margr1440\vieww10800\viewh8400\viewkind0
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural
\f0\fs24 \cf0 \
The Scripts are in the following order: \
\
1: QCfilter.R. This script was used in order to perform a preprocessing of gene expression microarray data downloaded from public repositories and \
prior to any normalization.\
\
2: Mutual_Exclusivity_Script.R. This script was used in order to run the mutual exclusivity analysis and in particular to calculate the \
significance of the enrichment of the mutually exclusive mutated genes in the CM-GRNs with respect to randomly generated networks. The script \
takes in input two files: "CM_Nets_TCGA_mutational_data_1.txt" that contains three columns: the list of the CM-GRNs in the first column, the number of \
mutually exclusive mutated genes in the network in the second column and the total set of patients that were considered (TCGA mutational data) in the third column. \
The second input file is called: "RANDOM_Nets_TCGA_mutational_data_1.txt". This file contains the same type of information of the previous one \
but it refers to randomly generated networks. The "*_1*" in both file names refers to the single gene mutual exclusivity. The mutual exclusivity \
was also investigated between groups of genes (e.g. couples of mutually exclusive mutated genes) although in the manuscript only single genes \
were considered. \
\
3: Parallel_Aracne script. This script was used to speed up and parallelise the ARACNE calculations on the computer cluster because of the high number of HUBS (1,516 genes against the whole transcriptome) to be tested in the reverse engineering analysis performed on: breast cancer and Glioblastoma, Lung (LUSC, LUAD) and Ovarian Cancer \
\
\
4: Concordance analysis.R. This script was used to calculate the Activation Score (AS) of the networks relative to the expression profile of breast cancer patients form the Metabric cohort. In particular, the script transforms the original continuous data matrix in a discrete 0-1 indexed matrix in order to calculate the AS i.e. the number of genes of each network having the same or the opposite expression pattern with respect to the transcriptional profile of the hub gene. The hub-gene neighbourhood of each GRNs, was matched with the genes of the Metabric cohort after the expression matrix transformation.\
\
}