如何下載TCGA臨床信息中的用藥信息

這篇文章主要介紹了如何下載TCGA臨床信息中的用藥信息的相關(guān)知識(shí)，內(nèi)容詳細(xì)易懂，操作簡(jiǎn)單快捷，具有一定借鑒價(jià)值，相信大家閱讀完這篇如何下載TCGA臨床信息中的用藥信息文章都會(huì)有所收獲，下面我們一起來(lái)看看吧。

以下載用藥信息為例：

#  加載需要的包
library(SummarizedExperiment)
library(TCGAbiolinks)

###########################################################
# GDC: https://portal.gdc.cancer.gov/
###########################################################

# 設(shè)置程序參數(shù)
work_dir <- "/Users/zhangqiuxue/Downloads" 

# 設(shè)置需要下載癌癥對(duì)應(yīng)的project 和數(shù)據(jù)類(lèi)型
project <- "TCGA-GBM"
data_category <- "Clinical"
data_type <- "Clinical Supplement"
legacy <- FALSE
file_type = "xml"

# 設(shè)置工作目錄
setwd(work_dir)

# 下載臨床數(shù)據(jù)的結(jié)果
DataDirectory <- paste0(work_dir,"/GDC/",gsub("-","_",project))

# 查詢可以下載的數(shù)據(jù)
query <- GDCquery(project = project,
                  data.category = data_category,
                  data.type = data_type, 
                  file.type = file_type,
                  legacy = legacy)

# 該癌癥總樣品數(shù)量
samplesDown <- getResults(query,cols=c("cases"))
cat("Total Clinical sample to down:", length(samplesDown))

# 下載數(shù)據(jù)
GDCdownload(query = query,
            directory = DataDirectory,files.per.chunk=6, method='client')

# 用專(zhuān)門(mén)的函數(shù)去整合下載好的數(shù)據(jù)
clinical <- GDCprepare_clinic(query, clinical.info = "drug",directory = DataDirectory)

# 將數(shù)據(jù)保存到文件，方便后面的進(jìn)一步分析
clinical_file <- paste0(DataDirectory, "_","clinical",".txt")
write.csv(clinical, file = clinical_file, row.names = F, quote = F)

其中的關(guān)鍵就是設(shè)置：

file_type = "xml"

clinical.info = "drug"

關(guān)于“如何下載TCGA臨床信息中的用藥信息”這篇文章的內(nèi)容就介紹到這里，感謝各位的閱讀！相信大家對(duì)“如何下載TCGA臨床信息中的用藥信息”知識(shí)都有一定的了解，大家如果還想學(xué)習(xí)更多知識(shí)，歡迎關(guān)注億速云行業(yè)資訊頻道。