R script (12/23)

R script