source('../../Data/get_vital_rates_data.R')
DF = get_vital_rates_data()

par(mfrow=c(2, 1))
m1 = lm(All_cause ~ Year, data=DF)
plot(DF$Year, DF$All_cause, type='p', pch=19, xlab='Year', ylab='Mortality (all causes)')
abline(m1, col='red')
grid()

m2 = lm(Respiratory_disease ~ Year, data=DF)
plot(DF$Year, DF$Respiratory_disease, type='p', pch=19, xlab='Year', ylab='Respirator disease')
abline(m2, col='red')
grid()
par(mfrow=c(1, 1))