search for: varimax

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 factanal function. Does anybody know which are the instructions into R to do exactly the same of what I h...

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: http://r.78...

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

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 help...

...ce be due to? Rounding in the initial data? Or implementation details of the R oblimin function in GPArotation? Jackson writes about 'raw oblimin', 'normal oblimin' and 'direct oblimin' and I do not know how that relates to the R oblimin implementation. I have also tried varimax on data and results given by Mardia in his 'Multivariate analysis' book Table 9.4.1. Mardia uses the communalities from the factor analysis in the expression for the varimax rotation. I dont see how the R varimax function can handle the communalities. I dont have the book right at hand,...

I am using R mainly for multiple linear regression and principal components analyses and I am quite happy with it. I think it is a worderful work you have done. But I am having a problem when using Varimax for rotating loadings obtaines from a princomp (I use the cor = TRUE option): the variance across the components in the variables does not mantain. Let me explain it with an example in which the variances, which should be the same, are listed in both the unrotated and rotated models: Unrotated mo...

Hello, recently I was looking at the implementation of the "varimax" rotation procedure from the "stats" package and to me it looks quite different from the algorithm originally suggested by Kaiser in 1958. The R procedure iteratively uses singular value decompositions of some matrices whereas Kaiser proposed to iteratively compute rotation matr...

...365 -0.2779157 DS140_06 0.2928942 0.1218132 DS140_07 0.2601528 -0.1162116 DS140_08 0.2737338 0.2811998 DS140_09 0.3012719 -0.1714994 DS140_10 0.2653410 0.3159160 DS140_11 0.2590944 0.2347922 DS140_12 0.2837112 -0.2653533 DS140_13 0.3246268 -0.2187217 DS140_14 0.2896170 0.2190227 > varimax(pc$loadings[,1:2]) $loadings Loadings: Comp.1 Comp.2 DS140_01 -0.424 DS140_02 0.390 DS140_03 0.146 -0.384 DS140_04 0.375 DS140_05 0.395 DS140_06 0.143 0.283 DS140_07 0.273 DS140_08 0.392 DS140_09 0.340 DS140_10 0.413 DS140_11...

Hi: How i can make Varimax and Promax Rotation in R 0.63 thanks Jorge M. A. Magalhães -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the &q...

Dear R users the factanal pacakge is always MLE, which package can do varimax rotation with the results from princomp ? thank you yong

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,...

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 stu...

Is it possible to retrieve sums of squared loadings after applying varimax rotation? Here's the setup to my problem: I ran PCA using prcomp(). I then applied the Kaiser criterion to retain only the components having eigenvalues >= 1. (I know there's debate about the wisdom of that criterion, but I don't want to get sucked into that.) I then fed the reduce...

I was wondering whether the following inconsistency of the promax rotation function with the results of a promax rotation using SAS should be considered a bug in the promax function of R. Any comments will be highly appreciated. The following is a loading matrix obtained from a varimax rotation in SAS: # Factor loadings after varimax rotation x <- t(array(c(0.78107, 0.35573, 0.85995, 0.07887, 0.79376, 0.22499, 0.26307, 0.85087, 0.29449, 0.72966, 0.06064, 0...

...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? Thank you very much. > pc <- principal(a,rotate="varimax",scores=TRUE) > pc Principal Components Analysis Call: principal(r = a, rotate = "varimax", scores = TRUE) Standardized loadings (pattern matrix) based upon correlation matrix PC1 h2 u2 V1 0.80 6.4e-01 0.36 V2 0.03 7.9e-04 1.00 V3 -0.92 8.4e-01 0.16 V4 0.00 2.0e-06...

...in terms of the cumulative % of the total variance, so when i go from total solution of 8 (8 variables in the data set), to a reduced number of components, i can evaluate % of variance explained, or am I missing something?? 8 variables in the data set > princ = prcomp(df[,-1],rotate="varimax",scale=TRUE) > summary(princ) Importance of components: PC1 PC2 PC3 PC4 PC5 PC6 PC7 PC8 Standard deviation 1.381 1.247 1.211 0.994 0.927 0.764 0.6708 0.4366 Proportion of Variance 0.238 0.194 0.183 0.124 0.107 0.073 0.0562 0.0238 Cumulative Prop...

...363925 0.4326113 #v3 0.5128866 0.4124441 1.0000000 0.8770750 0.5128866 0.4320310 #v4 0.4320310 0.4084281 0.8770750 1.0000000 0.4320310 0.4323259 #v5 0.4664948 0.4363925 0.5128866 0.4320310 1.0000000 0.9473451 #v6 0.4086076 0.4326113 0.4320310 0.4323259 0.9473451 1.0000000 factanal(m1, factors=3) # varimax is the default # Output from R #Call: #factanal(x = m1, factors = 3) #Uniquenesses: # v1 v2 v3 v4 v5 v6 #0.005 0.101 0.005 0.224 0.084 0.005 #Loadings: # Factor1 Factor2 Factor3 #v1 0.944 0.182 0.267 #v2 0.905 0.235 0.159 #v3 0.236 0.210 0.946 #v4 0.180 0.242...

...ch PC prop = pca$sdev^2/sum(pca$sdev^2) # cumulative proportion variance explained by each PC cumProp = cumsum(prop) # following print statements also can be obtained # from print(summary(pca)) #print(prop) #print(cumProp) print(summary(pca)) # Rotate the PCA loadings through PCs 1 and 2 using VARIMAX rotation rot = varimax(pca$rotation[,1:2], normalize = TRUE, eps = 1e-5) < end code > How can I calculate the new % variance explained by each PC after rotation ?????? Many thanks once more, -- Maurice J. McHugh, Ph.D. Assistant Professor Department of Geography and Anthropology 227...

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 US...

Dear Folks I need to rotate PCA loadings and scores using R. I have run a pca using princomp and I have rotated PCA results with varimax. Using varimax R gives me back just rotated PC loadings without rotated PC scores. Does anybody know how I can obtain/calculate rotated PC scores with R? Your kindly help is appreciated in advance Francesca [[alternative HTML version deleted]]