24小时热门版块排行榜    

查看: 1826  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

快乐君子

银虫 (小有名气)


[交流] 利用R对RNA-seq的reads 数据进行统计分析

我现在有一组RNA-seq的 reads count数据,这个实验有11种处理(1个control),每一处理3重复,所以用表格整合总的有33列reads count 数据,
我现在想用R中的DEseq2 package算出log2 fold change,我能否对平行实验的reads count取平均?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhkzhou

木虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
你的问题没有表述清楚,1个control 是指没有做重复?
我猜是(1个control + 10 个处理)*3。
一直强调 RNA-seq 要有重复,你若是取平均,则完全丧失了数据的优势。
DEseq2和edgeR的说明书里都非常详细的列举了计算过程。
16楼2015-07-24 20:31:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

whd2013

木虫 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
滥发金币的效果就是一堆表情加祝福了

[ 发自小木虫客户端 ]
14楼2015-07-20 19:52:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见