```
library(data.table)
setwd(dirname(parent.frame(2)$ofile))
set.seed(2716) # Keep at 2716. Or else HRs will be different every time.
# ---------------------------------------------------------------------------------|
# make new heat rates
# ---------------------------------------------------------------------------------|
final.fits = fread('Final_Fits.csv')
# calculate average heat rate at min stable level and three marginal heat rate points thereafter
final.fits[,load_33:=load_min + (1/3)*(load_max-load_min)]
final.fits[,load_66:=load_min + (2/3)*(load_max-load_min)]
final.fits[,HR_avg_0:=a0+a1*load_min+a2*load_min^2+a3*load_min^3+a4*load_min^4]
final.fits[,HR_incr_1:=(load_33*(a0+a1*load_33+a2*load_33^2+a3*load_33^3+a4*load_33^4)-
load_min*(a0+a1*load_min+a2*load_min^2+a3*load_min^3+a4*load_min^4))/
(load_33-load_min)]
```

This file has been truncated. show original