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进行相应匹配
##参考生信技能树提供的教程，由于很久之前写的可自行搜索到原文，本文仅供自己记录。