R script (13/23)

R script