Deseq2

mycounts<-read.csv("gene_count_matrix.csv",row.names=1)
head(mycounts)
dim(mycounts)
mycounts_1<-mycounts[rowSums(mycounts) !=0,]
dim(mycounts_1)
mymeta<-read.csv("mymeta.csv",stringsAsFactors =T)
mymeta
colnames(mycounts_1)==mymeta$id
library(DESeq2)
dds<-DESeqDataSetFromMatrix(countData=mycounts_1,
                            colData=mymeta,
                            design=~dex)
dds<-DESeq(dds)
res<-results(dds)
head(res)
class(res)
res_1<-data.frame(res)
library(dplyr)
res_1 %>%
  mutate(group = case_when(
    log2FoldChange >= 2 & padj <= 0.05 ~"up",
    log2FoldChange <= -2 & padj <=0.05 ~"DOWN",
    TRUE ~ "NOT_CHANGE"
  )) -> res_2
table(res_2$group)
write.csv(res_2,file="R_S_diff_expr_result.csv",
          quote = F)

需要文件：

gene_count_matrix.csv
mymeta.csv
file=“R_S_diff_expr_result.csv”