add R scripts for statistical analysis