R script (10/23)

R script