Trabalho 1

Use o mesmo banco de dados acima sobre renda. Rode um modelo estratificado por n?vel educacional:

No High School:NOHS

High School:HS

bachelor’s degree: BS

PHD

summary(income\(NOHS) summary(income\)HS) summary(income\(BS) summary(income\)PHD)

regress?o estratificada

mod.nohs <- lm(INCTOT ~ YEAR, data=income, subset=(NOHS==0)) mod.Withs <- lm(INCTOT ~ YEAR, data=income, subset=(NOHS==1)) mod.wouths <- lm(INCTOT ~ YEAR, data=income, subset=(HS==0)) mod.hs <- lm(INCTOT ~ YEAR, data=income, subset=(HS==1)) mod.woutbs <- lm(INCTOT ~ YEAR, data=income, subset=(BS==0)) mod.bs <- lm(INCTOT ~ YEAR, data=income, subset=(BS==1)) mod.woutphd <- lm(INCTOT ~ YEAR, data=income, subset=(PHD==0)) mod.phd <- lm(INCTOT ~ YEAR, data=income, subset=(PHD==1))

summary(mod.nohs) summary(mod.Withs) summary(mod.wouths) summary(mod.hs) summary(mod.woutbs) summary(mod.bs) summary(mod.woutphd) summary(mod.phd)

par(mfrow=c(4,2), mar=c(3,3,2,1))

1. NOHS = 1

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.Withs$coefficients, lwd=1.2) title(main=“No High School (YES)”)

2. NOHS = 0

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.nohs$coefficients, lwd=1.2) title(main=“No High School (NO)”)

3. HS = 1

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.hs$coefficients, lwd=1.2) title(main=“High School (YES)”)

4. HS = 0

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.wouths$coefficients, lwd=1.2) title(main=“High School (NO)”)

5. BS = 1

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.bs$coefficients, lwd=1.2) title(main=“Bachelor’s (YES)”)

6. BS = 0

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.woutbs$coefficients, lwd=1.2) title(main=“Bachelor’s (NO)”)

7. PHD = 1

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.phd$coefficients, lwd=1.2) title(main=“PhD (YES)”)

8. PHD = 0

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.woutphd$coefficients, lwd=1.2) title(main=“PhD (NO)”)

#SEGUNDA Opção para plot #high school par(mfrow=c(1,2))

1. NOHS = 1

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.Withs$coefficients, lwd=1.2) title(main=“No High School”)

2. HS = 1

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.hs$coefficients, lwd=1.2) title(main=“High School”)

bachelor’s degree:

par(mfrow=c(1,2))
# 3. BS = 1 plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.bs$coefficients, lwd=1.2) title(main=“Bachelor’s (YES)”)

4. BS = 0

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.woutbs$coefficients, lwd=1.2) title(main=“Bachelor’s (NO)”)

PHD

5. PHD = 1

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.phd$coefficients, lwd=1.2) title(main=“PhD (YES)”)

6. PHD = 0

plot(income\(YEAR, income\)INCTOT, xlab=“YEAR”, ylab=“Income”, xlim=c(1990,2010), ylim=c(0,400000), pch=’’) abline(mod.woutphd$coefficients, lwd=1.2) title(main=“PhD (NO)”)