Statistics programming using R