search for: prcomponents

The call to sweep in this function which was working in 0.62.2 is giving me trouble in 62.3: prcomponents <- function(x, center=TRUE, scale=TRUE, N=nrow(x)-1) {if (center) center <- apply(x,2,mean) else center <- rep(0, ncol(x)) if (scale) scale <- sqrt(apply(x,2,var)) else scale <- rep(1, ncol(x)) s <- svd(sweep(sweep(x,2, center),2, scale, FUN=&qu...

...oth Splus and R) and perhaps someone else would like to work on it. Of all the changes below, the only ones I consider really important are the documentation and the new version of prcomp which gives results like Splus. I have some misgivings about adding yet another principal components function (prcomponents below). However, I would like prcomp and princomp to return the same results as Splus, and yet I see that there is need for improvement. I would like this improvement to happen in a function which has a different name from the functions in Splus, and perhaps be moved into an expanded "compatib...