参照的是bioconductor workflow里面的RNAseq123教程。
输入class(x)的结果如下：
[1] "DGEList"
attr(,"package")
[1] "edgeR"
而x本身包括三个列表x$samples x$counts 和x$genes
[b]上一步操作为：keep.exprs <- rowSums(cpm>1)>=3 返回结果为TRUE或者FALSE
下一步操作为：[/b]
[b]x <- x[keep.exprs,, keep.lib.sizes=FALSE] #筛选符合上面条件的基因，可是这种筛选完全不知道该如何解读。[/b]
[b]求指导！[/b]
x的形成代码如下：
[attach]180[/attach]
x的数据格式如下，好像是是个多维数组，最后一个$genes没有展示出来
[attach]181[/attach]
阅读全文
收起全文
