PCA main code