我要经过Deseq2的手册,有一些地方我仍然不明白:
DDS < - desqdatasetfrommatrix(countdata = countdata,coldata = coldata,design =〜cell + dex)
在设计中,细胞是细胞型,DEX含有不属实和治疗。在结果提取部分中,手册中的示例是
结果(DDS,对比度= C(“Cell”,“N061011”,“N61311)
我的理解:这是2个细胞类型之间的差异表达的结果:N061011和N61311,无论治疗还是不属于特征?
Puzzle1:如果我想在某种细胞类型中比较治疗和不属于的差异快递,请说明,N061011我该怎么办?
有了这个,我在结果中指导了结果功能这里但它提出了更多问题:
- 设计=〜基因型+条件+基因型的差异是什么:条件和设计=〜基因型+条件
例如服用2个基因型和2个符合条件;在手册中,它是:
DDS < - MakeExampledEsqDataset(n = 100,m = 12)DDS $ Genotype < - 因子(rep(c(c(c),“ii”),每个= 3),2))设计(DDS)< - 〜基因型+条件+基因型:条件DDS < - DESQ(DDS)结果Name(DDS)#基因型I(主要效果)结果的条件效果(DDS,对比度= C(“条件”,“B”,“A”))# the condition effect for genotype II # this is, by definition, the main effect *plus* the interaction term # (the extra condition effect in genotype II compared to genotype I). results(dds, list( c("condition_B_vs_A","genotypeII.conditionB") )) # the interaction term, answering: is the condition effect *different* across genotypes? results(dds, name="genotypeII.conditionB")
阅读后,我仍然有一些问题来通过:
2.在我的真实实验中,我可以知道哪些是主要影响?3.为什么对基因型I和基因型II的条件效果的命令是如此不同?
最后的命令结果(DDS,Name =“GenotypeII.ConditionB),它只提到了Genotypeii和ConditionB,在他的评论中,它被解释为基因型中的条件效果”差异“,我有点困惑。所以我的下一个问题是:
4.如何解释结果名称(DDS)的含义?
太多问题。谢谢你的时间!对这些问题的任何建议或我如何通过它更好地欣赏。
不胜感激!将练习您的建议。