In kingsuching/bestpredictor: Predict It

knitr::opts_chunk$set(echo = TRUE)

library(githubinstall)
library(devtools)
library(tidyverse)
library(bestpredictor)
sf_salaries <- read.csv("/Users/sucheen/Documents/Cal Poly SLO/Academics/Spring 2021/STAT 331-01/Data/SF_Salaries_sub.csv")

sf_salaries <- sf_salaries[, unlist(lapply(sf_salaries, is.numeric))]
response = "TotalPay"
sf_salaries <- Filter(function(x)!all(is.na(x)), sf_salaries)
sf_salaries <- subset(sf_salaries, select = -c(factor(response)))
for(i in names(sf_salaries)) {
  f <- paste(response, "~", i)
  print(i)
  model <- lm(f, data = sf_salaries) %>% anova() %>%
    pluck("Sum Sq")
  print(model[1])
}

babies_ca <- read.csv("~/Documents/Cal Poly SLO/Academics/Spring 2021/STAT 331-01/Data/babies_ca.csv")

bestpredictor::predictors(babies_ca, "Count")

kingsuching/bestpredictor documentation built on Dec. 21, 2021, 6:42 a.m.