Adding principal component analysis code