# devtools::install_github("pbreheny/hdrm")
library(hdrm)

# Slide 12
Fig1.1()

# Slide 19
out <- Ex1.1()
Fig1.2(out)

# Slide 21
n <- 25
b <- out$results[,"Estimate"]
mean(b^2)
12/n
mean(b^2)/(12/n)

# Slide 23
mean(out$pred^2)

# Slide 24
median(out$results[,"Pr(>|t|)"])
range(out$results[,"Pr(>|t|)"])
mean(out$cover)