# #——包括= FALSE --------------------------------------------------------- knitr: opts_chunk美元集(崩溃= TRUE,发表评论 = "#>" ) ## 预先加载,以避免加载消息报告库(Glimma)库(limma)库(磨边机 ) ## ----------------------------------------------------------------------------- 库(Glimma)库(limma)图书馆(磨边机)dge < readRDS(执行(“RNAseq123 / dge。rds”,包= " Glimma ")) ## ----------------------------------------------------------------------------- glimmaMDS (dge ) ## ----------------------------------------------------------------------------- 设计< - readRDS(执行(“RNAseq123 /设计。rds"， package = "Glimma")) contre .matrix <- readRDS(system.file("RNAseq123/ contre .matrix. zip ")rds”,包= " Glimma ")) ## ----------------------------------------------------------------------------- v < -轰(dge、设计)vfit < - lmFit (v,设计)vfit < -对比。适合(vfit对比= contr.matrix) efit < - ebay (vfit ) ## ----------------------------------------------------------------------------- dge < - estimateDisp (dge设计)礼物< - glmFit (dge、设计)glrt < - glmLRT(礼物,设计,对比= contr.matrix ) ## ----------------------------------------------------------------------------- glimmaMA (efit dge = dge) # glimmaMA (glrt,dge = dge)使用磨边机的结果 ## ----------------------------------------------------------------------------- glimmaVolcano (efit dge = dge ) ## ---- eval = FALSE ------------------------------------------------------------ # htmlwidgets:: saveWidget (glimmaMA (efit dge = dge),“ma-plot.html”)# #你可以链接到使用[MA-plot] (ma-plot.html Rmarkdown ) ## ----------------------------------------------------------------------------- sessionInfo ()