add the R script used for the analysis:

1 download and read the raw data