R code from my computational statistics projects