[Rd] A bug in princomp(), perhaps?

[peter dalgaard]

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