##数据下载来源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进行相应匹配
##参考生信技能树提供的教程,由于很久之前写的可自行搜索到原文,本文仅供自己记录。
TCGA突变数据提取分组
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 从果子老师的从TCGA数据中提取lncRNA并进行下游分析一文中,我学到了如何从TCGA表达谱数据中提取RNA。老...
- 页面解析和数据提取 一般来讲对我们而言,需要抓取的是某个网站或者某个应用的内容,提取有用的价值。内容一般分为两部分...
- 本人长期出售超大量微博数据、旅游网站评论数据,并提供各种指定数据爬取服务,Message to YuboonaZh...
- 今天在完成pm的时候,发现了这样的一条需求。我的层次表中有三条数分别是专升本,高起本,专科三条数据。但是查询分组的...