Add PCA scripts from Mingfeng