similar to: body of non-visible function

Hi. I have been looking at the PCAproj function in package pcaPP (R 2.4.1) for robust principal components, and I'm trying to interpret the results. I started with a data matrix of dimensions RxC (R is the number of rows / observations, C the number of columns / variables). PCAproj returns a list of class princomp, similar to the output of the function princomp. In a case where I can

I'm attempting to modify how biplot draws its red vectors (among other things). This is how I've started: Biplot <- function(xx, comps = c(1, 2), cex = c(.6, .4)) { ## Purpose: Makes a biplot with princomp() object to not show arrows ## ---------------------------------------------------------------------- ## Arguments: xx is an object made using princomp() ##

Hi all, I have a difficulty to calculate the PCA scores. The PCA scores I calculated doesn't match with the scores generated by R, mypca<-princomp(mymatrix, cor=T) myscore<-as.matrix(mymatrix)%*%as.matrix(mypca$loadings) Does anybody know how the mypca$scores were calculated? Is my formula not correct? Thanks a lot! Phoebe [[alternative HTML version deleted]]

Dear All, It might be a simple question. But I could not find the answer from function “prcomp” or “princomp”. Does anyone know what are the codes to get coefficient and scores of Principal component analysis in R? Your reply will be appreciated! Best Zunqiu [[alternative HTML version deleted]]

R version: 1.7.1 OS: Red Hat Linux 7.2 When "covmat" is supplied in princomp(), the output value "center" is all NA's, even though the input matrix was indeed centered. I haven't read anything about this in the help file for princomp(). See code below for an example: pc2$center is all NA's. Jerome Asselin x <- rnorm(6) y <- rnorm(6) X <- cbind(x,y)

I start a clean session of R 1.9.1 on Windows and I run the following code: > library(MASS) > data(painters) > pca.painters <- princomp(painters[ ,1:4]) > loadings(pca.painters) Loadings: Comp.1 Comp.2 Comp.3 Comp.4 Composition 0.484 -0.376 0.784 -0.101 Drawing 0.424 0.187 -0.280 -0.841 Colour -0.381 -0.845 -0.211 -0.310 Expression 0.664 -0.330 -0.513

Dear List, How can I convert a list with elements being character strings, like: "c(1,2,3,4)", ?c(1,3,4,2) ? to a list with elements as numerical vectors: c(1,2,3,4), c(1,3,4,2)?? Thanks! Anna

Hi all, I am using the vegan package to run a prcincipal components analysis on forest structural variables (tree density, basal area, average height, regeneration density) in R. However, I could not find out how to extract factor loadings (correlations of each variable with each pca axis), as is straightforwar in princomp. Do anyone know how to do that? Moreover, do anyone knows

Dear all, i wonder if there?s a command to obtain the actual values of a principal component or a factor (not as.factor, but factanal) . test=princomp(USArrests, cor = TRUE) summary(test) just outputs, standard deviation, Prop of Variance and cumulative proportion of variance. test$loadings offers yet another proportion of variance scheme. why is that? Apart from that: Is there a

Hello, I am trying to run a PCA in R and I cannot get the PC scores for each of the values. Using pcX <- princomp(X) then loadings(pcX) I can get a listing of the eigenvectors but not the actual PC scores for each value in the dataset. I greatly appreciate any help anyone can offer Thanks Ken

Dear all I'm unable to find an example of extracting the rotated scores of a principal components analysis. I can do this easily for the un-rotated version. data(mtcars) .PC <- princomp(~am+carb+cyl+disp+drat+gear+hp+mpg, cor=TRUE, data=mtcars) unclass(loadings(.PC)) # component loadings summary(.PC) # proportions of variance mtcars$PC1 <- .PC$scores[,1] # extract un-rotated scores of

G'Day R users! Following an ordination using prcomp, I'd like to test which variables singnificantly contribute to a principal component. There is a method suggested by Peres-Neto and al. 2003. Ecology 84:2347-2363 called "bootstrapped eigenvector". It was asked for that in this forum in January 2005 by J?r?me Lema?tre: "1) Resample 1000 times with replacement entire

Hi all! I am working with R package cluster and I have a little problem: let's say I have two datasets...first one ("A") is divided into 4 clusters by means of Pam algorythm. Let's say I want to project the second database ("B") onto the Comp.1 X Comp.2 graph, and see where its elements are placed. The two datasets are made of different dim (54x19 and 28x19). I tried

Hi everyone, I am relatively new to R, and I need to perform the principal components analysis of a data matrix. I know that there are a bunch of methods to do it (dudi.pca, princomp, prcomp...) but I have not managed to find a method that can return the residuals obtained by retaining X principal components of the original data, as this MATLAB function can do: http://is.gd/6WeUFF Suggestions?

Hi everyone, I am trying to get the factor score for each individual case from a principal component analysis, as I understand, both princomp() and prcomp() can not produce this factor score, the principal() in psych package has this option: scores=T, but after running the code, I could not figure out how to show the factor score results. Here is my code, could anyone give me some advice please?

Dear all, I'm using R 2.8.1 under Vista. I programmed a Simulation with the code enclosed at the end of the eMail. After the simulation I want to analyse the columns of the single simulation-runs, i.e. e.g. Simulation[[1]][,1] sth. like that but I cannot address these columns... Can anybody please help? Best, Thomas ############################ CODE ############################

Hi- I'm trying to make my PCA (princomp) colored. In my csv excel sheet, I have the first column numbered according to the groupings I want to assign to the PCA. I've played around with trying to set this first column as the color vector, but haven't had any luck. Any suggestions? Thanks, Hillary [[alternative HTML version deleted]]

Hello, Could somebody point me, is there any function in R which returns expected values of order statistics for normal distribution? I have been looking and couldn't find it. Thanks! Anna

Hi All, There were lot of diffrences in the R and SPSS results for Exploratory Factor Analysis.why is it so ?I used standard factor analysis functions like:-- factanal(m, factors=3, rotation="varimax") princomp(m, cor = FALSE, scores = TRUE, subset = rep(TRUE, nrow(as.matrix(m)))) print(summary(princomp(m, cor=TRUE),loadings = TRUE, cutoff = 0.2), digits = 2) prcomp(m, scale = TRUE)

$ R --version R 1.6.1 (2002-11-01). So I would like to perform principal components analysis on a 16X16 correlation matrix, [princomp(cov.mat=x) where x is correlation matrix], the problem is princomp complains that it is not non-negative definite. I called eigen() on the correlation matrix and found that one of the eigenvectors is close to zero & negative (-0.001832311). Is there any way