DEGseq 流程

Degseq 用法
读入文件需要是表达矩阵

library(DESeq2)
library(tidyverse)

# https://www.sohu.com/a/400667791_120736615
# https://www.jianshu.com/p/d473470e6d0c

exp <- read.table('./all.count.txt',header = TRUE, row.names = 1)
exp <- exp[,-c(1:5)]

col_n <- colnames(exp)
col_n <- grep()

design <- factor(c(rep("control",22), rep('test',19)))
dds <- DESeqDataSetFromMatrix(exp, DataFrame(design), design= ~ design )
dds <- DESeq(dds)
results <- results(dds)

# Results with mean,change,p-value
results <- as.data.frame(results)

after_deg <- cbind(exp,results)