clear
clear matrix
set matsize 800
set mem 500m
cd [your directory]


**team w/l records and runs scored

use [team_w_l_r_ra.csv]


*** merge team-level starter and reliever runs saved 


    merge 1:1 teamID yearID using [efip_team_starters_relievers_rs.csv],nogen
	merge 1:1 teamID yearID using [fg_team_starter_value.csv,nogen


    keep yearID teamID G W L R ip prsRP z_prsSP z_fgsrar

*** create season standardized runs scoredd

bys yearID: egen z_runs=std(R)


*** create season-standardized reliever runs saved per inning
** don't use reliever runs saved: endogenous to starter performance

    gen prip=prsRP/ip



* 1. Summarize ip within each year
egen total_ip = sum(ip), by(year)

* 2. Compute weighted mean of prip
generate weighted_prip_product = prip * ip
egen total_weighted_prip = sum(weighted_prip_product), by(year)
generate weighted_mean_prip = total_weighted_prip / total_ip

* 3. Compute weighted sum of squared deviations
*    (prip - weighted_mean_prip)^2 * ip
generate wssd = (prip - weighted_mean_prip)^2 * ip
egen total_wssd = sum(wssd), by(year)

* 4. population variance: divide by total weight
generate popvar_prip = total_wssd / total_ip

* 5. Population weighted standard deviation
generate ppripd_prip = sqrt(popvar_prip)


drop wssd weighted_prip_product total_weighted_prip


rename weighted_mean_prip wmu_prip
rename ppripd_prip wsd_prip
drop total_wssd

gen z_prip=(prip-wmu_prip)/wsd_prip


*** merge gridwar 

merge 1:1 yearID teamID using  [gridwar_team_coded.csv; ceated with GWAR_raw_code_prep.txt],nogen

*** merge team fielding runs saved ****
merge m:1 teamID yearID using [team_fielding_rs.csv], nogen

drop if year<1900



gen wp=W/G // winning percentage outcome variable


**** nested regressions

nestreg: regress wp (z_runs z_field_composite z_prip) (z_gwar) [aweight=G] if year>1951

nestreg: regress wp (z_runs z_field_composite z_prip) (z_fgsrar) [aweight=G] if year>1951

nestreg: regress wp (z_runs z_field_composite z_prip) (z_prsSP) [aweight=G] if year>1951

**note: gridwar only since 1952; drop yr restriction if you want to test FG & efip since 1900