R script (17/23)

R script