# #——包括= FALSE --------------------------------------------------------- knitr: opts_chunk美元集(崩溃= TRUE,发表评论 = "#>" ) ## 预先加载,以避免加载消息报告库(Glimma)库(磨边机)库(DESeq2 ) ## ----------------------------------------------------------------------------- 库(Glimma)库(磨边机)图书馆(DESeq2) dge < readRDS(执行(“RNAseq123 / dge。rds”,包= " Glimma ")) dds < DESeqDataSetFromMatrix (countData = dge数美元,美元colData = dge样本,rowData = dge $基因,设计= ~ ) ## ----------------------------------------------------------------------------- glimmaMDS (dds ) ## ----------------------------------------------------------------------------- dds < - DESeq (dds,安静= TRUE ) ## ----------------------------------------------------------------------------- glimmaMA (dds ) ## ---- eval = FALSE ------------------------------------------------------------ # # 在工作目录中创建ma-plot.html # #在Rmarkdown链接使用[MA-plot] (ma-plot.html) # htmlwidgets:: saveWidget (glimmaMA (dds)、“ma-plot.html ") ## ----------------------------------------------------------------------------- sessionInfo ()