p1 <- ggplot(evidence,aes(x= evidence$'Mass.Error..ppm.'))+
geom_histogram(binwidth = 150,fill="steelblue",color="white",size=0.1)+
theme_classic2(base_size = 18)+xlim(-5,5)+
labs(title = "Peptide mass tolerance distribution",x="Mass.Error..ppm.")
missed.cleavages <- data.frame(table(evidence$Missed.cleavages))
colnames(missed.cleavages)=c("Missed\ncleavages","Count")
missed.cleavages$Perentage=sprintf("%.1f%%",missed.cleavages$Count/sum(missed.cleavages$Count)*100)
evidence$'Missed.cleavages'=as.factor(evidence$'Missed.cleavages')
evidence$Length=as.factor(evidence$Length)
p1
Warning messages:
1: Use of `evidence$Mass.Error..ppm.` is discouraged. Use `Mass.Error..ppm.` instead.
2: Removed 117997 rows containing non-finite values (stat_bin).


雷达卡






京公网安备 11010802022788号







