From pgilbert at bank-banque-canada.ca  Mon May  5 21:14:23 2003
From: pgilbert at bank-banque-canada.ca (pgilbert@bank-banque-canada.ca)
Date: Mon May  5 20:14:35 2003
Subject: [Rd] prcomp need Conj(t(s$vt)) (PR#2924)
Message-ID: <200305051814.h45IEN2b016635@pubhealth.ku.dk>

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


[Rd] prcomp need Conj(t(s$vt)) (PR#2924)
pgilbert@bank-banque-canada.ca