### R code from vignette source '/home/zeileis/svn/projects/AER/Course/Ch-Intro.Rnw' ################################################### ### code chunk number 1: setup ################################################### options(prompt = "R> ", continue = "+ ", width = 64, digits = 4, show.signif.stars = FALSE, useFancyQuotes = FALSE) library("AER") set.seed(1071) ################################################### ### code chunk number 2: journals-data ################################################### library("AER") data("Journals", package = "AER") dim(Journals) names(Journals) ################################################### ### code chunk number 3: journals-plot (eval = FALSE) ################################################### ## plot(log(subs) ~ log(price/citations), data = Journals) ################################################### ### code chunk number 4: journals-lm (eval = FALSE) ################################################### ## j_lm <- lm(log(subs) ~ log(price/citations), data = Journals) ## abline(j_lm) ################################################### ### code chunk number 5: journals-lmplot1 ################################################### plot(log(subs) ~ log(price/citations), data = Journals) ################################################### ### code chunk number 6: journals-lmplot2 ################################################### plot(log(subs) ~ log(price/citations), data = Journals) j_lm <- lm(log(subs) ~ log(price/citations), data = Journals) abline(j_lm) ################################################### ### code chunk number 7: journals-lmsummary ################################################### summary(j_lm) ################################################### ### code chunk number 8: CPS1985-data ################################################### data("CPS1985", package = "AER") cps <- CPS1985 ################################################### ### code chunk number 9: CPS1985-lm ################################################### cps_lm <- lm(log(wage) ~ experience + I(experience^2) + education, data = cps) ################################################### ### code chunk number 10: CPS1985-fitted ################################################### cps2 <- data.frame(education = mean(cps$education), experience = min(cps$experience):max(cps$experience)) cps2 <- cbind(cps2, predict(cps_lm, newdata = cps2, interval = "prediction")) ################################################### ### code chunk number 11: CPS1985-plot (eval = FALSE) ################################################### ## plot(log(wage) ~ experience, data = cps) ## lines(fit ~ experience, data = cps2, col = 2) ################################################### ### code chunk number 12: CPS1985-plot1 ################################################### plot(log(wage) ~ experience, data = cps) lines(fit ~ experience, data = cps2, col = 2) ################################################### ### code chunk number 13: CPS1985-quantreg ################################################### library("quantreg") cps_rq <- rq(log(wage) ~ experience + I(experience^2) + education, data = cps, tau = seq(0.2, 0.8, by = 0.15)) ################################################### ### code chunk number 14: CPS1985-fitted2 ################################################### cps2 <- cbind(cps2, predict(cps_rq, newdata = cps2)) ################################################### ### code chunk number 15: CPS1985-rq-plot (eval = FALSE) ################################################### ## plot(log(wage) ~ experience, data = cps) ## for(i in 6:10) lines(cps2[,i] ~ experience, ## data = cps2, col = 2) ################################################### ### code chunk number 16: CPS1985-rq-summary (eval = FALSE) ################################################### ## plot(summary(cps_rq)) ################################################### ### code chunk number 17: CPS1985-rq-plot1 ################################################### plot(log(wage) ~ experience, data = cps) for(i in 6:10) lines(cps2[,i] ~ experience, data = cps2, col = 2) ################################################### ### code chunk number 18: CPS1985-rq-summary1 ################################################### plot(summary(cps_rq)) ################################################### ### code chunk number 19: CPS1985-bkde ################################################### library("KernSmooth") cps_bkde <- bkde2D(cbind(cps$experience, log(cps$wage)), bandwidth = c(3.5, 0.5), gridsize = c(200, 200)) ################################################### ### code chunk number 20: CPS1985-bkde-plot (eval = FALSE) ################################################### ## image(cps_bkde$x1, cps_bkde$x2, cps_bkde$fhat, ## col = rev(gray.colors(10, gamma = 1)), ## xlab = "experience", ylab = "log(wage)") ## box() ## lines(fit ~ experience, data = cps2) ## lines(lwr ~ experience, data = cps2, lty = 2) ## lines(upr ~ experience, data = cps2, lty = 2) ################################################### ### code chunk number 21: CPS1985-bkde-plot1 ################################################### png(file = "Ch-Intro-CPS1985-bkde.png", height = 500, width = 650) par(mar = c(5, 5, 2, 4)) image(cps_bkde$x1, cps_bkde$x2, cps_bkde$fhat, col = rev(gray.colors(10, gamma = 1)), xlab = "experience", ylab = "log(wage)") box() lines(fit ~ experience, data = cps2) lines(lwr ~ experience, data = cps2, lty = 2) lines(upr ~ experience, data = cps2, lty = 2) dev.off() ################################################### ### code chunk number 22: objects ################################################### objects() ################################################### ### code chunk number 23: search ################################################### search() ################################################### ### code chunk number 24: objects2 (eval = FALSE) ################################################### ## objects("package:base") ################################################### ### code chunk number 25: create_objects ################################################### x <- 2 x objects() ################################################### ### code chunk number 26: remove_objects ################################################### remove(x) objects() ################################################### ### code chunk number 27: log1 ################################################### log(x = 16, base = 2) ################################################### ### code chunk number 28: log2 (eval = FALSE) ################################################### ## log(16, 2) ## log(x = 16, 2) ## log(16, base = 2) ## log(base = 2, x = 16) ################################################### ### code chunk number 29: quit (eval = FALSE) ################################################### ## q() ################################################### ### code chunk number 30: options (eval = FALSE) ################################################### ## ?options ## help("options") ################################################### ### code chunk number 31: example (eval = FALSE) ################################################### ## example("options") ## example("lm") ################################################### ### code chunk number 32: apropos ################################################### apropos("help") ################################################### ### code chunk number 33: vignette (eval = FALSE) ################################################### ## vignette("strucchange-intro", package = "strucchange") ################################################### ### code chunk number 34: demo (eval = FALSE) ################################################### ## demo() ################################################### ### code chunk number 35: demo2 (eval = FALSE) ################################################### ## demo("graphics")