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数据是完全一致的,所以不用过多纠结。。。。