R script (8/23)

R script