DESeq分析基因的差异表达以及安装中出现的问题

65 阅读 0 评论 43 点赞

我是靠谱客的博主妩媚项链，最近开发中收集的这篇文章主要介绍DESeq分析基因的差异表达以及安装中出现的问题，觉得挺不错的，现在分享给大家，希望可以做个参考。

概述

DESeq采用NB（负二项分布检验的方式）对reads数进行差异显著性检验，同时还增加了矫正由于长度引起的误差，估算基因表达量的方式采用basemean值来估算表达量（标准化以后）这里，我在R下安装DESeq包出现了一些问题帮助总结一下。

首先，我的R是15.0版本

打开R
source("http://bioconductor.org/biocLite.R")
biocLite（‘DEseq’）
会出现安装成功
loading的时候缺少locfit包则需要再输入（locfit包只在15.1版本下编译的，不过貌似没什么问题）
biocLite（‘locfit’）
再loading会出现缺少xtable包
输入
biocLite（‘xtable’）
在loading会出现XML报错
输入
biocLite（‘XML’）
这样DESeq就运行没有问题了。

注中间可能会出现某个包无法安装的情况，注意检查网络连接

下面具体的如何进行差异显著性检验，我参考别人的数据，这里计算出reads在每个样本中各个基因上的counts数文件，下面couns_...分别表示6个样本（三个重复对照三个重复处理组）

library(DESeq)

setwd("wherethe data is")