dreamfly1993 2021-01-12 19:31
浏览 565

R语言gsub函数处理TCGA数据库样品名的问题

TCGA样品名:TCGA-H6-8124-11A-01R-2404-07;TCGA-H6-8124-01A-01R-2404-07

按“-”分隔,第四组11A,1开头为正常样品,01A,0开头为肿瘤样品。选择0开头样品。

#按“-”分隔,选第四组

group=sapply(strsplit(rownames(countst),"\\-"),"[",4)

#选第一个数字

group=sapply(strsplit(group,""),"[",1)

#这步是什么意思?gsub是替换函数,我运行了一下,不加这步后面group==0就运行不下去,提示:Error in data[, group == 0] : (下标)逻辑下标太长。但加了这步得到的还是1和0的结果,对数据类别和值都没有影响,这是为啥。
group=gsub("2","1",group)

#选择等于0的项
data=data[,group==0]

group=sapply(strsplit(rownames(countst),"\\."),"[",4)
group=sapply(strsplit(group,""),"[",1)
group=gsub("2","1",group)
data=data[,group==0]

谢谢~

  • 写回答

0条回答 默认 最新

    报告相同问题?

    悬赏问题

    • ¥15 R语言Rstudio突然无法启动
    • ¥15 关于#matlab#的问题:提取2个图像的变量作为另外一个图像像元的移动量,计算新的位置创建新的图像并提取第二个图像的变量到新的图像
    • ¥15 改算法,照着压缩包里边,参考其他代码封装的格式 写到main函数里
    • ¥15 用windows做服务的同志有吗
    • ¥60 求一个简单的网页(标签-安全|关键词-上传)
    • ¥35 lstm时间序列共享单车预测,loss值优化,参数优化算法
    • ¥15 Python中的request,如何使用ssr节点,通过代理requests网页。本人在泰国,需要用大陆ip才能玩网页游戏,合法合规。
    • ¥100 为什么这个恒流源电路不能恒流?
    • ¥15 有偿求跨组件数据流路径图
    • ¥15 写一个方法checkPerson,入参实体类Person,出参布尔值