Changes between Version 12 and Version 13 of SOPs/go_annotation
- Timestamp:
- 09/28/20 09:04:08 (4 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/go_annotation
v12 v13 13 13 [[http://www.broadinstitute.org/gsea/index.jsp|Broad GSEA]] 14 14 15 [[http://www.broadinstitute.org/cancer/software/gsea/wiki/index.php/Main_Page | GSEA Wiki]] 15 [[https://www.gsea-msigdb.org/gsea/login.jsp|Download the GSEA software and additional resources to analyze, annotate and interpret enrichment results.]] 16 17 [[https://www.gsea-msigdb.org/gsea/msigdb/index.jsp|Explore the Molecular Signatures Database (MSigDB), ]]a collection of annotated gene sets for use with GSEA software. 18 19 [[http://software.broadinstitute.org/cancer/software/gsea/wiki/index.php/Main_Page | GSEA and MSigDB documentation]] 20 21 [[http://software.broadinstitute.org/cancer/software/gsea/wiki/index.php/Using_RNA-seq_Datasets_with_GSEA | Guidelines for using RNA-seq datasets with GSEA]] 16 22 17 23 ==== Ranked List ==== 18 1. Create a two column file with gene names as first column and numeric values for second column (eg. weight, p-value, etc), does not need to be sorted.24 1. Create a two column file with gene names as first column and numeric values for second column (eg. log FC), does not need to be sorted and it should have extension ".rnk". 19 25 * Assigning weights: There is no standard way to assign weights, however, it should reflect some logical order. GSEA uses the correlation between expression and phenotype to assign weights, if the list is not pre-ordered or ranked. A similar scheme can be used to rank genes by log2 ratio, t-statistic, or a scoring scheme that takes into account both log ratio and p-value. 20 26 * If a gene list is not unique, duplicate genes can be given a //shared// weight, for eg. if a gene occurs four times in the list it is given a weight of 0.25, if it is unique a weight of 1 is given. … … 30 36 GSEA will rank the genes 31 37 32 1. Create necessary files in correct format for expression, phenotype and chip annotation ( see GSEA wiki)38 1. Create necessary files in correct format for expression, phenotype and chip annotation ([[http://software.broadinstitute.org/cancer/software/gsea/wiki/index.php/Using_RNA-seq_Datasets_with_GSEA | see GSEA wiki]]) 33 39 1. Use MSigDB for gene sets or create custom gene sets in correct format 34 40 1. Run GSEA, use default options to start