為什么GEO2R/GEOquery的結(jié)果可能是錯的？

在準備GEO/TCGA培訓(xùn)時，因為數(shù)據(jù)部分來自NCBI GEO數(shù)據(jù)庫，就系統(tǒng)梳理了GEO數(shù)據(jù)庫的結(jié)構(gòu)和存儲方式，順便也看了GEO2R的使用。然后就發(fā)現(xiàn)了一個問題，可能直接使用GEO2R的分析結(jié)果是錯誤的。

原因見下圖 （官網(wǎng)介紹截圖）

而GEO2R的R代碼使用的是GEOquery獲取GSEMatrix，所以如果我們自己寫代碼這樣獲取非原始數(shù)據(jù)時，也有可能得到的數(shù)據(jù)是不可比的。

所以，還是建議使用工具或命令之前，好好讀一下幫助文檔，做到心中有數(shù)，知其所以然。如果能基于原始數(shù)據(jù)進行分析，可以做更多的質(zhì)控和更深入的比較。如果不能，分析之前看下數(shù)據(jù)的分布是否均一 (median-centered)。

所以生信寶典的原理類文章，值得好好閱讀。

典型醫(yī)學(xué)設(shè)計實驗GEO數(shù)據(jù)分析 (step-by-step) - Limma差異分析、火山圖、功能富集

典型醫(yī)學(xué)設(shè)計實驗GEO數(shù)據(jù)分析 (step-by-step) - 數(shù)據(jù)獲取到標準化

GO、GSEA富集分析一網(wǎng)打進

GSEA富集分析 - 界面操作

WGCNA分析，簡單全面的最新教程

PCA主成分分析實戰(zhàn)和可視化 附R代碼和測試數(shù)據(jù)

一文掌握GSEA，超詳細教程

從課件中截取GEO簡介部分分享如下 （只占課件很小一部分）：

總結(jié)

以上是生活随笔為你收集整理的为什么GEO2R/GEOquery的结果可能是错的？的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。