Change pca calculation by svd to avoid covariance calculation. Use ml covariance estimator.

These changes were motivated by the question of Julien Delafontaine on the julia-stats mailing list.