My R programming code