(PR#6452) Re: [Rd] predict.princomp fails if princomp was called
Yes, and it should fail or return NA, as the required information is not available. The problem is that it had logical NAs not numeric NAs for the center component. #3682 points out something that is on the help page: if you supply covmat as a matrix (and not a covariance list) there is no information on the centring done and x is ignored.
On Thu, 22 Jan 2004 joehl@gmx.de wrote:
Dear all, predict.princomp seems to expect non-NA component $center which is not set by princomp.default. Best regards Jens Oehlschl?gel
x <- matrix(rnorm(400), ncol=4) predict(princomp(x, covmat=cov(x)), newdata=x[1:10, ])
Error in if (center) { : missing value where TRUE/FALSE needed
In addition: Warning message:
the condition has length > 1 and only the first element will be used in: if
(center) {
version
_ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 8.1 year 2003 month 11 day 21 language R
Brian D. Ripley, ripley@stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595