similar to: Ignoring Errors in Simulations

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

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

Hi Can anybody please suggest me about the documentation of "factanal.fit.mle()" (Not factanal()------ searching factanal.fit.mle() in R always leads to factanal()). Is there any function for doing principal component factor analysis in R. Regards Souvik Bandyopadhyay JRF, Dept Of Statistics Calcutta University [[alternative HTML version deleted]]

Thanks to Brian Ripley for the time and effort put into developing the factanal package for R. I have found it very useful. However, in trying to replicate some results from previous research, I have run into the need for a couple of extensions and was wondering if they might find their way into future implementations. 1. Though maximum likelihood estimates might be more rigorous and

Hi there, Does anyone know how to explicitly refer to the p-value of thet test that the chosen number of factors is significant in a factor analysis. It's not in the list of values for the factanal command output yet it is printed out with the results. Thanks in advance. Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James's Buildings 79 Oxford Street

Hello! I have a question about the factanal function. This function returns at the end test statistics like this: Test of the hypothesis that 4 factors are sufficient. The chi square statistic is 4.63 on 2 degrees of freedom. The p-value is 0.0988 Is it possible to get the chi square statistic and the p-value as variables, not the text on the screen? An object of class "factanal"

Dear R-help, Using R 1.2.1 on Windows98 to run a factor analysis on a 64x150 matrix of data generated from a simulation model, factanal() reported that it failed to find a solution. Looking at the factanal code, I see the immediate condition that triggered the result: if (best == Inf) stop("Unable to optimize from these starting value(s)") So I am sure factanal() is giving

# 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

Dear expe-R-ts, I try some test data for a factorAnalysis (resp. pca) in the sense of Prof. Ripley's MASS ? 11.1, p. 330 ff., just to prepare myself for an analysis of my own empirical data using R (instead of SPSS). 1. the data. ## The test data is (from the book of Backhaus et al.: Multivariate ## Analysemethoden. Springer 2000 [9th ed.], p. 300 ff):

I have put a preliminary version of factanal in a package at http://www.stat.ox.ac.uk/pub/R/fa_0.1-1.tar.gz and would appreciate feedback. In particular 1) What else you might want. Please give precise references with algorithmic details, as that has been a major problem for me. 2) Some convincing examples of factor analysis. Most (literally) of the ones I have looked at are local maxima, not

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

Dear R users: I'm wondering if it's possible to get the residual correlation matrix when using factanal. Since factanal assumes that the errors are normally distributed and independent (provided the factor model fits the data) this would be useful. Of course you would need to submit the data to the function to get the residuals (not just their correlation matrix), but it should be possible

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

Dear list members, I apologize for putting this (probably) very basic question on the mailing list. I have scanned through the R website (using search) but did not found an answer. (code included below) A factor matrix is simply extracted (which can then subsequently be exported using write.table) by FACT$loadings[1:6,]. I would also like to specifically extract and export

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

Hi. It seems that nearly every time I try to use factanal I get the following response: >faa2db1<-factanal(mretdb1,factors=2,method="mle",control=list(nstart=25)) Error in factanal(mretdb1, factors = 2, method = "mle", control = list(nstart = 25)) : unable to optimize from these starting value(s) > In the case cited above, mretdb1 is synthetic data created

I am sure there is a simple way to do the following, but i haven't been able to find it. I am hoping a merciful soul on R-help could point me in the right direction. I am doing a factor analysis on survey data with missing values. to do this, I run: FA1<-factanal(na.omit(DATA), factors = X, rotation = 'oblimin', scores = 'regression') Now that I have my factors and

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

> Does factanal() force the user to use the formula interface if they > wish to specify an na.action? Short answer: yes. Long answer: The handling of na.action is a built in part of the formula processing functions, so it's automatic when dealing with a formula. There are also downstream effects on predict() and resid() that are worked out for the formula case, but aren't clear