varianceStabilizingTransformation

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

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

對(duì)于這個(gè)函數(shù)，一直在糾結(jié)varianceStabilizingTransformation（） 括號(hào)里面是加dds還是dds_2。運(yùn)行上述代碼發(fā)現(xiàn)vsd與vsd_2數(shù)據(jù)是完全一致的，所以不用過多糾結(jié)。。。。

image.png
image.png
最后編輯于 ：
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點(diǎn)，簡書系信息發(fā)布平臺(tái)，僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

  • References： 1.RNA-seq(4)：Hisat2+FeatureCounts+DESeq2流程+作圖...
    CuteCurse閱讀 2,068評(píng)論 0贊 2
  • 寫在前面的話 代碼中的# > 表示的是輸出結(jié)果 輸入 使用input()函數(shù) 用法 注意input函數(shù)輸出的均是字...
    FlyingLittlePG閱讀 3,213評(píng)論 0贊 9
  • 官網(wǎng) 中文版本 好的網(wǎng)站 Content-type: text/htmlBASH Section: User ...
    不排版閱讀 4,711評(píng)論 0贊 5
  • need:表達(dá)矩陣 分組矩陣，值要是整數(shù) DESeq2和EdgeR都可用于做基因差異表達(dá)分析，主要也是用于RNA-...
    白云夢(mèng)_7閱讀 46,539評(píng)論 1贊 24
  • library(pheatmap) data<- read.table("new 1.txt",header = ...
    Weiyx閱讀 718評(píng)論 0贊 0

友情鏈接更多精彩內(nèi)容