# United States Department of Labor
# U.S. Bureau of Labor Statistics 
# Median usual weekly earnings (second quartile), Employed full time, Wage and salary workers

names=c("White Man","White Woman","Black Man","Black Woman","Latino Man","Latino Woman")

data = read.csv("earnings.csv",header=TRUE)

pdf(file="earnings.pdf",height=8,width=10)

par(mfrow=c(1,1))
plot(data[,1],data[,4],xlab="Year",ylab="Median earnigns per week",ylim=range(data[,4:9]),type="l",lwd=3)
lines(data[,1],data[,5],col=2,lwd=3)
lines(data[,1],data[,6],col=3,lwd=3)
lines(data[,1],data[,7],col=4,lwd=3)
lines(data[,1],data[,8],col=5,lwd=3)
lines(data[,1],data[,9],col=6,lwd=3)
legend("topleft",legend=names,col=1:6,lwd=3)

par(mfrow=c(1,1))
plot(data[,4],data[,5],xlim=range(data[,4:9]),ylim=range(data[,4:9]),xlab="Man",ylab="Woman",pch=16)
title("White")
abline(0,1)

par(mfrow=c(1,1))
plot(data[,6],data[,7],xlim=range(data[,4:9]),ylim=range(data[,4:9]),xlab="Man",ylab="Woman",pch=16,col=2)
title("Black")
abline(0,1)

par(mfrow=c(1,1))
plot(data[,8],data[,9],xlim=range(data[,4:9]),ylim=range(data[,4:9]),xlab="Man",ylab="Woman",pch=16,col=3)
title("Hispanic/Latino")
abline(0,1)

par(mfrow=c(1,1))
plot(data[,4],data[,5],xlim=range(data[,4:9]),ylim=range(data[,4:9]),xlab="Men",ylab="Women",pch=16)
abline(0,1,lty=3)
points(data[,6],data[,7],col=2,pch=16)
points(data[,8],data[,9],col=3,pch=16)
legend("topleft",legend=c("White","Black","Hispanic"),col=1:3,pch=16)

abline(lm(data[,5]~data[,4]))
abline(lm(data[,7]~data[,6]),col=2)
abline(lm(data[,9]~data[,8]),col=3)

lm(data[,5]~data[,4])$coef
lm(data[,7]~data[,6])$coef
lm(data[,9]~data[,8])$coef


par(mfrow=c(1,1))
plot(whiteman,blackman,xlim=range(data[,4:9]),ylim=range(data[,4:9]),xlab="",ylab="",pch=16)
abline(0,1)
points(whiteman,latinoman,col=2,pch=16)
points(blackman,latinoman,col=3,pch=16)
title("Man")
legend("topleft",legend=c("white x black","white x latino","black x latino"),col=1:3,pch=16)

par(mfrow=c(1,1))
plot(whitewoman,blackwoman,xlim=range(data[,4:9]),ylim=range(data[,4:9]),xlab="",ylab="",pch=16)
abline(0,1)
points(whitewoman,latinowoman,col=2,pch=16)
points(blackwoman,latinowoman,col=3,pch=16)
title("Woman")
legend("topleft",legend=c("white x black","white x latino","black x latino"),col=1:3,pch=16)


dev.off()