發(fā)簡信
IP屬地：江蘇
-
關(guān)于heatmap一、熱圖的介紹 （1）什么是熱圖 熱圖是heatmap的直譯，用暖色表示數(shù)值大，冷色表示數(shù)值小。行為基因，列為樣本。 （2）熱圖的作用 數(shù)據(jù)質(zhì)量控制：通過熱圖，我們能很直觀的...
-
-
努力學(xué)習(xí)
-
道可拾遺 ·文章寫得不錯，但是有個問題
res = results(dds2, contrast=c("condition", "control", "treat")) 這個地方應(yīng)該是存在問題的
下面是help文件中result中contrast參數(shù)的描述
this argument specifies what comparison to extract from the object to build a results table. one of either:
a character vector with exactly three elements: the name of a factor in the design formula, the name of the numerator level for the fold change, and the name of the denominator level for the fold change (simplest case)
a list of 2 character vectors: the names of the fold changes for the numerator, and the names of the fold changes for the denominator. these names should be elements of resultsNames(object). if the list is length 1, a second element is added which is the empty character vector, character(). (more general case, can be to combine interaction terms and main effects)
a numeric contrast vector with one element for each element in resultsNames(object) (most general case)
If specified, the name argument is ignored.
可見，當(dāng)contrast參數(shù)有三個元素時依次代表 數(shù)據(jù)中作為treat/control標(biāo)簽的列的名稱, fold change的分子， 以及fold change的分母。
fold change 一般算的時treat/control，作者也在他的描述中說他算的是treat versus control，所以這里正常應(yīng)該是
res = results(dds, contrast=c("condition", "treat", "control"))
在作者后面分析部分
> diff_gene_deseq2 <-subset(res, padj < 0.05 & abs(log2FoldChange) > 1)
> head(diff_gene_deseq2)
log2 fold change (MLE): condition control vs treat
也展示了他這里算的log2 fold change是control vs treat -
-
"#"后面就是樓主的注釋
RNAseq分析（5）：差異分析工具 GDCRNATools 安裝及測試
前言 GDCRNATools 是一個用于下載、整理和綜合分析GDC中IncRNA、mRNA和miRNA數(shù)據(jù)的R/Bioconductor包。主要功能包括：差異基因分析、生存分...
-
喜歡寫代碼的外科醫(yī)生。
-
