Revision 616fe292ea4c66f8db35d04c9fa045fea56eea02 authored by Roger Koenker on 04 September 2016, 13:02:15 UTC, committed by cran-robot on 04 September 2016, 13:02:15 UTC 
version 5.29
1 parent db6345a
Raw File
engel1.R 
## Demo for a Plot of Engel curve in sample space

data(engel)

plot(foodexp ~ income, data = engel, cex= .5, col = "blue",
     xlab = "Household Income", ylab = "Food Expenditure")

z <- rq(foodexp ~ income, tau= .50, data = engel)# "median line": L1 - regression
abline(z, col = "dark blue")
abline(lm(foodexp ~ income, data = engel), lty=2, col="red") #the dreaded ols line

taus <- c(.05,.1,.25,.75,.90,.95)
nt <- length(taus)

for( i in 1:length(taus)) {
    abline(rq(foodexp~income, tau=taus[i], data = engel), col="gray")
}

legend("bottomright",
       c("L1 (tau = .50)", "OLS", paste("tau= ", formatC(rev(taus)))),
       col = c("dark blue", "red", rep("gray", nt)),
       lty = c(1,2, rep(1, nt)),
       inset = 0.03)
back to top