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)