added pca_analysis.R