add principal component analysis (calcpca.m)