百度360必应搜狗淘宝本站头条
当前位置:网站首页 > 热门文章 > 正文

R语言rownames函数报错:生信差异分析DEG时报错

bigegpt 2024-12-25 10:23 7 浏览


本文较长,将遇到问题过程写了进来,请耐心观看。
此类问题百度了很久,暂无详细教程,因此做了。

【一】

以下报错你肯定很纠结吧:
【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。需要去除后再分析。

温馨提示:表达矩阵的第一列是样品还是基因名,如果第一列是基因名,要把第一列设置成为行名,不然只有的差异分析会出错。

因此,学会了吗?解决不容易,请点赞,收藏,评论等,文明三连,拒绝白嫖。

相关推荐

得物可观测平台架构升级:基于GreptimeDB的全新监控体系实践

一、摘要在前端可观测分析场景中,需要实时观测并处理多地、多环境的运行情况,以保障Web应用和移动端的可用性与性能。传统方案往往依赖代理Agent→消息队列→流计算引擎→OLAP存储...

warm-flow新春版:网关直连和流程图重构

本期主要解决了网关直连和流程图重构,可以自此之后可支持各种复杂的网关混合、多网关直连使用。-新增Ruoyi-Vue-Plus优秀开源集成案例更新日志[feat]导入、导出和保存等新增json格式支持...

扣子空间体验报告

在数字化时代,智能工具的应用正不断拓展到我们工作和生活的各个角落。从任务规划到项目执行,再到任务管理,作者深入探讨了这款工具在不同场景下的表现和潜力。通过具体的应用实例,文章展示了扣子空间如何帮助用户...

spider-flow:开源的可视化方式定义爬虫方案

spider-flow简介spider-flow是一个爬虫平台,以可视化推拽方式定义爬取流程,无需代码即可实现一个爬虫服务。spider-flow特性支持css选择器、正则提取支持JSON/XML格式...

solon-flow 你好世界!

solon-flow是一个基础级的流处理引擎(可用于业务规则、决策处理、计算编排、流程审批等......)。提供有“开放式”驱动定制支持,像jdbc有mysql或pgsql等驱动,可...

新一代开源爬虫平台:SpiderFlow

SpiderFlow:新一代爬虫平台,以图形化方式定义爬虫流程,不写代码即可完成爬虫。-精选真开源,释放新价值。概览Spider-Flow是一个开源的、面向所有用户的Web端爬虫构建平台,它使用Ja...

通过 SQL 训练机器学习模型的引擎

关注薪资待遇的同学应该知道,机器学习相关的岗位工资普遍偏高啊。同时随着各种通用机器学习框架的出现,机器学习的门槛也在逐渐降低,训练一个简单的机器学习模型变得不那么难。但是不得不承认对于一些数据相关的工...

鼠须管输入法rime for Mac

鼠须管输入法forMac是一款十分新颖的跨平台输入法软件,全名是中州韵输入法引擎,鼠须管输入法mac版不仅仅是一个输入法,而是一个输入法算法框架。Rime的基础架构十分精良,一套算法支持了拼音、...

Go语言 1.20 版本正式发布:新版详细介绍

Go1.20简介最新的Go版本1.20在Go1.19发布六个月后发布。它的大部分更改都在工具链、运行时和库的实现中。一如既往,该版本保持了Go1的兼容性承诺。我们期望几乎所...

iOS 10平台SpriteKit新特性之Tile Maps(上)

简介苹果公司在WWDC2016大会上向人们展示了一大批新的好东西。其中之一就是SpriteKitTileEditor。这款工具易于上手,而且看起来速度特别快。在本教程中,你将了解关于TileE...

程序员简历例句—范例Java、Python、C++模板

个人简介通用简介:有良好的代码风格,通过添加注释提高代码可读性,注重代码质量,研读过XXX,XXX等多个开源项目源码从而学习增强代码的健壮性与扩展性。具备良好的代码编程习惯及文档编写能力,参与多个高...

Telerik UI for iOS Q3 2015正式发布

近日,TelerikUIforiOS正式发布了Q32015。新版本新增对XCode7、Swift2.0和iOS9的支持,同时还新增了对数轴、不连续的日期时间轴等;改进TKDataPoin...

ios使用ijkplayer+nginx进行视频直播

上两节,我们讲到使用nginx和ngixn的rtmp模块搭建直播的服务器,接着我们讲解了在Android使用ijkplayer来作为我们的视频直播播放器,整个过程中,需要注意的就是ijlplayer编...

IOS技术分享|iOS快速生成开发文档(一)

前言对于开发人员而言,文档的作用不言而喻。文档不仅可以提高软件开发效率,还能便于以后的软件开发、使用和维护。本文主要讲述Objective-C快速生成开发文档工具appledoc。简介apple...

macOS下配置VS Code C++开发环境

本文介绍在苹果macOS操作系统下,配置VisualStudioCode的C/C++开发环境的过程,本环境使用Clang/LLVM编译器和调试器。一、前置条件本文默认前置条件是,您的开发设备已...