update R code to do the PCA...