R script (9/23)

R script