DESeq2 결과

R로 하는 Bioconductor 기반 RNA-Seq

Mary Piper

Bioinformatics Consultant and Trainer

DESeq2 결과 테이블

mcols(wt_res)

DESeq2 결과 열 설명

R로 하는 Bioconductor 기반 RNA-Seq
head(wt_res, n=10)

log2 fold change (MAP): condition fibrosis vs normal 
Wald test p-value: condition fibrosis vs normal 
data frame with 6 rows and 6 columns
                           baseMean   log2FoldChange             lfcSE             stat
                          <numeric>        <numeric>         <numeric>        <numeric>
ENSMUSG00000102693                0               NA                NA               NA
ENSMUSG00000064842                0               NA                NA               NA
ENSMUSG00000051951 19.5084656230804 3.55089043143673 0.648400500074659 4.66871842838828
ENSMUSG00000102851                0               NA                NA               NA
ENSMUSG00000103377                0               NA                NA               NA
ENSMUSG00000104017                0               NA                NA               NA
                                 pvalue                 padj
                              <numeric>            <numeric>
ENSMUSG00000102693                   NA                   NA
ENSMUSG00000064842                   NA                   NA
ENSMUSG00000051951 3.03084428526558e-06 1.93776447202312e-05
ENSMUSG00000102851                   NA                   NA
ENSMUSG00000103377                   NA                   NA
ENSMUSG00000104017                   NA                   NA
R로 하는 Bioconductor 기반 RNA-Seq

유의미한 DE 유전자 - 요약

summary(wt_res)

DESeq2 요약

R로 하는 Bioconductor 기반 RNA-Seq

유의미한 DE 유전자 - 배수 변화 임계값

wt_res <- results(dds_wt, 
                  contrast = c("condition", "fibrosis", "normal"), 
                  alpha = 0.05, 
                  lfcThreshold = 0.32)

wt_res <- lfcShrink(dds_wt, 
                    contrast=c("condition", "fibrosis", "normal"), 
                    res=wt_res)
R로 하는 Bioconductor 기반 RNA-Seq

유의미한 DE 유전자 - 요약

summary(wt_res)

DESeq2 요약 lfc 임계값

R로 하는 Bioconductor 기반 RNA-Seq

결과 - 주석 추가

library(annotables)

grcm38

annotables_grcm38 테이블

R로 하는 Bioconductor 기반 RNA-Seq

결과 - 추출

wt_res_all <- data.frame(wt_res) %>% 
              rownames_to_column(var = "ensgene") %>%
              left_join(x = wt_res_all, 
                        y = grcm38[, c("ensgene", "symbol", "description")], 
                        by = "ensgene")                    
View(wt_res_all)

병합된 전체 결과

R로 하는 Bioconductor 기반 RNA-Seq

유의미한 DE 유전자 - 정렬

wt_res_sig <- subset(wt_res_all, padj < 0.05)

wt_res_sig  <- wt_res_sig %>%  
        arrange(padj)

View(wt_res_all)

병합된 전체 결과

R로 하는 Bioconductor 기반 RNA-Seq

유의미한 유전자

R로 하는 Bioconductor 기반 RNA-Seq

연습해 봅시다!

R로 하는 Bioconductor 기반 RNA-Seq

Preparing Video For Download...