如何下载TCGA临床信息中的用药信息
这篇文章主要介绍了如何下载TCGA临床信息中的用药信息的相关知识,内容详细易懂,操作简单快捷,具有一定借鉴价值,相信大家阅读完这篇如何下载TCGA临床信息中的用药信息文章都会有所收获,下面我们一起来看看吧。
以下载用药信息为例:
#加载需要的包library(SummarizedExperiment)library(TCGAbiolinks)############################################################GDC:https://portal.gdc.cancer.gov/############################################################设置程序参数work_dir<-"/Users/zhangqiuxue/Downloads"#设置需要下载癌症对应的project和数据类型project<-"TCGA-GBM"data_category<-"Clinical"data_type<-"ClinicalSupplement"legacy<-FALSEfile_type="xml"#设置工作目录setwd(work_dir)#下载临床数据的结果DataDirectory<-paste0(work_dir,"/GDC/",gsub("-","_",project))#查询可以下载的数据query<-GDCquery(project=project,data.category=data_category,data.type=data_type,file.type=file_type,legacy=legacy)#该癌症总样品数量samplesDown<-getResults(query,cols=c("cases"))cat("TotalClinicalsampletodown:",length(samplesDown))#下载数据GDCdownload(query=query,directory=DataDirectory,files.per.chunk=6,method='client')#用专门的函数去整合下载好的数据clinical<-GDCprepare_clinic(query,clinical.info="drug",directory=DataDirectory)#将数据保存到文件,方便后面的进一步分析clinical_file<-paste0(DataDirectory,"_","clinical",".txt")write.csv(clinical,file=clinical_file,row.names=F,quote=F)
其中的关键就是设置:
file_type = "xml"
clinical.info = "drug"
关于“如何下载TCGA临床信息中的用药信息”这篇文章的内容就介绍到这里,感谢各位的阅读!相信大家对“如何下载TCGA临床信息中的用药信息”知识都有一定的了解,大家如果还想学习更多知识,欢迎关注恰卡编程网行业资讯频道。