add R script to fit polynomial regression with data.