setwd('~/mlb')
require("lme4")
require("dplyr")
schmidt <- "schmj001"
pedro <- "martp001"

yr <- 2004
facs <- c("inning", "catcher", "outs_ct", 
          "bat_home_id","umpire","year_id", 
          "pitcher", "bats", "batter", 
          "stadium", "role", "start_bases_cd")

ifile <- sprintf("draFiles/dra.in.%d.csv", yr)
din.dra <- read.csv(ifile, stringsAsFactors=TRUE)
for ( fac in facs) {
    din.dra[fac] <- as.factor(din.dra[[fac]])
}

din.dra$lwts <- (din.dra$woba_pts-din.dra$woba_mean)/din.dra$wobascale
din.dra$temp_log <- log((din.dra$temp-32)*5/9 + 273) # log of temperture in Kelvin


load(sprintf("csaaFiles/csaa.%d.model.output.R", yr))

din.dra.csaa <- inner_join(din.dra, df.csaa.full, by='catcher')

#value.model <- lmer(lwts ~ bats + inning*score_diff + start_bases_cd*outs_ct + csaa + temp + stadium + PF + role + fraa*bat_home_id + inning*bat_home_id + (1|batter) + (1|pitcher) + (1|catcher) + (1|umpire), data=din)

summary(din.dra.csaa)

value.dra.full <- lmer(lwts ~ inning*score_diff + 
                   start_bases_cd*outs_ct + 
                   csaa + 
                   temp_log + 
                   bats*stadium + 
                   role + 
                   fraa*bat_home_id + 
                   inning*bat_home_id + 
                   (1|batter) + (1|pitcher) + (1|catcher) + (1|umpire), 
                   data=din.dra.csaa)


rr.dra <- ranef(value.dra.full)

value.dra.0 <- lmer(lwts ~ (1|pitcher), data=din.dra.csaa)
value.dra.inn.scorediff <- lmer(lwts ~ (1|pitcher) + inning*score_diff, data=din.dra.csaa)
value.dra.bases.outs <- lmer(lwts ~ (1|pitcher) + start_bases_cd*outs_ct, data=din.dra.csaa)
value.dra.csaa <-  lmer(lwts ~ (1|pitcher) + csaa, data=din.dra.csaa)
value.dra.temp_log <-  lmer(lwts ~ (1|pitcher) + temp_log, data=din.dra.csaa)
value.dra.bats.stadium <-  lmer(lwts ~ (1|pitcher) + bats*stadium, data=din.dra.csaa)
value.dra.role <-  lmer(lwts ~ (1|pitcher) + role, data=din.dra.csaa)
value.dra.fraa.bat_home_id <-  lmer(lwts ~ (1|pitcher) + fraa*bat_home_id, data=din.dra.csaa)
value.dra.inn.bat_home_id <-  lmer(lwts ~ (1|pitcher) + inning*bat_home_id, data=din.dra.csaa)
value.dra.batter <-  lmer(lwts ~ (1|pitcher) + (1|batter), data=din.dra.csaa)
value.dra.catcher <-  lmer(lwts ~ (1|pitcher) + (1|catcher), data=din.dra.csaa)
value.dra.umpire <-  lmer(lwts ~ (1|pitcher) + (1|umpire), data=din.dra.csaa)

save(value.dra.full
     , value.dra.0
     , value.dra.inn.scorediff
     , value.dra.bases.outs
     , value.dra.csaa
     , value.dra.temp_log
     , value.dra.bats.stadium
     , value.dra.role
     , value.dra.fraa.bat_home_id
     , value.dra.inn.bat_home_id
     , value.dra.batter
     , value.dra.catcher
     , value.dra.umpire
     , file=sprintf("dra.%d.model.output.R", yr)
     )