R语言GSEA分析(三)
GAEA
df_all_sort <- df_all[order(df_all$logFC, decreasing = T),]#先按照logFC降序排序 gene_fc = df_all_sort$logFC #把foldchange按照从大到小提取出来 head(gene_fc) names(gene_fc) <- df_all_sort$ENTREZID #给上面提取的foldchange加上对应上ENTREZID head(gene_fc) > head(gene_fc) [1] 41.99218 35.00852 22.78417 16.82801 16.64222 15.33221 > head(gene_fc) 972 126868 10984 201853 378938 3514 41.99218 35.00852 22.78417 16.82801 16.64222 15.33221
赞 (0)