similar to: factanal: error message

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

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

Greetings, all: I am using factanal in R. When I enter a matrix or a formula, the print method winds up with something like this: Test of the hypothesis that 6 factors are sufficient. The chi square statistic is 28.1 on 22 degrees of freedom. The p-value is 0.172 But when I enter a covmat, the print method winds up with something like this: The degrees of freedom for the model is 22 and the

Dear list members, Could any expert on factor analysis be so kind to explain how to calculate AIC on the output of factanal. Do I calculate AIC wrong or is factanal$criteria["objective"] not a negative log-likelihood? Best regards Jens Oehlschl?gel The AIC calculated using summary.factanal below don't appear correct to me: n items factors total.df rest.df model.df

Dear list members, Could any expert on factor analysis be so kind to explain how to calculate AIC on the output of factanal. Do I calculate AIC wrong or is factanal$criteria["objective"] not a negative log-likelihood? Best regards Jens Oehlschl?gel The AIC calculated using summary.factanal below don't appear correct to me: n items factors total.df rest.df model.df

Hi, I am currently trying to execute the following command: f<-factanal(factors=k$Components$nparallel,covmat=m,n.obs=2287,rotation="varimax",control=list(opt=list(method=c("BFGS")))) but keep getting the error: L-BFGS-B needs finite values of 'fn' I can't figure out what I am doing wrong here, why isn't optim being told to use BFGS instead of L-BFGS-B...

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

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: https://stat.ethz.ch/pipermail/r-help/attachments/20070513/f07553e9/attachment.pl

Can someone please direct me to how to run a factor analysis in R by first inputting a correlation matrix? Does the function "factanal" allow one to read a correlation matrix instead of data vectors? Thanks, Matt. [[alternative HTML version deleted]]

Hi, I had a look at the MASS4 scripts in the MASS package, in Ch 11.3 Factor Analysis, there is a section of codes like: data(ability.cov) ability.FA <- factanal(covmat = ability.cov, factors = 1) ability.FA (ability.FA <- update(ability.FA, factors = 2)) #summary(ability.FA) round(loadings(ability.FA) %*% t(loadings(ability.FA)) + diag(ability.FA$uniq), 3)

I try to code the ULS factor analysis descrbied in ftp://ftp.spss.com/pub/spss/statistics/spss/algorithms/ factor.pdf # see PP5-6 factanal.fit.uls <- function(cmat, factors, start=NULL, lower = 0.005, control = NULL, ...) { FAfn <- function(Psi, S, q) { Sstar <- S - diag(Psi) E <- eigen(Sstar, symmetric = TRUE, only.values = TRUE) e <- E$values[-(1:q)] e <-

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

Dear all, Has someone implemented in R (or any other language) Knol DL, ten Berge JMF. Least-squares approximation of an improper correlation matrix by a proper one. Psychometrika, 1989, 54, 53-61. or any other similar algorithm? Best regards Jens Oehlschl?gel Background: I want to factanal() matrices of polychoric correlations which have negative eigenvalue. I coded Highham 2002

Hi, I wish to compute multivariate test statistics for a within-subjects repeated measures design with anova.mlm. This works great if I only have two factors, but I don't know how to compute interactions with more than two factors. I suspect, I have to create a new "grouping" factor and then test with this factor to get these interactions (as it is hinted in R News 2007/2), but

Dear all, The 'make check' step fails for the pacakge mva on IBM AIX. The tail of the Rout log file looks like: > for(factors in 2:4) print(update(Harman23.FA, factors = factors)) Call: factanal(factors = factors, covmat = Harman23.cor) Uniquenesses: height arm.span forearm lower.leg weight 0.170 0.107 0.166

I am having difficulty obtaining the scores from my principal component analysis. I have used this method before and have had no problems. The data set that I am using this time is similar to what I have used in the past. What do I need to do to my dataset in order for me to obtain these scores? R screen says the following message Error in factanal(covmat = pasa.cov, factors = 4) :

Hi! Suppose I have a factor: > fac1 <- factor(rep(c("dog","cat","tree"),c(2,3,4))) > fac1 [1] dog dog cat cat cat tree tree tree tree Levels: cat dog tree > length(fac1) [1] 9 Now I want to get rid of the dogs: > fac1 <- fac1[3:9] > fac1 [1] cat cat cat tree tree tree tree Levels: cat dog tree > How should I remove the dog from the

Let's say I have a within-subject experiment with 2 observables, obs1 and ob2 and 2 independent factors, fac1 and fac2. I can do summary( aov( obs1~fac1*fac2 + Error(Subject/(fac1*fac2)) ) ) summary( aov( obs2~fac1*fac2 + Error(Subject/(fac1*fac2)) ) ) to test the 2 observables separately. > summary( fit<-manova( cbind(obs1,obs2)~fac1*fac2 + Error(Subject/(fac1*fac2)) ) ) gives

Hey, hopefully there is an easy way to solve my problem. All that i think off is lengthy and clumsy. Given a data.frame d with columns VALUE, FAC1, FAC2, FAC3. Let FAC1 be something like experiment number, so that there are exactly the same number of rows for each level of FAC1 in the data.frame. Now i would like to scale all values according to the center of its experiment. So i can apply s