[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