add R script for correlation statistics calculation