RE: [R] R conversion

2004-09-09 Thread John Fox
Dear Mark,

See ?factanal for a factor-analysis function, in the standard stats package
(though factanal does ML factor analysis and not principal axes). Note that
help.search(factor analysis) turns this up.

With respect to your last question, it's not possible to know the source of
the difference without knowing what the difference is. A guess is that
specifying proportion=0.9 in SAS causes the program to use communality
estimates rather than 1's on the diagonal of the correlation matrix.

I hope this helps,
 John

 -Original Message-
 From: [EMAIL PROTECTED] 
 [mailto:[EMAIL PROTECTED] On Behalf Of Mark Strivens
 Sent: Thursday, September 09, 2004 7:26 PM
 To: [EMAIL PROTECTED]
 Subject: [R] R conversion
 
 I am a newcomer to R trying to convert a SAS program to R.
 Does anyone know if there is a functional equivalent of the 
 SAS 'Factor' procedure?
 
 For example in SAS:
 
 proc factor DATA=cor method=principal rotate=varimax 
 proportion=0.9 scree
 
 where 'cor' is a correlation matrix (as in the R 'cor' function)
 
 This should get you a list of eigen values as well as a 
 factor pattern matrix.
 
 Also why when I use the 'eigen' function in R does it seem to 
 give a subtly different answer to the eigen values generated 
 by the above program?
 
 Many thanks for any help
 
 Mark Strivens

__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] R conversion

2004-09-09 Thread Mark Strivens
Thanks John for the post - I had found factanal
However factanal it seems you have specify the number of factors to be
fitted
up front, whereas with the SAS procedure you don't - this is apparently
important
to the analysis!

Perhaps I could emulate this function by breaking it down using factor
extraction by principal component analysis and the doing the varimax
rotation second? If that makes
any sense?

__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html