## ---- para，echo = false，结果='hide'--------------------------------------------------------------------------- biocStyle :: markdown（）knitr :: opts_chunk $ set（dev =“ png”，fig.show =“ hold”，fig.width = 4，fic.height = 4.5，fig.align =“中心“，消息= false，collapse = true）## ----库------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------库（RNASEQCOMP）## ----数据-------------------------------------------------------------------------------------------------------------＃将数据集加载到此软件包数据（SIMDATA）类（SIMDATA）名称（SIMDATA）## ---- META --------------------------------------------------------------------------------------------------------------------------------------------------- condinfo < - 因子（simdata $ samp $条件）repinfo < - 因子（simdata $ samp $ replicate）estaruationFeature <- rep（true，nrow（simdata $ meta））calibrationFeature <-simdata $meta $ house＆simdata $ meta $ chr =='chr1'umitReference <-1 ## ----过滤器--------------------------------------------------------------------------------------------------------------------------------------------------------- dat <-simailcalibrate（simdata $ quant），condinfo，repinfo，评估表演，校准图，统一，校准ature2 = calibrationFeature）类（dat）显示（dat）## -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------1.4））## --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- plotne（dat，xlim = c（0.5,1））## -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- plot2TX(dat,genes=simdata$meta$gene,ylim=c(0,0.6)) ## ----diffroc---------------------------------------------------------------------------------------------------- plotRoc（dat，simdata $ meta $阳性，simdata $ meta $ fcsign，ylim = c（0,0.8））## ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- simdata$meta$fcsign[simdata$meta$fcstatus == "off.on"] <- NA plotFC(dat,simdata$meta$positive,simdata$meta$fcsign,ylim=c(0,1.2))