This is R code, I copy to here