零代码学生信第八天-什么是差异分析?怎么做差异分析?

提到差异分析,大家都不陌生,什么是差异分析呢?差异分析的目的是为了什么呢?怎么做差异分析呢?

首先差异分析是什么?

差异分析又可以被称为差距分析,主要就是为了查看两组之间影响两组之间差距较大的因素有哪些。

        对于基因芯片的差异分析而言,由于芯片数据基本都是服从正态分布。由于一次性找的基因多,因此需要对多重试验进行校正,控制假阳性,目前基因芯片数据分析用的最多的就是limma包。limma包使用是广义线性模型的一种,主要就是通过对每一个基因的表达量拟合成一个线性方程。limma的分析包括ANOVA分析、线性回归等。该方法的核心步骤包括voom、fit、eBays等步骤。

        而对于高通量测序的差异分析分析而言,由于高通量的read count 基本都是服从泊松分布,不能直接使用limma包筛选差异。一般情况下,使用DEseq2包和edgeR包对read count的高通量测序数据进行差异分析,这两个包都是基于负二项分布的模型,将数据进行转化,让其满足正态分布。

总之,不管是使用什么方法进行差异分析,核心点都是分为如下几步:

1、预处理,不管是什么样的表达数据,数据的预处理是非常重要的,过滤/取对数,过滤的目的是为了去除某一些基因在所有样本的表达变化较小的基因,或者去除表达量较低的基因;而取对数的目的其实是为了伸缩数据,让数据更符合正态分布。

2、差异计算,即计算两组比较中,不同基因在两组的表达的差距,如用A组的基因比B组的基因的表达量,得到差异倍数,差异倍数越大,说明该基因在A组和B组的表达差别越大,同样的也说明该基因也就越重要。

3、p值的计算,或者是FDR计算。在计算差异后,同样的也需要计算该基因在两组中的表达是否存在偶然性,因此会根据数据分布的类型选择合适的统计方法。如数据符合正态分布,样本量小于30使用t检验,样本量大于30使用u检验,如果数据符合偏态分布,使用秩和检验。为了控制假阳性的存在,使用多重检验对p值进行校正,方法有很多,如BH等等。

目前做差异分析的工具很多,很多都是需要用到编程的,这对于不会代码的人来说,很不友好。有没有一种工具,可以做到直接输入表达谱矩阵,直接作出差异分析的结果呢?

这里推荐生信豆芽菜平台


如果是芯片数据,直接用limma包做差异,如果是高通量测序数据,count数据则使用edgeR包或者是deseq2包做差异分析,如果是tpm/fpkm则先取log,在进行差异分析。

上传准备的文件都是一样的,需要准备一个基因的表达谱数据,和一个分组文件,记得比较组和被比较组尽量不要出现中文,或者空格这些


表达谱数据


行名不能重复,如果出现duplicate 'row.names' are not allowed,就要用去重工具去重

分组数据:有两列,第一列为样本,第二列为分组(最好只写两组),不要出现空行,如果不确定可以选择下载notepad进行查看。


运行后等到成功即可。

如果不清楚可以查看左侧的说明文档,写的很详细。

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

推荐阅读更多精彩内容