similar to: spatial evolution and variance after rotation of Principal components

Dear R users, I have been doing some eof analysis using princomp function, then the eof results were rotated with varimax and promax functions. Those functions are working fine. However, after their uses, I cannot find how to obtain the temporal evolution and the variance of the obtain components, as varimax and promax does not have the arguments $scores and $sdev. If any body, know how to

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,

# I have a question about how factanal is calculating the regression factor # scores based on an oblique rotation (promax) of the factors. # # As is explained in the help file, regression factor scores are # obtained as # # hat f = Lambda' Sigma^-1 x # # However, according to Harman's "Modern Factor Analysis" (e.g. second # edition, pp. 351-352) the formula is # # hat f = Phi

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

Hi, everyone, This is a little silly, but I cant figure out the algorithm behind lm.fit function used in the context of promax rotation algorithm: The promax function is: promax <- function(x, m = 4) { if(ncol(x) < 2) return(x) dn <- dimnames(x) xx <- varimax(x) x <- xx$loadings Q <- x * abs(x)^(m-1) U <- lm.fit(x, Q)$coefficients d <-

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 "body", not the subject !) To: r-help-request at

Hello I seem to find only two types of rotation for the factanal function in R, the Varimax and Promax, but is it possible to run a orthogonal and oblique rotations in R? Thanks in advance Rosario

Hello all, I'm trying to understand when exactly the variable substitution happens in Dovecot, as I get different results in the LDA and in the POP/IMAP server. Here's my general setup: - Postfix MTA, delivers mail via Dovecot LDA - users list sits in a database - mail directory is stored in the database, but with Dovecot variables inside Delivery works perfectly fine,

Hi, I found a discrepancy between results in R and Stata for a factor analysis with a promax rotation. For Stata: . *rotate, factor(2) promax* (promax rotation) Rotated Factor Loadings Variable | 1 2 Uniqueness -------------+-------------------------------- pfq_amanag~y | -0.17802 0.64161 0.70698 pfq_bwalk_~ø | 0.72569 0.05570

--=-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 =

Hi list, (Please Cc: me in your replies because I'm not subscribed to this list.) While trying to build lukemftpd staticaly on FreeBSD, I got a link-time error. Libssh.a indeed provides the "xmalloc" symbol (I suppose there are more). I wonder if this is whether intentional or not. It is a very common function name, and I think it would be worth renaming it to something like

Full_Name: conor dolan Version: 1.4.0 OS: windows98 Submission from: (NULL) (146.50.170.247) On the basis of Lawley and Maxwell's explanation in Factor Analysis as a Statistical Method (London: Butterworth, 1971), the following line in the promax routine (library, mva): attr(z, "covariance") <- crossprod(U) should be attr(z, "covariance") <- solve(crossprod(U))

Hi, In a discussion of factor analysis in "Using Multivariate Statistics" by Tabachnick and Fidell, two matrices are singled out as important for interpreting an exploratory factor analysis (EFA) with an oblique promax rotation. One is the "structure matrix". The structure matrix contains the correlations between variables and factors. However, these correlations may be

FATNA FATTACH 23 ANS Hay El FALAH RUE 1 N?56 RESIDENCE BEN AMOUR 2 ETAGE 2 Casablanca 0653185512 E-mail : jalilafattach at hotmail . Fr. TECHNICIEN SPECIALISE EN RESEAUX INFORMATIQUE ETUDES ET FORMATION : -2013-2014 : Dipl?me en technicien sp?cialis? en r?seaux informatiques. -Niveau d?instruction : Bac science de la vie

I am performing a large (105 variable) factor analysis with factanal, specifying promax rotation. I kow that some loadings over 1.0 are not unsual with that rotation, but I have some as large as 1.2, which seems extreme. I am skirting the assumptions of the model by using responses on a 7-point rating scale as data; I may have to go back and compute polychoric correlations instead of product

Dear R-devel list members, Ben Fairbank draw it to my attention that factanal() (in the stats package) doesn't report factor correlations for oblique rotations. Looking at the source, I see that factanal also doesn't save the factor-transformation (rotation) matrix from which these correlations can be computed. I've modified the source, attached below, so that the transformation

I have data set with binary responses. I would like to conduct polychoric principal component analysis (pPCA). I know there are several packages used in PCA but I could not find one that directly estimate pPCA and graph the individuals and variables maps. I will appreciate any help that expand these reproducible scripts. #How to conduct polychoric principal component analysis pPCA using #either

Hi all, I need to compute some "occurence matrix": given a zero matrix and a set of paired indexes, I want to store the number of occurences of each paired index in a matrix. The paired indexes are stores as an index matrix. I prefere not to use loops for performances purpose. Here follows a dummy example: > occurence <- matrix(0, 2, 2); data [,1] [,2] [1,] 0 0

Dear Tobias, I finally succeeded in exporting my dataframes from R into Access with this script: library(RODBC) canal <- odbcConnectAccess("D:/Floristique.mdb") sqlSave(channel=canal, dat=flore, tablename="Floristique", rownames=F, safer=F, fast=F, varTypes=c(dates="Date")) odbcClose(canal) My problem in exporting my dataframe "flore" seems to

hello all, I would like to introduce a summary table into the pdf along with the plots (in order to archive my data into single files automatically). Similarly, It would be great to have the result of the statistical analysis (for instance anova) in the same file. Is there a way to do that? example: pdf("example.pdf") layout (matrix(1:2,1,2)) plot (groups, scores) Result <-