""" Performs PCA PCA(W)""" function PCA(W) W0 = mean(W,1); W = W-repmat(W0,size(W,1),1); (score,latent,C) = svd(W) score = score*diagm(latent) C,score,latent,W0 end