R/ridge-summary.R

In pbreheny/hdrm: High-Dimensional Regression Modeling

#' Reproduce Table 1.1
#'
#' Reproduces Table 1.1 from the book.  If you specify any options, your results may look different.
#'
#' @examples Tab1.1()
#'
#' @export

Tab1.1 <- function() {
  Data <- readData("pollution")
  XX <- std(Data$X)
  y <- Data$y
  fit.ridge <- ridge(XX, y)
  t.ridge <- summary(fit.ridge, which=which.min(fit.ridge$GCV))[,3]
  fit.ols <- lm(y ~ ., data=as.data.frame(XX))
  t.ols <- summary(fit.ols)$coef[,3]
  out <- cbind(t.ridge, t.ols)[-1,]
  ord <- order(out[,1], decreasing=TRUE)
  out[ord,]
}