Assignment 3: Normalization and OPS+
For assignment three, we worked on understanding normalization and OPS+. For OPS+, it was important to calculate each year's average on base percentage and slugging average for the league. Then using the OPS+ we are able to create an overall function, making this usable to any time frame.
Also, this assignment required us to use the same data and compare it to the winning percentage of the teams. The data was displayed in a scatter plot and a linear model with R and R-Squared was created. While the plot gives a visual comparison, the linear model allows us to see if there is a correlation between the winning percentage and the OPS+. It also allows us to determine if the model is predictive.
Download
library(tidyverse)
library(ggpubr)
library(gt)
library(readr)
#Importing Team file
teamsA3 <- read_csv("MATH 494/Teams (for Assignment 3)(in).csv")
#Creating subset for years 2000-2019
teamSub <- filter(teamsA3, yearID %in% 2000:2019)
teamSub <- select(teamSub, yearID, teamID, franchID, G, W, L, H, BB, HBP,
AB, SF, `2B`, `3B`, HR, BPF)
#Get OBP and SLG for each year and by team
teamFun <- function(df){
OBP <- with(df, round((H+BB+HBP)/(AB+BB+HBP+SF),3))
TB <- with(df, (H+`2B`+2*`3B`+3*HR))
SLG <- with(df, round(TB/AB,3))
df <- cbind(df, OBP, SLG)
}
teamSub <- teamFun(teamSub)
#Get Averages for OBP and SLG for each team
teamSub |>
group_by(yearID) |>
summarise(across(c(OBP, SLG), mean),
.groups="keep") -> sumOPS
names(sumOPS)[2] <- "avgOBP"
names(sumOPS)[3] <- "avgSLG"
teamSub |>
left_join(sumOPS |>
select(yearID, avgOBP, avgSLG),
join_by(yearID)) -> teamSub
#Get OPS+ for each year and by team
teamOPSF <- function(df){
tOBP <- with(df, OBP/avgOBP)
tSLG <- with(df, SLG/avgSLG)
dOPS <- with(df, (tOBP+tSLG-1)/(BPF))
OPSP <- with(df, round(dOPS*100, 3))
df <- cbind(df, OPSP)
}
teamSub <- teamOPSF(teamSub)
#Create table for the team with the best OPS+ by year
tableSub <- select(teamSub, yearID, teamID, OPSP)
tableSub |>
group_by(yearID) |>
top_n(1, OPSP) |>
ungroup(yearID)-> tableSub
tableSub <- gt(select(tableSub, yearID, teamID, OPSP),
rowname_col = "yearID") |>
tab_header(title = "Highest Team OPS+ by Year") |>
tab_stubhead(label="Year") |>
cols_label(teamID = "Team", OPSP = "OPS+")
tableSub
#Export for Problem 2
write.csv(teamSub,
"C:\\Users\\Anastashia Pelletier\\Documents\\MATH 494\\teamSub.csv",
row.names = TRUE)
