这篇文章主要介绍了怎么用limma包进行芯片数据差异分析的相关知识,内容详细易懂,操作简单快捷,具有一定借鉴价值,相信大家阅读完这篇怎么用limma包进行芯片数据差异分析文章都会有所收获,下面我们一起来看看吧。

利用limma 包进行芯片数据差异分析,需要按照分组信息构建分组矩阵,基于model.matrix函数完成,获得的分组矩阵结构类似下方显示结果:

> design
           A B
GSM1626001 1 0
GSM1626002 1 0
GSM1626003 1 0
GSM1626010 0 1
GSM1626011 0 1
GSM1626012 0 1
attr(,"assign")
[1] 1 1
attr(,"contrasts")
attr(,"contrasts")$`description`
[1] "contr.treatment"

其0 1表示有无,可以理解为A组存在GSM1626001 GSM1626002 GSM1626003 三个样品,B组理解依然,整个矩阵描述了将要比较的两组样品情况,之后利用makeContrasts构建差异对比矩阵 以便继续进行后续分析

      Contrasts
Levels B - A
     A    -1
     B     1

表示开始进行B相对于A组的差异表达分析