In prcomp s <- La.svd(x, nu = 0) s$v <- t(s$vt) the second above line should be s$v <- Conj(t(s$vt)) (to cover complex cases). Paul Gilbert