R script (21/23)

R script