varianceStabilizingTransformatio

2019-06-12 本文已影响0人装在瓶子里的阳光_2a70

dds <- DESeqDataSetFromMatrix(countData = exprSet,colData = colData, design = ~ group_list)
dds_2 <- DESeq(dds)

vsd <- assay(varianceStabilizingTransformation(dds))
vsd_2 <- assay(varianceStabilizingTransformation(dds_2))

对于这个函数，一直在纠结varianceStabilizingTransformation（）括号里面是加dds还是dds_2。运行上述代码发现vsd与vsd_2数据是完全一致的，所以不用过多纠结。。。。

image.png