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

m1 = lm(Percent_in_2005 ~ Percent_in_1985, data=DF)
print(summary(m1))

if( T ){
    par(mfrow=c(2, 1))
    plot(DF$Percent_in_1985, DF$Percent_in_2005, type='p', pch=19, col='black', xlab='Percent in 1985', ylab='Percent in 2005')
    abline(m1, col='blue')
    grid()

    plot(fitted(m1), residuals(m1), type='p', pch=19, col='black', xlab='fitted values', ylab='residuals')
    abline(h=0, col='black')
    grid()
    par(mfrow=c(1, 1))
}


DF$Percent_Change = DF$Percent_in_2005 / DF$Percent_in_1985

m2 = lm(Percent_Change ~ Percent_in_1985, data=DF)
print(summary(m2))

if( T ){
    par(mfrow=c(2, 1))
    plot(DF$Percent_in_1985, DF$Percent_Change, type='p', pch=19, col='black', xlab='Percent in 1985', ylab='Percent Change 1985-2005')
    abline(m2, col='blue')
    grid()

    plot(fitted(m2), residuals(m2), type='p', pch=19, col='black', xlab='fitted values', ylab='residuals')
    abline(h=0, col='black')
    grid()
    par(mfrow=c(1, 1))
}