similar to: Rotating pca scores

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

I am running a PCA, but would like to rotate my data and limit the number of factors that are analyzed. I can do this using the "principal" command from the psych package [principal(my.data, nfactors=3,rotate="varimax")], but the issue is that this does not report scores for the Principal Components the way "princomp" does. My question is: Can you get an

A bit of a newbee to R and factor rotation I am trying to understand factor rotations and their implementation in R, particularly the GPArotation library. I have tried to reproduce some of the examples that I have found, e.g., I have taken the values from Jacksons example in "Oblimin Rotation", Encyclopedia of Biostatistics

Haj i try to perform a principal component analysis by using a tetrachoric correlation matrix as data input tetra <- tetrachoric (image_na, correct=TRUE) t_matrix <- tetra$rho pca.tetra <- principal(t_matrix, nfactors = 10, n.obs = nrow(image_na), rotate="varimax", scores=TRUE) the problem i have is to compute the individual factor scores from the pca. the code runs perfect

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?

Hello all, I've run a principal component regression using the PLS package. I then applied varimax rotation (i.e., using http://stat.ethz.ch/R-manual/R-patched/library/stats/html/varimax.html). I cannot figure out how to extract the factor loadings post-varimax. Is there a command to do this? scores(x) does not do it. Thanks and happy holidays -- View this message in context:

Dear R Helper, I am writing because I seek to perform a varimax rotation on my Principal Components Analysis (PCA) solution. (I have been performing PCA's using the eigen command in R.) If you can tell me how to perform this rotation when I use the eigen command (or the princomp command) I would be thrilled. Thanks so much! Wendy Treynor Ann Arbor, MI USA

Hello, I need to do a principal component analysis with EQUAMAX-rotation. Unfortunately the function principal() I use normally for PCA does not offer this rotation specification. I could find out that this might be possible somehow with the package GPArotation but until now I could not figure out how to use this in the principal component analysis. Maybe someone can give an example on how to do

--=-hiYzUeWcRJ/+kx41aPIZ Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: 8bit Hi, on March 10 I filed a wishlist bug report asking for the inclusion of some changes to factanal() and the associated print method. The changes were originally proposed by John Fox in 2005; they make print.factanal() display factor correlations if factanal() is called with rotation =

Hello everybody, I would need some help from you. I am trying to fit a logistic model to some presence absence data of animals living on river islands. I have got 12 predictor variables and I am trying to use a stepwise forward method to fit the best logistic model to my data. I am using the function STEP (stats). I have a question for you. Can I use step function if my variables have a

hi, I have two questions: #first (SPSS vs. R): I just compared the output of different PCA routines in R (pca, prcomp, princomp) with results from SPSS. the loadings of the variables differ vastly! in SPSS the variables load constantly higher than in R. I made sure that both progr. use the correlation matrix as basis. I found the same problem with rotated values (varimax rotation and rtex=T

Hi! I am trying to analyse with R a database that I have previously analysed with SPSS. Steps with SPSS: Factorial analysis Extraction options : I select = Principal component analysis Rotation: varimax Steps with R: I have tried it with varimax function with factanal or with princomp...and the results are different of what I have with SPSS. I think that varimax function is incorporated in

In some documentation for a package I am working on I have > \code{\link[stats]{varimax}} > \code{\link[stats]{promax}} The link to varimax works, but not the one to promax. Promax is an alias under \name{varimax}. This kind of link works within a package, but I'm not sure if it is suppose to work when it is a link to another package. Is this a known limitation or bug, or

Thanks to Mr.Liviu Androvic and Mr.Richard Rowe helped me in PCA. Because I have just learn R language in a few day so I have many problem. 1) I don't know why PCA rotation function not run although I try many times. Would you please hepl me and explain how to read the PCA map (both of rotated and unrotated) in a concrete example. 2) Where I can find document relate: Plan S(A), S(A*B),

hello, subsequently to a NMDS analysis (performed with metaMDS or isoMDS) is it possible to rotate the axis through a varimax-rotation? Thanks in advance. Bernd Panassiti

I'm using R for PCA and? factor analysis. I want to create biplots of varimax rotated factors that color-code points by their classification. My research is on streams that are urban and rural. So, I want to color code them by this classification. If you just do a biplot from prcomp or princomp, you cannot add this color. So, I have used some code developed by a graduate student in our

Hi everyone- I have used several methods to calculate principal components rotated using the varimax procedure. This is simple enough. But I would like to calculate the % of variance explained associated with each PC before and after rotation. factanal returns the % of variance explained associated with each PC but I cannot seem to get it to change after rotation. Many thanks for your

Hello again- Thanks to Prof. Ripley for responding to my previous question. I would like to clarify my question using sample code. I will use some sample code taken from ?prcomp Again, I would like to compare the % variance explained by each PC before and after rotation. < code follows > data(USArrests) pca = prcomp(USArrests, scale = TRUE) # proportion variance explained by each

Dear Users, I ran factor analysis using R and SAS. However, I had different outputs from R and SAS. Why they provide different outputs? Especially, the factor loadings are different. I did real dataset(n=264), however, I had an extremely different from R and SAS. Why this things happened? Which software is correct on? Thanks in advance, - TY #R code with example data # A little

Hi, I am puzzeled with a differing result of princomp in R and FACTOR in SPSS. Regarding the amount of explained Variance, the two results are the same. However, the loadings differ substantially, in the unrotated as well as in the rotated form. In both cases correlation matrices are analyzed. The sums of the squared components is one in both programs. Maybe there is an obvious reason, but I