TCGA突变数据提取分组

##数据下载来源UCSC xena
mutype_file <- read.table('TCGA-BRCA.mutect2_snv.tsv',header = T,sep = '\t',quote = '')
dim(mutype_file)
mutype_file[1:5,1:5]
save(mutype_file,file = "BRCA_mutfile.Rdata")

##挑选突变数据
BRCA1 <- mutype_file[mutype_file$gene=='BRCA1',]
BRCA1[1:5,1:5]
dim(BRCA1)##27
##后可据此进行突变的分组，自行对sample进行相应匹配
##参考生信技能树提供的教程，由于很久之前写的可自行搜索到原文，本文仅供自己记录。
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容