TeamStats

library(Lahman)
library(tidyverse)
library(plotly)
data(Pitching)
data("Teams")
data("Batting")

df <- Teams %>% 
  filter(yearID >= 2000) %>% 
  mutate(BA = H/AB, 
         singles = H - (X2B+X3B+HR), 
         OBP=(H+BB+HBP)/(AB+BB+SF+HBP), 
         SLG=((1*singles)+(2*X2B)+(3*X3B)+(4*HR))/(AB), 
         OPS = OBP+SLG,
         pOPS = (3*OBP)+SLG,
         winPerc = W/G)