网站首页 > 基础教程 正文
本文较长,将遇到问题过程写了进来,请耐心观看。
此类问题百度了很久,暂无详细教程,因此做了。
【一】
以下报错你肯定很纠结吧:
【01】rownames(design)=colnames(mRNAdata)
Error in dimnames(x) <- dn : 'dimnames'的长度[1]必需与陈列范围相等
【02】row.names(mRNAdata) <- mRNAdata[, 1]
Error in .rowNamesDF<-(x, value = value) : 不允许有重复的'row.names'
此外: Warning message:
non-unique values when setting 'row.names': ‘MATR3’, ‘PINX1’, ‘SIGLEC5’, ‘TMSB15B’
【03】rownames(mRNAdata)<- mRNAdata[,1]
Error in .rowNamesDF<-(x, value = value) : 'row.names'里不允许有遗漏值
【二】
一般的生信分析步骤是这样的:
数据库下载→数据整理【ID转换,GTF文件筛选出需要的基因等】→差异分析DEG→差异分析可视化【PCA,箱式图,火山图,热图......】→各种操作。。。。。。
巧了,今天又碰到钉子了。
进行DEG分析的时候,发现常规代码不能运行。
group_list=ifelse(as.numeric(substr(colnames(mRNAdata),14,15)) < 10,'tumor','normal')
design <- model.matrix(~0+factor(group_list))
colnames(design)=levels(factor(group_list))
rownames(design)=colnames(mRNAdata)
design
rownames(design)=colnames(mRNAdata):【此代码运行出错】
以为是运行的代码有问题,对比其他代码,大同小异。
那就是数据处理有问题。
仔细对比自己数据和教程数据:
教程图片如下:
自己图片如下:
看出了区别吗?区别就在红色方框部分。
找到问题那就简单了:
row.names(mRNAdata) <- mRNAdata[, 1]
mRNAdata <- mRNAdata[, -1]
再次报错:
这里是需要去除重复值吗?:
#去除重复着代码:
index <- duplicated(d$gene_name)
table(index)
d <- d[!index,]
rownames(d) <- d$gene_name
d$gene_name <- NULL
#或者这样的代码:
mRNAdata=distinct(d,gene_name,.keep_all = T)
可是还是报错。
回到点,我们是要干嘛?把gene_name变成行名。
以前使用代码时候没有出错。
难道是DEG分析包错误?我一般习惯用limma包。没办法,换成了熟悉的DESeq包。再次分析。
rm(list = ls())
load("mRNA_exprSet_dds_sample.Rdata")
res <- results(dds, tidy=TRUE) #获取结果
res <- as_tibble(res)
require(dplyr)
res <- res %>%
separate(row,c("symbol","ensemble","genetype"),sep = " \\| ") %>%
dplyr::select(- c(ensemble,genetype)) %>%
arrange(desc(abs(log2FoldChange))) %>% #排序。为了去重
distinct(symbol,.keep_all = TRUE) %>%
arrange(desc(log2FoldChange))#再次按照log2FoldChange从大到小排序
save(res,file = "result.Rda")
res
然而,还是出问题了。每次都是在【把gene_name变成行名。】这个地方出问题。
看来跟包没关系。
图片展示:
查看此处错误,再次百度。
row.names(res) <- res[, 1]
Error in `.rowNamesDF<-`(x, value = value) : 'row.names'的长度不对
此外: Warning message:
Setting row names on a tibble is deprecated.
【三】
tibble is deprecated. 这是个啥?百度。
离真相越来越近了,自己却浑然不知。
百度看了几个答案后,大概明白了,是tidyverse这个包出问题了。查看
https://github.com/tidyverse/tibble/issues/288这个网址,原来是开发包的人在搞事情。他把这个包的部分功能给GG了。
事已至此,找到真凶。可是既然不能用,怎么才能实现我们想要的功能呢?
再次百度。某乎有个回答。大概内容就是,先将数据保存为“.CVS”格式,然后Excel打开,然后借助Excel的功能,直接将“symbol”设置成行名。
https://zhuanlan.zhihu.com/p/404403216。
数据稍小可以,大数据的时候崩了怎么办?
还得R出手吧!
再次百度寻求答案。
https://github.com/tidyverse/tibble/issues/288里面有答案。
最终解决方法如下:【需手动设置】
rownames(res)=as.vector(as.matrix(res[,1]))
自己把res$symbol这一行去掉就行了。
需要注意的是,res$symbol这一行由数据框格式变为:
class(res$symbol)
[1] "character"
【03】rownames(mRNAdata)<- mRNAdata[,1]
Error in .rowNamesDF<-(x, value = value) : 'row.names'里不允许有遗漏值
这个解决办法:
#这里报错,有遗漏值。因此进行一下操作。
sum(is.na(mRNAdata$gene_name))
which(is.na(mRNAdata$gene_name))
mRNAdata<-mRNAdata[-203,]
sum(is.na(mRNAdata$gene_name))
大功告成。
【四】
总结一下:
由于包的缘故,导致rownames调用失败。所以这个地方报错。其实仔细观察报错截图,每次都有rownames他的身影。
gene_name有时候会出现NA。需要去除后再分析。
温馨提示:表达矩阵的第一列是样品还是基因名,如果第一列是基因名,要把第一列设置成为行名,不然只有的差异分析会出错。
因此,学会了吗?解决不容易,请点赞,收藏,评论等,文明三连,拒绝白嫖。
猜你喜欢
- 2024-10-12 R数据分析:倾向性评分匹配实例操作
- 2024-10-12 要为学习神经网络奠定基础,你需要认真读读R深度学习
- 2024-10-12 怎样快速入门Arduino?(二十三)—TCS3200颜色传感器
- 2024-10-12 把数据输入R之后,如何进行简单的操作(一)
- 2024-10-12 聚类分析5—物种集合-数量生态学:R语言的应用 第四章
- 2024-10-12 一文读懂R中的探索性数据分析(附R代码)
- 2024-10-12 从0开始自制解释器——添加对括号的支持
- 2024-10-12 常见编程语言:Go:Go语言函数与方法
- 2024-10-12 Python基础之变量、循环、函数(二)
- 2024-10-12 R语言中使用scan函数从键盘获取数据的方法
- 最近发表
- 标签列表
-
- gitpush (61)
- pythonif (68)
- location.href (57)
- tail-f (57)
- pythonifelse (59)
- deletesql (62)
- c++模板 (62)
- css3动画 (57)
- c#event (59)
- linuxgzip (68)
- 字符串连接 (73)
- nginx配置文件详解 (61)
- html标签 (69)
- c++初始化列表 (64)
- exec命令 (59)
- canvasfilltext (58)
- mysqlinnodbmyisam区别 (63)
- arraylistadd (66)
- node教程 (59)
- console.table (62)
- c++time_t (58)
- phpcookie (58)
- mysqldatesub函数 (63)
- window10java环境变量设置 (66)
- c++虚函数和纯虚函数的区别 (66)