yr <- 2004 setwd('~/mlb') require("dplyr") require("lme4") require('ggplot2') schmidt <- "schmj001" pedro <- "martp001" bigunit <- "johnr005" load(sprintf('draFiles/dra.%d.model.output.R', yr)) din <- read.csv(sprintf('draFiles/dra.in.%d.csv', yr), stringsAsFactors=T) objects() summary(din) value.dra.full df.small <- data.frame(sapply(1:nrow(din), function(x) {return(1)})) names(df.small) <- "npa" df.small$npa <- as.integer(df.small$npa) df.small$pitcher <- din$pitcher head(df.small) df.agg <- aggregate(npa ~ pitcher, df.small, sum) df.agg$year <- yr head(df.agg) objects() rr.dra.0 <- ranef(value.dra.0)$pitcher rr.dra.bases.outs <- ranef(value.dra.bases.outs)$pitcher rr.dra.bats.stadium <- ranef(value.dra.bats.stadium)$pitcher rr.dra.batter <- ranef(value.dra.batter)$pitcher rr.dra.catcher <- ranef(value.dra.catcher)$pitcher rr.dra.csaa <- ranef(value.dra.csaa)$pitcher rr.dra.fraa.bat_home_id <- ranef(value.dra.fraa.bat_home_id)$pitcher rr.dra.full <- ranef(value.dra.full)$pitcher rr.dra.inn.bat_home_id <- ranef(value.dra.inn.bat_home_id)$pitcher rr.dra.inn.scorediff <- ranef(value.dra.inn.scorediff)$pitcher rr.dra.role <- ranef(value.dra.role)$pitcher rr.dra.temp_log <- ranef(value.dra.temp_log)$pitcher rr.dra.umpire <- ranef(value.dra.umpire)$pitcher names(rr.dra.0) <- "icpt.dra.0" names(rr.dra.bases.outs) <- "icpt.dra.bases.outs" names(rr.dra.bats.stadium) <- "icpt.dra.bats.stadium" names(rr.dra.batter) <- "icpt.dra.batter" names(rr.dra.catcher) <- "icpt.dra.catcher" names(rr.dra.csaa) <- "icpt.dra.csaa" names(rr.dra.fraa.bat_home_id) <- "icpt.dra.fraa.bat_home_id" names(rr.dra.full) <- "icpt.dra.full" names(rr.dra.inn.bat_home_id) <- "icpt.dra.inn.bat_home_id" names(rr.dra.inn.scorediff) <- "icpt.dra.inn.scorediff" names(rr.dra.role) <- "icpt.dra.role" names(rr.dra.temp_log) <- "icpt.dra.temp_log" names(rr.dra.umpire) <- "icpt.dra.umpire" ndf <- data.frame(cbind(df.agg, rr.dra.0, rr.dra.bases.outs, rr.dra.bats.stadium, rr.dra.batter, rr.dra.catcher, rr.dra.csaa, rr.dra.fraa.bat_home_id, rr.dra.full, rr.dra.inn.bat_home_id, rr.dra.inn.scorediff, rr.dra.role, rr.dra.temp_log, rr.dra.umpire)) head(ndf) write.table(ndf, file = sprintf("dra.value.%d.out.csv", yr), row.names=FALSE, na="",col.names=TRUE, sep=",")