plot.R

x<-c(1:6)
y<-c(542, 182, 22, 7, 21, 66)
z<-c(330, 136, 42, 8, 39, 144)
a<-c(12, 3, 3, 13, 125, 267)
b<-c(15, 21, 3, 1, 77, 174)
c<-c(0,0,0,0,1,19)
plot(x, y, type="o",xaxt="n", col="red", pch=0, lwd=3, xlab="Blast e-value", ylab="Number of hits")
axis(side=1,at=c(1:6),labels=c("1E-51~1E-50", "1E-50~1E-40", "1E-40~1E-30","1E-30~1E-20", "1E-20~1E-10", "1E-10~1E-5"))
lines(x, z, type="o", col="green", pch=5, lwd=3)
lines(x, a, type="o", col="blue", pch=2, lwd=3)
lines(x, b, type="o", col="purple", pch=4, lwd=3)
lines(x, c, type="o", col="#FF9900FF", pch=16, lwd=3)
legend(locator(1), title="Type of Human papillomavirus", c("type 34", "type 18", "type 39", "type 56", "type 101"), pch=c(0, 5, 2, 4, 16), col=c("red", "green", "blue", "purple", "#FF9900FF"), lwd=c(3,3,3,3,3))