两条命令就跑完RNA差异分析

    数据分析实在太难了,差异分析到底怎么跑啊,怎么才不出错呀!coldata到底怎么构建呀……!别慌,傻瓜包来了!它就是simpleRNA,居然只要两条命令就能跑完复杂的RNA差异分析,并生成结果矩阵和结果图!

simpleRNA —— 简单易用的R包,不需要懂分析就能跑

    为了方便各位不熟悉数据分析的小白们能够自己做好RNA差异分析,小编联合另一位小伙伴Skye直接写了一个傻瓜包——simpleRNA。10秒就能跑完

    你需要准备的只有2个:

1> 你的表达矩阵

2> 一个分组矩阵

    运行流程也非常简单,只要你电脑装了R,在R中指定这两个矩阵的路径,然后直接复制2行代码就可以直接跑出结果啦!

test <- runDiff(matdir = matdir,condi = condi)comp <- compareSmp(test$dds,group1 = "CV",group2 = "normal")

结果会直接输出标准化过的矩阵两两比较的结果、以及表达火山图和热图

数据来源于【GSE81455】

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE81455

是不是很心动!

接下来,小编就教大家怎么用。

使用教程

【一、安装】(高阶玩家可以直接跳过)

1> 首先你的电脑要安装R,然后安装依赖包(已经装了这些的话请跳过):

      依赖包:DESeq2, edgeR, pheatmap, ggplot2, ggrepel【强者跳过】

      没装过的直接复制以下代码

install.packages("ggplot2")install.packages("ggrepel")install.packages("pheatmap")if (!requireNamespace("BiocManager", quietly = TRUE))  install.packages("BiocManager")BiocManager::install("DESeq2")BiocManager::install("edgeR")

2> 安装simpleRNA,下载simpleRNA安装包【获取安装包请拉到最底下

      在Rstudio中右下角点开packages选项卡,点击install

install from 选择第二个package Archieve file,在弹出的对话框中找到刚刚下载的package文件

点击install即可--安装完成。

【二、准备输入文件】

1> 表达矩阵

仅支持txt文件,如果你的表达矩阵是excel文件,请打开该文件,然后执行 【CTRL+A】,【CTRL+C】, 然后新建一个txt文件,并【CTRL+V】,然后直接保存。 (其实就是全选然后复制粘贴过去)

【example】:

第一列第一行必须是空

【全选】——【复制】——【粘贴】后如下:


2> 分组矩阵

分组矩阵的建立步骤:

    a> 在EXCEL中建立分组矩阵

【注意分组矩阵的列名和表达矩阵的行名一致】

【第一列第一行必须是空,第二列第一行必须是condition

    b> 复制粘贴到txt文件

【三、运行】

#导入所有需要的包library(DESeq2)library(edgeR)library(ggplot2)library(ggrepel)library(pheatmap)library(simpleRNA)matdir <- "C:/R/wd/expression_Matrix.txt"  [路径]请换成你自己的路径condi <- "C:/R/wd/condition.txt"  [路径]同样需要更换成你自己的路径test <- runDiff(matdir = matdir,condi = condi)comp <- compareSmp(test$dds,group1 = "CV",group2 = "normal")

    记得把路径换成你自己的路径哦。

    这样就完成了,接下来会在你的工作路径里面生成上调结果和下调结果,以及一个plot文件包含一个【火山图】和一个【热图】,如果你不知道工作目录在哪里,没关系,运行完会有提示的。

提示我已经生成了结果,来看看吧。

自动生成的火山图:

以及热图:

数据来源于【GSE81455】经过了修改,将第一行第一列改成了空

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE81455

还不错!省了很多事!

各位小伙伴如果有什么别的希望能做成傻瓜包的,欢迎给我留言哦~

直接戳公众号回复即可!

【simpleRNA】下载:回复公众号关键词simpleRNA,即可获得下载地址啦!转发让更多的小伙伴用上~

特别鸣谢:

simpleRNA共同作者Skye

隔壁桌师妹不厌其烦的调试

转载自公众号文章两条命令就跑完RNA差异分析

本文使用 文章同步助手 同步

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 193,968评论 5 459
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 81,682评论 2 371
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 141,254评论 0 319
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 52,074评论 1 263
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 60,964评论 4 355
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 46,055评论 1 272
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 36,484评论 3 381
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 35,170评论 0 253
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 39,433评论 1 290
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 34,512评论 2 308
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 36,296评论 1 325
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 32,184评论 3 312
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 37,545评论 3 298
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 28,880评论 0 17
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,150评论 1 250
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 41,437评论 2 341
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 40,630评论 2 335

推荐阅读更多精彩内容