sleuth基于TPM值sleuth進行的差異分析是怎樣的

這期內(nèi)容當中小編將會給大家?guī)碛嘘Psleuth基于TPM值sleuth進行的差異分析是怎樣的，文章內(nèi)容豐富且以專業(yè)的角度為大家分析和敘述，閱讀完這篇文章希望大家可以有所收獲。

創(chuàng)新互聯(lián)建站服務項目包括相城網(wǎng)站建設、相城網(wǎng)站制作、相城網(wǎng)頁制作以及相城網(wǎng)絡營銷策劃等。多年來，我們專注于互聯(lián)網(wǎng)行業(yè)，利用自身積累的技術優(yōu)勢、行業(yè)經(jīng)驗、深度合作伙伴關系等，向廣大中小型企業(yè)、政府機構等提供互聯(lián)網(wǎng)行業(yè)的解決方案，相城網(wǎng)站推廣取得了明顯的社會效益與經(jīng)濟效益。目前，我們服務的客戶以成都為中心已經(jīng)輻射到相城省份的部分城市，未來相信會繼續(xù)擴大服務區(qū)域并繼續(xù)獲得客戶的支持與信任！

kallisto等alignment-free轉錄本定量軟件，會給出TPM值的定量結果?；谶@種類型的結果進行差異分析時，有兩種策略可以選擇。

第一種是采用tximportR包，將結果導入到DESeq2種進行分析；第二種是直接采用sleuthR包進行差異分析。本章主要介紹sleuth的使用。

這個包的源代碼存放在github上，鏈接如下

https://github.com/pachterlab/sleuth

github上的R包其安裝方式比較特殊, 具體過程如下

source("http://bioconductor.org/biocLite.R")
biocLite("rhdf5")
library(devtools)
install_github("pachterlab/sleuth")

首先從Bioconductor上安裝依賴的rhdf5包，因為kallisto的定量結果為HDF5格式，這個R包用來讀取數(shù)據(jù)，然后采用devtools這個R包，自動從github的源代碼進行安裝。

所有差異分析需要的都是定量結果和樣本分組這兩個基本元素，只不過不同的R包要求的格式不同。在sleuth中，將這兩種信息存儲在一個三列的數(shù)據(jù)框中，示例如下

> s2c
    samples   group              paths
1 control-1 control kallisto/control-1
2 control-2 control kallisto/control-2
3 control-3 control kallisto/control-3
4    case-1    case    kallisto/case-1
5    case-2    case    kallisto/case-2
6    case-3    case    kallisto/case-3

第一列為樣本名稱，第二列為樣本對應的分組信息，第三列為每個樣本kallisto定量結果的文件夾。通過這樣的一個數(shù)據(jù)框，就包含了差異分析所需的所有信息。

假定有6個樣本，分成control,case 兩組， 每組3個生物學重復，可以通過以下代碼構建上述的數(shù)據(jù)框

samples = c(
"control-1",
"control-2",
"control-3",
"case-1",
"case-2",
"case-3")

s2c <- data.frame(
samples = samples,
group   = rep(c("control", "case"), each = 3),
paths   = paste("kallisto", samples, sep = "/")
)

上述代碼要求將所有樣本的定量結果放在同一個文件夾下，目錄結構如下

kallisto/
├── control-1
├── control-2
├── control-3
├── case-1
├── case-2
└── case-3

上述數(shù)據(jù)框準備好之后，就可以讀取數(shù)據(jù)進行差異分析了，完整的代碼如下

library(sleuth)
so <- sleuth_prep(s2c, extra_bootstrap_summary = TRUE)
so <- sleuth_fit(so, ~condition, 'full')
so <- sleuth_fit(so, ~1, 'reduced')
so <- sleuth_lrt(so, 'reduced', 'full')
sleuth_table <- sleuth_results(so, 'reduced:full', 'lrt', show_all = FALSE)

上述就是小編為大家分享的sleuth基于TPM值sleuth進行的差異分析是怎樣的了，如果剛好有類似的疑惑，不妨參照上述分析進行理解。如果想知道更多相關知識，歡迎關注創(chuàng)新互聯(lián)行業(yè)資訊頻道。

文章標題：sleuth基于TPM值sleuth進行的差異分析是怎樣的
標題路徑：http://muchs.cn/article32/ihdspc.html

成都網(wǎng)站建設公司_創(chuàng)新互聯(lián)，為您提供定制網(wǎng)站、網(wǎng)站策劃、云服務器、網(wǎng)站制作、自適應網(wǎng)站、外貿(mào)建站

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉載內(nèi)容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉載，或轉載時需注明來源： 創(chuàng)新互聯(lián)