R语言绘制累计积分密度图（笔记）

绘制累积积分密度图code

######
######数据sum1_excel.xls
install.packages("readxl")
library(readxl)
sum_excel_data<-read_excel("路径.xls",col_names=TRUE)

temp_df<-as.data.frame(sum_excel_data)
temp_mt<-as.matrix(sum_excel_data)
sum_excel_data
xxx_data<-temp_df[,c(8,9,28,29)]
colnames(xxx_data)<-c("mapped reads number","mapped reads rate","unmapped reads number","unmapped reads rate")
ggplot(xxx_data, aes(x = xxx_data[,4])) + geom_density(colour = "cadetblue3")
group_mapped<-rep(c("Mapped"),185)
xxx_combine_mapped<-group_mapped
xxx_combine_mapped<-cbind(xxx_combine_mapped,temp_df[,c(8,9)])
group_unmapped<-rep(c("Unmapped"),185)
xxx_combine_unmapped<-group_unmapped
xxx_combine_unmapped<-cbind(group_unmapped,temp_df[,c(28,29)])
xxx_temp<-xxx_combine_mapped
all<-cbind(xxx_combine_mapped,xxx_combine_unmapped)
colnames(xxx_temp)<-c('group','number','rate')
colnames(xxx_combine_unmapped)<-c('group','number','rate')
ggplot(xxx,mapping = aes(x=Tatal_Reads/1e7,y=rate))+
  ######图形标题名称
  ggtitle("Mapped/Unmapped Reads")+
  ######累积积分密度图绘制
  stat_ecdf(aes(color=group),size=0.8)+
  ######x、y轴名称
  xlab("Total Reads[×1e7]")+ylab("Frequency")+
  ######x、y轴范围
  coord_cartesian(ylim = c(0, 1),xlim=c(0,3))+
  ######分组颜色标签名称
  scale_color_manual("group",                 # 修改图例名
                     values = c("red", "green"),     #manual: 直接指定分组使用的颜色
                     # breaks = c("4", "f", "r"),
                     labels = c(" Mapped", " Unmapped"))+
  
  #####刻度尺的步长
  scale_x_continuous(breaks=seq(0,4,1),minor_breaks=seq(0,5,by=1))+
  #scale_y_continuous(breaks=seq(0,1,0.15))+
  #scale_y_continuous(expand = c(0, 0))  使得坐标轴在原点上
  #删除网格
  theme_classic()+
  #####设置属性
  theme(  
      #坐标轴颜色、刻度线长度
      axis.text=element_text(colour="#36648B"),axis.ticks.length=unit(.4,"lines"),axis.ticks=element_line(#坐标轴刻度线的设置
    colour="#36648B",size=0.5,linetype=1,
    lineend=1),
      #坐标轴颜色的设置
    axis.line=element_line(colour="#36648B"),
      #x、y轴标题颜色、大小、字体、位置
    plot.title = element_text(color="#36648B", size=15, face="bold",vjust=0.6,hjust=0.6),
      #x、y轴标题颜色、大小、字体、位置
    axis.title.x = element_text(color="#36648B", size=12, face="plain"),
    axis.title.y = element_text(color="#36648B", size=12, face="plain")
    )+
  #####调用工具包调色板
  scale_color_brewer(type="qual", palette=6)