## ----样式,echo = false,结果='asis'---------------------------------------------------------------------------图书馆(BiocStyle)## --------------------------------------------------------------------------------------------------------------------------------------------库(Hilda)图书馆(TIDYR)库(GGPLOT2)库(GGPLOT2)库(dplyr)inputfile <-system.file(“ extdata/entdata/ayophageal.mp.txt.gz”,package =“ hilda”)g <-hildareadmpfile(inputfile,numbases = 5,trdir = 5,trdir = true)## -----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------库(SelectKsigs)load(System.File(“ Extdata/sample.rdata”,package =“ selectKsigs”))## --include = false --------------------------------------------------------------------------------------------------------- set.seed(5)结果<-cv_pmsignature(g,kfold = 3,nrep = 3,klimit = 7)print(结果)## -------------------------------------------------------------------------------------------------------------------------------------------------------结果$ kvalue <-seq_len(nrow(ress)) + 1结果_df < - 收集(结果,方法,方法,值,-kvalue)%>%group_by(方法)%>%突变(xmin = wher.min(value) + 1-0.1,xmax = wher.min(value) + 1 + 0.1)ggplot(resuction_df) + geom_point(aes(x = kvalue,y = value,color =方法),size = 2) + facet_wrap(〜方法,scales =“ free”) + geom_rect(mapping = aes = aes(xmin = xmin,xmax,xmax = xmax = xmax = xmax = xmax = xmax = xmax = xmax =xmax,ymin = -inf,ymax = inf),填充='灰色',alpha = 0.05) + them_bw() + xlab(“签名数”)## ------------------------------------------------------------------------------------------------------------------------------------------------------------- sessioninfo()