added R code for hill function plots