similar to: temporal evolution and variance after rotation of eof

Dear R users, I have been doing rotation on Principal components analyse, with varimax function and promax. Following those changes, I cannot find now how to get the spatial evolution and the variance. Indeed, with the function princomp that has been used to get the principal components, arguments such as $scores and $sdev were available to get the spatial evolution and the variance, with

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

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

DiBcom (http://www.dibcom.com) is a company designing chipsets in the field of digital TV reception. The "DSP" team, in charge of the LLVM/CLANG based compilation chain, as well as the cpu design, is hiring. Beside participating to the toolchain development and deployment, you will also be part of the improvement process, by proposing efficient coding style for the dsp firmware or

Le 29/11/2017 ? 08:26, Sorin Srbu a ?crit?: > If updated, the instrument software will break, like in an atomic mushroom > cloud, rendering critical hardware non-working and lab people standing > outside my office with torches, pitch-forks and all the shebang. > Yupp, unfortunately that's the current state of some lab equipment > manufacturers that use RHEL as a base...