R devel - May 2014 - A bug in princomp(), perhaps?

Hi,
It may be my misunderstanding, but it seems that the "na.action" in
the princomp() function for principal components analysis does not work.  Please
see this simple example:

u <- matrix(rnorm(75), ncol=1)
v <- matrix(rnorm(20), ncol=1)
x <- u%*%t(v) + matrix(rnorm(20*75),ncol=20)
x[1,1] <- NA
pc.out <- princomp(x, na.action=na.exclude)
Error in cov.wt(z) : 'x' must contain finite values
only
>
Note, I have:
> options("na.action")
$na.action
[1] "na.omit"

Thanks,
Ravi

	[[alternative HTML version deleted]]

It's only documented to work for princomp.formula; other methods do not know
about na.action.

-pd

On 29 May 2014, at 22:10 , Ravi Varadhan <ravi.varadhan at jhu.edu> wrote:
> Hi,
> It may be my misunderstanding, but it seems that the "na.action"
in the princomp() function for principal components analysis does not work. 
Please see this simple example:
> 
> u <- matrix(rnorm(75), ncol=1)
> v <- matrix(rnorm(20), ncol=1)
> x <- u%*%t(v) + matrix(rnorm(20*75),ncol=20)
> x[1,1] <- NA
> pc.out <- princomp(x, na.action=na.exclude)
> Error in cov.wt(z) : 'x' must contain finite values only
>> 
> 
> Note, I have:
>> options("na.action")
> $na.action
> [1] "na.omit"
> 
> Thanks,
> Ravi
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
-- 
Peter Dalgaard, Professor,
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: pd.mes at cbs.dk  Priv: PDalgd at gmail.com