 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (499)
- >考研 (51)
- >蟲(chóng)友互識(shí) (34)
- >導(dǎo)師招生 (15)
- >教師之家 (8)
- >考博 (8)
- >碩博家園 (7)
- >基金申請(qǐng) (6)
- >文獻(xiàn)求助 (5)
- >無(wú)機(jī)/物化 (3)
- >論文道賀祈福 (2)
- >文獻(xiàn)求助完結(jié)子版 (2)
- >找工作 (2)
- >公派出國(guó) (2)
- >論文投稿 (2)
- >復(fù)合材料 (1)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 20694 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
QQ894064647新蟲(chóng) (初入文壇)
|
[交流]
TCGA(癌癥和腫瘤基因圖譜)數(shù)據(jù)下載和處理(TCGA-Assembler) 已有2人參與
|
||
|
國(guó)政府發(fā)起的癌癥和腫瘤基因圖譜（Cancer Genome Atlas，TCGA）計(jì)劃，試圖通過(guò)應(yīng)用基因組分析技術(shù)，特別是采用大規(guī)模的基因組測(cè)序，將人類全部癌癥（近期目標(biāo)為50種包括亞型在內(nèi)的腫瘤）的基因組變異圖譜繪制出來(lái)，并進(jìn)行系統(tǒng)分析，旨在找到所有致癌和抑癌基因的微小變異，了解癌細(xì)胞發(fā)生、發(fā)展的機(jī)制，在此基礎(chǔ)上取得新的診斷和治療方法，最后可以勾畫(huà)出整個(gè)新型“預(yù)防癌癥的策略”。 TCGA 使命:提高人們對(duì)癌癥發(fā)病分子基礎(chǔ)的科學(xué)認(rèn)識(shí)及提高我們?cè)\斷、治療和預(yù)防癌癥的能力 TCGA 目標(biāo):完成一套完整的與所有癌癥基因組改變相關(guān)的“圖譜”。 圖1.png TCGA數(shù)據(jù)源大部分都是公開(kāi)的，如何有效的進(jìn)行收集和預(yù)處理是一個(gè)頭疼的問(wèn)題。今天我們講解下怎么將TCGA的數(shù)據(jù)轉(zhuǎn)化成癌癥類型的二維數(shù)據(jù)矩陣（例如基因?yàn)閞ows，樣本為columns）。得到這個(gè)矩陣之后，后面的事情就好辦了，我們可以做差異表達(dá)，共表達(dá)網(wǎng)絡(luò)，生存分析等。今天我們主要講解如何下載TCGA的數(shù)據(jù)，大家對(duì)后續(xù)分析感興趣的話，可以在加“生物信息培訓(xùn)+視頻”裙，或者大家可以在掏寶搜索“生物信息視頻”，跟我們聯(lián)系。 我們開(kāi)始吧，我們可以使用TCGA-Assembler這軟件去下載TCGA的數(shù)據(jù)http://www.compgenome.org/TCGA-Assembler/。TCGA-Assembler不但可以很方便的下載數(shù)據(jù)，還能對(duì)數(shù)據(jù)進(jìn)行初始化處理，非常方便。下載完后，我們使用首先要安裝一些依賴包。通過(guò)下面的命令： install.packages(c("HGNChelper", "RCurl", "httr", "stringr", "digest", "bitops"  , dependencies=T)安裝完了依賴包，我們進(jìn)入剛才下載的TCGA-Assembler的目錄，使用setwd(C:/Users/cloud/Desktop/TCGA-Assembler)設(shè)置TCGA-Assembler的目錄為工作目錄，接下來(lái)，我們就可以下載數(shù)據(jù)了。我們需要下載什么數(shù)據(jù)，就選擇相應(yīng)的腳本。具體腳本如下： # Load module A functions. source("Module_A.r" ;# Download level-3 miRNA-seq data of six rectum adenocarcinoma (READ) samples miRNASeqRawData = DownloadmiRNASeqData(traverseResultFile = "./DirectoryTraverseResult_Jul-08-2014.rda", saveFolderName = "./QuickStartGuide_Results/RawData/", cancerType = "READ", assayPlatform = "miRNASeq", inputPatientIDs = c("TCGA-EI-6884-01", "TCGA-DC-5869-01", "TCGA-G5-6572-01", "TCGA-F5-6812-01", "TCGA-AF-2689-11", "TCGA-AF-2691-11" ); # Download level-3 DNA copy number data of six READ samples CNARawData = DownloadCNAData(traverseResultFile = "./DirectoryTraverseResult_Jul-08-2014.rda", saveFolderName = "./QuickStartGuide_Results/RawData/", cancerType = "READ", assayPlatform = "genome_wide_snp_6", inputPatientIDs = c("TCGA-EI-6884-01", "TCGA-DC-5869-01", "TCGA-G5-6572-01", "TCGA-F5-6812-01", "TCGA-AF-2692-10", "TCGA-AG-4021-10" ); # Download level-3 RNASeqV2 gene expression and exon expression data of six READ samples RNASeqRawData = DownloadRNASeqData(traverseResultFile = "./DirectoryTraverseResult_Jul-08-2014.rda", saveFolderName = "./QuickStartGuide_Results/RawData/", cancerType = "READ", assayPlatform = "RNASeqV2", dataType = c("rsem.genes.normalized_results", "exon_quantification" , inputPatientIDs = c("TCGA-EI-6884-01", "TCGA-DC-5869-01", "TCGA-G5-6572-01", "TCGA-F5-6812-01", "TCGA-AG-3732-11", "TCGA-AG-3742-11" ); # Download level-3 HumanMethylation27 data of six READ samples Methylation27RawData = DownloadMethylationData(traverseResultFile = "./DirectoryTraverseResult_Jul-08-2014.rda", saveFolderName = "./QuickStartGuide_Results/RawData/", cancerType = "READ", assayPlatform = "humanmethylation27", inputPatientIDs = c("TCGA-AG-3583-01", "TCGA-AG-A032-01", "TCGA-AF-2692-11", "TCGA-AG-4001-01", "TCGA-AG-3608-01", "TCGA-AG-3574-01" ); # Download level-3 HumanMethylation450 data of six READ samples Methylation450RawData = DownloadMethylationData(traverseResultFile = "./DirectoryTraverseResult_Jul-08-2014.rda", saveFolderName = "./QuickStartGuide_Results/RawData", cancerType = "READ", assayPlatform = "humanmethylation450", inputPatientIDs = c("TCGA-EI-6884-01", "TCGA-DC-5869-01", "TCGA-G5-6572-01", "TCGA-F5-6812-01", "TCGA-AG-A01W-11", "TCGA-AG-3731-11" ); # Download level-3 RPPA protein expression data of six READ samples RPPARawData = DownloadRPPAData(traverseResultFile = "./DirectoryTraverseResult_Jul-08-2014.rda", saveFolderName = "./QuickStartGuide_Results/RawData", cancerType = "READ", assayPlatform = "mda_rppa_core", inputPatientIDs = c("TCGA-EI-6884-01", "TCGA-DC-5869-01", "TCGA-G5-6572-01", "TCGA-F5-6812-01", "TCGA-AG-3582-01", "TCGA-AG-4001-01" ); # Download de-identified clinical information of READ patients DownloadClinicalData(traverseResultFile = "./DirectoryTraverseResult_Jul-08-2014.rda", saveFolderName = "./QuickStartGuide_Results/RawData", cancerType = "READ", clinicalDataType = c("patient", "drug", "follow_up" );運(yùn)行上面的腳本，我們就能得到我們想要的結(jié)果了，假如我們需要下載adenocarcinoma的miRNA數(shù)據(jù)，我們可以使用。下載完后，我們就得到了adenocarcinoma的矩陣了（基因?yàn)閞ows，樣本為columns）。 setwd(C:/Users/cloud/Desktop/TCGA-Assembler) source("Module_A.r" ;miRNASeqRawData = DownloadmiRNASeqData(traverseResultFile = "./DirectoryTraverseResult_Jul-08-2014.rda", saveFolderName = "./QuickStartGuide_Results/RawData/", cancerType = "READ", assayPlatform = "miRNASeq" ; |
回復(fù)此樓» 猜你喜歡
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
-
327求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 急求！核苷酸序列翻譯成氨基酸后，怎么回到對(duì)應(yīng)的核苷酸序列？
已經(jīng)有4人回復(fù)
- 請(qǐng)問(wèn)TCGA的level2數(shù)據(jù)是通過(guò)哪種方法標(biāo)準(zhǔn)化的？
已經(jīng)有5人回復(fù)
- 從網(wǎng)上下載了cel文件，但不知怎么下對(duì)應(yīng)的實(shí)驗(yàn)條件信息
已經(jīng)有3人回復(fù)
- 求氫氣分子的鍵長(zhǎng)數(shù)據(jù)
已經(jīng)有3人回復(fù)
- TCGA數(shù)據(jù)庫(kù)的使用交流
已經(jīng)有12人回復(fù)
- 有人做過(guò)或接觸過(guò)microRNA的paired-end測(cè)序數(shù)據(jù)分析嗎
已經(jīng)有18人回復(fù)
- 如何找到DNA甲基化的位點(diǎn)
已經(jīng)有3人回復(fù)
- The Cancer Genome Atlas (TCGA)數(shù)據(jù)庫(kù)的使用？
已經(jīng)有6人回復(fù)
- 求編程，求應(yīng)助
已經(jīng)有3人回復(fù)
- 求TCGA 數(shù)據(jù)庫(kù)的使用說(shuō)明，關(guān)于數(shù)據(jù)的。
已經(jīng)有7人回復(fù)
- 【視頻】【教程】PDB數(shù)據(jù)庫(kù)的使用——視頻教程
已經(jīng)有254人回復(fù)
- 有關(guān)PDB數(shù)據(jù)庫(kù)的使用問(wèn)題
已經(jīng)有11人回復(fù)
xingzhou823
木蟲(chóng) (正式寫(xiě)手)
五道杠
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 2391.6
- 散金: 200
- 紅花: 7
- 帖子: 789
- 在線: 170.7小時(shí)
- 蟲(chóng)號(hào): 2070126
- 注冊(cè): 2012-10-18
- 性別: GG
- 專業(yè): 生物信息學(xué)
★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
3樓2015-11-12 11:48:08
htt12119
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 9
- 帖子: 3
- 在線: 1.6小時(shí)
- 蟲(chóng)號(hào): 3992325
- 注冊(cè): 2015-07-27
2樓2015-10-17 09:29:57
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。 +4 | XX小鄧 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-25 | 7/350 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開(kāi)開(kāi)心心沒(méi)煩惱 2026-03-23 | 4/200 |
|