R script (7/23)

R script