首页游戏资讯生信实操 丨 不花钱就能验证数据,来看看这个数据库的“骚操作”

生信实操 丨 不花钱就能验证数据,来看看这个数据库的“骚操作”

misa2 03-06 3次浏览 0条评论

撰写: 米妮 来源:小张聊科研平台的“ i生信”公众号,微信公众号搜索“ i生信”即可关注/扫描关注见文末

#生信分析# #生信发文# #数据库#

文章看得多的小伙伴都知道TCGA数据库,大多研究者在文章里都利用过TCGA的数据来验证他们的数据分析结果,所以TCGA数据的下载和分析,作为一名生信人员也必会,虽然我也下过,但是每次都忘记自己之前怎么下,都得重新研究,所以这次记下来和有需要的小伙伴一起分享~

TCGA (The cancer genome atlas, 癌症基因组图谱),由美国National Cancer Institute(NCI)和National Human Genome Research Institute(NHGRI)于2006年联合启动的项目,目前TCGA已经收录了各种人类癌症包括每个癌症亚型在内的临床数据,基因组,转录组,甲基化等数据,是癌症研究者验证数据的重要数据来源。

官网会由于国内网问题,需要耐烦等加载,主页如图1所示:

图 1

1.主要按project来觅觅数据

2.主要按样本的实验信息和临床信息来找觅数据

3.网页分析功能

4.整合样本临床信息和样本数据信息来筛选数据

5.按样本的组织类型查询数据点击后与2的页面一致

()

下载manifest 文件

示例带大家一起下载下肾癌的表达矩阵抉择准则化后的FPKM表达矩阵,其它类似的数据下载类似,点击‘Repository’进进图2

展开全文

图 2

1.Files抉择样本得数据类型

2.抉择样本的临床特征

3.下载抉择好的manifast,利用此mnifast为后续数据下载做预备,此案例演

示下载kidney的表达矩阵,数据为准则化的FPKM值,样本包括

Kidney cancer的KIRC和KIRP

下载manifest选中样本的元数据

下载manifest 选中样本的原数据,一般包括样本信息和实验测序样本编号的信息,在刚才的抉择页面上点击如图3:

图 3

先把选好的样本添加到 Cart,然后下载metadata,如图4:

图 4

读进json在R中可以看到:元数据中有它的文件名,如图5:

图 5

数据下载工具——gdc-client的下载

数据的下载要利用gdc-client,下载地址:

下载过程看图6:

图 6

给大家演示的是在Windows上的,下载后解压,配置环境变量后就可以使用了,如图78:

图 7

图 8

现在指定刚才从TCGA下载的 manifest文件开始下载文件,下载在当前路径下,提前进进数据下载目录:

gdc-client download -m gdc_manifest.2021-03-16_kidney.txt

同样也是有时网不好会Connection error,需要多试几次,等数据下载好就可以了~

合并所有样本数据

下载好的文件是长这样的,如图910,也就是一个样本在一个文件夹,一个文件夹有一个txt文件,二列数据,整理起来要是不用代码,得发疯,下面利用代码R,就可以自动合并样本得到所有样本的表达矩阵。

图 9

图 10

利用代码实现自动读取每个文件夹下的下载数据:

setwd('D:/bioTree/800-emt-m6a/')

# count_files = dir("kidney_fpkm/",pattern = "*.htseq.counts.gz$",recursive = T)

# 1. kidney cancer 表达数据的整理

fpkm_files = dir("kidney_KIRC_KIRP/",pattern = "*.FPKM.txt.gz$",recursive = T)

ex = function(x){

result - read.table(file.path("kidney_KIRC_KIRP/",x),row.names = 1,sep = "\t")

return(result)

exp = lapply(fpkm_files,ex)

exp - do.call(cbind,exp)

dim(exp)

把样本数据对应到元数据中的样本名

上一步得到的exp是没有样本列名,需要读进meta信息,利用样本合并的文件顺序来到元数据中匹配样本名,把样本ID重新作为列名。

meta - jsonlite::fromJSON("metadata.cart.2021-03-18_KIRP_FPKM.json")

colnames(meta)

ids - meta$associated_entities;class(ids)

ids[[1]]

class(ids[[1]][,1])

ID = sapply(ids,function(x){x[,1]}) ## 文件id

file2id = data.frame(file_name = meta$file_name,

ID = ID)

head(file2id$file_name)

head(fpkm_files)

fpkm_files2 = stringr::str_split(fpkm_files,"/",simplify = T)[,2]

fpkm_files2[1] %in% file2id$file_name

file2id = file2id[match(fpkm_files2,file2id$file_name),]

colnames(exp) = file2id$ID

exp[1:4,1:4]

dim(exp)

得到的矩阵如下:

依据样本ID区分样本是癌症还是对比

Barcode含义介绍:

详尽介绍参见官网:

/

看样本类型主要看Sample类的编码:

接下来对表达矩阵进行样本组的区分就可以进行下游分析了,不花钱的验证数据就到手了~

原文链接:

趣书网txt官网
装修活在哪里找?装修师傅接单平台 P站突然公开!网友竟爱看这些片?!
相关内容
发表评论

游客 回复需填写必要信息