fitted.factanal(object) residuals.factanal(object)
residuals.factanal returns a matrix which is the difference between the actual correlation matrix and the fitted matrix returned by fitted.factanal.
prim9.fa3 <- factanal(prim9, 3)fitted(prim9.fa3) residuals(prim9.fa3)