

I have a dataset, d1, on which I'd like to perform a factor analysis. 
Because there are some missing values that I don't want to lose or 
interpolate, and because factanal is less tolerant of these in a dataset 
than it is in a correlation matrix (it requires na.action=na.omit), I 
calculated a correlation matrix as follows:

r.matrix<-cor(d1, use="pairwise.complete.obs")
There are 44 cases in the dataset

However, when I type the command for factor analysis

Factor2<-factanal(d1,2,covmat=r.matrix,n.obs=44,scores="regression")

I get the error message

Error in factanal(d1, 2, covmat = r.matrix, n.obs = 44, scores = 
"regression") : requested scores without an 'x' matrix

What am I doing wrong?

Thanks,  Robert Seyfarth


	[[alternative HTML version deleted]]

