similar to: factanal

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

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

## Hello.. I have asked a similar question, but this is not fixed as before. ## I am running the following using Ubuntu OS: R version 2.13.1 (2011-07-08) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-pc-linux-gnu (64-bit) ## when I do this: > install.packages("mapproj", dependencies=T) ## I get this: Installing package(s)

Hi, I installed the maps and mapdata libraries on my R-1.8.0 on MacOSX 10.2.8 (jaguar on a powerbook G4), and failed to make the map function work properly: R : Copyright 2003, The R Development Core Team Version 1.8.0 (2003-10-08) [...] > library(maps) > map() Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file

Greetings folks, I'm currently dealing with a company to let me set Caller-ID-Name on outbound calls. So far pretty happy with their services. The basic service works like this: * CLEC sets Point Code to point to this company * CLEC has to sign LOA saying they give me permission to set the Caller-ID-Name through this company. * I go into web interface and set name. However, the CLEC is

In comparing the results of princomp and prcomp I find: 1. The reported standard deviations are similar but about 1% from each other, which seems well above round-off error. 2. princomp returns what I understand are variances and cumulative variances accounted for by each principal component which are all equal. "SS loadings" is always 1. 3. Same happens

Hi, I have a set of covariance matrices but not the original data. I want to carry out some exploratory factor analysis. So, I am trying to construct a covariance matrix list as the input for factanal. I can construct a list which includes the cov, the centers, and the n.obs. But it doesn't work. I get an error that says "Error in sqrt(diag(cv)) : Non-numeric argument to mathematical

Dear friends, I'm using R to produce the following Factor Analysis: > matriz.cor<-hetcor(matrix(as.factor(data), ncol=variables, byrow=T))$correlations > factanal(x=data, factors=2, covmat=matriz.cor, scores='regression') Then the screen output shows the following message: Error en factanal(x = data, factors = 2, covmat = matrix, : requested scores without

Hi everyone, I'm running a factor analysis on a correlation matrix with 32 rows and columns. I get the following error when I issue the command sequence mich.fac1 <- factanal(michcor, factor=1) Error in solve.default(cv) : system is computationally singular: reciprocal condition number = 3.24729e-18 I'd really appreciate an explanation for this error and a solution to the problem

I understand how to control the first local channel, but an having trouble getting the second local channel to enter stasis. I setup have the following extensions.conf to handle 1000 (basically had it setup so if first stasis not there try second, but believe second channel never processes the dial plan so even if second line was hello-world2 it would not matter. [mycontext] exten =>

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

In R factanal output, I can't find a function to give me residuals e. I mannually got it by using x -lamda1*f1 -lamda2*f2 - ... -lamdan*fn, but the e I got are not uncorrelated with all the f's. What did I do wrong? Please help. Yijun ____________________________________________________________________________________ Be a better friend, newshound, and

Alle, I'm following the instructions in section 19.3.3.2 of the docs @ http://www.centos.org/docs/5/html/5.1/Deployment_Guide/s2-nfs-config-autofs-LDAP.html, but I cannot add the following entry in LDAP: dn: automountMapName=auto.home,dc=subaru,dc=nao,dc=ac,dc=jp objectClass: top objectClass: automountMap automountMapName: auto.home After looking at the schemas in /etc/openldap/schema,

What is the deal with "CSR's"? TWTelecom is telling me that I can't port a number to their service without a Customer Service Record. Apparently this is easy with Verizon, and not so easy with some other companies. Basically I'm at a brick wall with a couple of ports because TWTelecom is telling me I HAVE to get a CSR and certain other providers (Time Warner Cable for

I need some help understanding how on of the example data sets is formatted in the basic R installation. If I load the Mona Loa CO2 data, with the command: > data(co2) I can view the data with: > co2 And the data are in the form of 11 rows labeled as years (1994-2004) and 12 columns labeled (Jan - Dec). This structure appears to be a dataframe, however, if I type the command

Hallo, can I use pxelinux for booting a notebook via WLAN? Viele Gruesse! Helmut

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,