PCA.jl

""" 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