各位老师好,我这遇到deseq包问题,中间不报错而卡住。

attachments-2017-07-GwOX3Xtz597219422a8f这是卡住的截图

attachments-2017-07-I4aMWuek59721965872e这是代码,求解,谢谢

请先 登录 后评论

1 个回答

liuzhe

从你的代码和运行上面看,并不是对应的啊~~~

你试试betaPrior=FALSE,具体为什么原因:

(整理来自SEQanswers)

I have a dataset where the condition of interest is a factor with 3 levels. DESeq2 issues me a warning to use betaPrior=False otherwise the log2 FoldChange will be affected by the base level. I thought the fold-change is always calculated w.r.t to a base level.

I tried to look in the DESeq2 paper at bioRXiv (http://biorxiv.org/content/early/2014/02/19/002832). In the paper, we say "logarithmic fold change prior", and as logarithmic fold changes are written as beta in the GLM, in the software we use the shorter "betaPrior". 

请先 登录 后评论