# NOT RUN {
api_key <- "abcdefghijklmnopqrstuvwxyz123456"
# Naive Phillips Curve Estimation
df <- fred(p = "CPIAUCSL", u = "u6rate", key=api_key)
head(df)
df$pi <- log(df$p) - log(c(NA, df$p[2:nrow(df)]))
reg_pc <- lm(pi ~ u, data = df)
summary(reg_pc)
# Different variations of the CPI
cpis <- c("CPIMEDSL", "CPIFABSL", "CPIHOSSL")
names(cpis) <- c("P_Med", "P_Food", "P_House")
df_cpi <- fred(cpis, key=api_key, long=TRUE)
head(df)
# }
Run the code above in your browser using DataLab