打开APP
userphoto
未登录

开通VIP,畅享免费电子书等14项超值服

开通VIP
万事开头难

基因集富集分析(GSEA)中的排序指标:它们重要吗?

费九牛二虎之力也无法重现的GSEA图

以下是中文

    前几天看了生信技能树老师的推文,关于曹雪涛老师课题组在Nature Medcine杂志上发表的题为“Tumor-educated B cells selectively promote breast cancer lymph node metastasis by HSPA4- targeting lgG”的最新研究成果,由于之前看过这篇新鲜出炉的文章,不久前还听过曹老师的学术报告,所以对这篇推文的关注度很高,再加上推文里有老师重复文献时整理的一整套R代码,如果至宝,于是就操刀实干起来。

万万没想到,第一行代码就报错了,果真是万事开头难!!!

具体的报错如下:

提示我第19行没有3个元素,不服气的我在GEO官网上手动下载文件,打开发现文件正常,不存在第19行没有3个元素的问题。于是我上网搜索了很久,发现大家都是在read.table 或者类似的读取函数中出现了类似的报错,而且提出的解决办法对我没有帮助,我只好向撰写推文的老师发了求助邮件,没想到老师很快就回复了我,老师建议我找一个网络速度快的地方运行代码,或者翻墙下载。我照做了,可是还是报错,我又和老师沟通,老师给我发来了他运行程序的最原始R代码文件,可是结果还是不可以。在老师的再一次建议下,我换了一台电脑运行,令人心碎的是,还是报错。再换别的电脑结果肯定还是如此,于是我及时止损。最后我尝试在实验室的服务器上运行代码,结果跑通了。于是我在想,在普通R上是仅仅下载这个GEO数据跑不通,还是都不行?于是我又尝试了别的GEO文件,结果一样的语句,只是文件编号不同,程序跑通了。可是对于这个GSE113249为什么在服务器上能跑通,在普通的R上就不行呢?难道是这个文件本身某处藏着玄机,使得它更适合在服务器上运行?

万事开头难,但是办法总会有的,看看程序通顺了之后制作出来的热图,如十里春风,暖人心房。

后记

小作者其实并没有彻底解决问题,亲爱的读者们,你能帮帮她吗?

本站仅提供存储服务,所有内容均由用户发布,如发现有害或侵权内容,请点击举报
打开APP,阅读全文并永久保存 查看更多类似文章
猜你喜欢
类似文章
【热】打开小程序,算一算2024你的财运
2021第一期生信入门微信群答疑精选200题
Python IDLE(python GUI)与python(comand line)有什么区别
PermissionError: [WinError 5] 拒绝访问。: 'C:\\Users\\Administrator\\Desktop\\001'
Tip:outlook无法通过CAS服务器更新脱机通讯薄
《Android模块化项目实战》08 服务器部署与报错处理写字字帖
关于Maven+SVN/Git+Dubbo项目开发实用的小技巧
更多类似文章 >>
生活服务
热点新闻
分享 收藏 导长图 关注 下载文章
绑定账号成功
后续可登录账号畅享VIP特权!
如果VIP功能使用有故障,
可点击这里联系客服!

联系客服