add PCA R code