分组信息是在这个id的第14-15位，01-09是tumor，10-29是normal。我拿了一个示例数据，请在生信星球公众号回复0129获取。
根据这个生成一个分组信息它是一个向量形如
如果在id的14-15位在1-9之间就标记tumor，10-29之间就标记normal。这里面涉及到字符串截取、数据类型转换，%in%函数以及ifelse函数（if-else循环的变体）
将示例数据放到你的工作目录
load(file="id.Rdata")
table(substring(id,14,15)) #table看有多少重复值
num <- as.numeric(substring(id,14,15)) #截取出来是字符串，要转为数字
#屡试不爽的ifelse
group_list=ifelse(num %in% 1:9,"Tumor","Normal")
如果你要用自己的数据试试，就把样本编号命名为id，来跑一跑就ok。这一列信息后面作图会用。
非常开心公众号粉丝到了2000
0.背景 关于样本id和分组信息：TCGA的样本id里藏着分组信息[https://mp.weixin.qq.co...
分组信息是在这个id的第14-15位，01-09是tumor，10-29是normal。我拿了一个示例数据，请在生...
做差异分析需要的数据：表达矩阵和分组信息TCGA的数据只要表达矩阵就够了，因为其TCGA的样本ID比较特殊，样本I...
在分析TCGA数据库里的RNA-seq数据之前，先了解一下TCGA样本的id名称里的小秘密：参考文章：TCGA的样...
经常有人会问小编，TCGA下载的数据中，样本的ID号究竟代表什么意思。从样本ID号上能看出样本类型吗？首先小编先给...
本文适用于下载受保护的数据的文件，从文件id转换为TCGA样本id官网教程[https://docs.gdc.ca...
字段 分组id 分组名 分组描述名 分组描述信息 权限list 懒
TCGA数据下载方法：gdc-client，Xena和gdcRNAtools 需要下载的数据组学信息（样本）：（存...
在TCGA中，一个患者可能会对应多个样本，如TCGA-A6-6650可以得到3个样本数据： TCGA-A6-665...
我遇到的问题：从TCGA GDC把表达谱的数据下载到服务器后，发现表达矩阵有两列，一列是ENSG号，一列是表达值，...
本文标题：TCGA的样本id里藏着分组信息
本文链接：https://www.haomeiwen.com/subject/zosayqtx.html
网友评论