[Rd] prcomp need Conj(t(s$vt)) (PR#2924)
pgilbert at bank-banque-canada.ca
pgilbert at bank-banque-canada.ca
Mon May 5 21:14:23 MEST 2003
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
More information about the R-devel
mailing list