def drawer_curve(R_result): ''' input 1: set_minus_set_minus.xlsx output 1: set minus curve plot input and output are all in R_result ''' R_code_draw_curve = ''' drawer_curve=function(result_R){ require(tidyverse) require(readxl) in_fl=read_xlsx( paste(result_R,"set_minus_set_minus.xlsx",sep = "") ) in_fl$species=factor(in_fl$species,levels = in_fl$species) plot_line=ggplot( data = in_fl, mapping=aes( x=as.numeric(row.names(in_fl)), # y=minus group=1 ) )+ geom_line( aes( y=minus ), size=1 )+ geom_line( aes( y=add ), size=1 )+ # scale_y_continuous( # breaks = c(8511,9000,10000,11000,12000,12110), # labels = c(8511,9000,10000,11000,12000,12110) # )+ theme( panel.background = element_blank(), axis.line = element_line(size=1, colour = "black"), axis.title = element_text(size=20), axis.text = element_text(size = 20) )+ xlab(label = "Number of strain")+ ylab(label = "Number of gene")+ geom_hline( yintercept = c(max(in_fl$add),min(in_fl$minus)),lty=3,lwd=1,alpha=0.8 )+ annotate( "text", x = 110, y = 9000, label="Core", size=10 )+ annotate( "text", x = 110, y = 14500, label="Pan", size=10 ) ggsave( paste(result_R,"curve.png",sep = ""), plot = plot_line, width = 8, height = 5 ) } ''' R_darw_curve = SignatureTranslatedAnonymousPackage(R_code_draw_curve, "R_darw_curve") R_darw_curve.drawer_curve(R_result)