Displaying 20 results from an estimated 400 matches similar to: "na.action in stats::factanal() must be using formula interface and dataframe input to specify na.action?"
2013 Jan 02
1
Extracting factors from "factanal"
Dear R users
Happy New year to all for a start. Below is some data that I ran a factor
analysis on. Using $score prints the scores for each of the three factors.
However, I would like to access those factors as variable for new
computations. How do I do that? In SPSS we just call fact1_1, fact2_1 etc..
Thanks for your suggestions. V
============
v1 <- c(1,1,1,1,1,1,1,1,1,1,3,3,3,3,3,4,5,6)
v2
2006 Aug 11
1
- factanal scores correlated?
Hi,
I wonder why factor scores produced by factanal are correlated, and I'd
appreciate any hints from people that may help me to get a deeper
understanding why that's the case. By the way: I'm a psychologist used
to SPSS, so that question my sound a little silly to your ears.
Here's my minimal example:
***********************************************
v1 <-
2009 Mar 31
3
Factor Analysis Output from R and SAS
Dear Users,
I ran factor analysis using R and SAS. However, I had different outputs from
R and SAS.
Why they provide different outputs? Especially, the factor loadings are
different.
I did real dataset(n=264), however, I had an extremely different from R and
SAS.
Why this things happened? Which software is correct on?
Thanks in advance,
- TY
#R code with example data
# A little
2007 May 03
3
factanal AIC?
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
2007 May 03
3
factanal AIC?
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
2008 Sep 09
1
Addendum to wishlist bug report #10931 (factanal) (PR#12754)
--=-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 =
2009 Aug 11
3
loadings function (PR#13886)
Full_Name: Mike Ulrich
Version: 2.9
OS: Mac OSX
Submission from: (NULL) (69.169.178.34)
The help documentation for loadings() lists more then one parameter. The
function call only expects one parameter. The digits, cutoff, and sort
parameters are not used in the function.
## S3 method for class 'loadings':
print(x, digits = 3, cutoff = 0.1, sort = FALSE, ...)
## S3 method for class
2007 Jul 09
1
factanal frustration!
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
2004 Nov 30
1
Info
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) :
2008 Dec 01
1
factanal question
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
2009 Jan 30
1
Factor Analysis-factanal function
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
2006 Mar 15
3
Help on factanal.fit.mle
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]]
2009 Nov 03
1
Passing Command to Optim in factanal
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...
2007 May 13
1
factanal
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
2009 Apr 20
1
factanal error
Hi. I tried running the following code. I don't understand the error.
mydata <- read.table("C:/dataForR/radiology/wordFrequencies.csv",
header=TRUE, sep=",")
> dim(mydata)
[1] 982 924
mydataN<-mydata[,-923]
> dim(mydataN)
[1] 982 923
cor(mydataN)
factanal(mydataN, factors=3)
Error in solve.default(cv) :
Lapack routine dgesv: system is exactly singular
2007 Oct 06
2
factanal: error message
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
2012 Feb 07
1
na.action in stats::factanal()
> 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
2010 May 06
1
how to get components / factors in factanal / princomp not loadings
Dear all,
i wonder if there?s a command to obtain the actual values of a principal component or a factor (not as.factor, but factanal) .
test=princomp(USArrests, cor = TRUE)
summary(test)
just outputs, standard deviation, Prop of Variance and cumulative proportion of variance.
test$loadings offers yet another proportion of variance scheme. why is that?
Apart from that:
Is there a
2008 Jan 06
2
how to get residuals in factanal
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
2006 Feb 24
1
Extracting information from factanal()
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