R script (3/23)

R script