绘制累积积分密度图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_