Bar and Line Plot

R Code used for this Plot:

year <- c(1980, 1981, 1982, 1983, 1984, 1985, 1986, 1987, 1988, 1989, 1990, 1991, 1992, 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013, 2014, 2015, 2016, 2017) 
measles.cases <- c(3852242, 4078455, 3623758, 3586101, 3026973, 2819553, 2078696, 1682847, 1663940, 1822390, 1325074, 1370243, 1481192, 1143945, 926808, 720391, 870218, 825673, 694466, 752407, 853479, 849173, 586471, 679856, 513406, 585701, 377576, 285031, 278751, 278637, 343806, 359332, 212376, 275307, 282078, 214816, 132328, 173330) 
vaccine.coverage <- c(17, 20, 21, 38, 42, 48, 47, 54, 63, 68, 73, 69, 69, 70, 71, 73, 73, 71, 71, 71, 72, 73, 72, 74, 76, 77, 79, 80, 81, 83, 84, 85, 84, 84, 84, 85, 85, 85) 

dat <- data.frame(year, measles.cases, vaccine.coverage) 
b <- barplot(height = dat$measles.cases/1000, names.arg = dat$year, xlab="Year", ylab="Mealses Cases", col=colors()[c(123)], main = "Annual Measles Cases and Vaccine Coverage between 1980-2017")

b
##       [,1]
##  [1,]  0.7
##  [2,]  1.9
##  [3,]  3.1
##  [4,]  4.3
##  [5,]  5.5
##  [6,]  6.7
##  [7,]  7.9
##  [8,]  9.1
##  [9,] 10.3
## [10,] 11.5
## [11,] 12.7
## [12,] 13.9
## [13,] 15.1
## [14,] 16.3
## [15,] 17.5
## [16,] 18.7
## [17,] 19.9
## [18,] 21.1
## [19,] 22.3
## [20,] 23.5
## [21,] 24.7
## [22,] 25.9
## [23,] 27.1
## [24,] 28.3
## [25,] 29.5
## [26,] 30.7
## [27,] 31.9
## [28,] 33.1
## [29,] 34.3
## [30,] 35.5
## [31,] 36.7
## [32,] 37.9
## [33,] 39.1
## [34,] 40.3
## [35,] 41.5
## [36,] 42.7
## [37,] 43.9
## [38,] 45.1
points(x=b[,1], y=dat$vaccine.coverage*40, type="o", col="red")
legend("topright", legend=c("Vaccine Coverage", "Measles Cases"), cex=0.5, fill=c("red", colors()[c(123)]))
axis(4, at=seq(0, 100, 10)*40, labels=seq(0, 100, 10), las=1, col="red", col.axis="red")
mtext("Population", side=4, line=3, col="red")