R script of analysis