Als toevoeging: ik gebruik zelf altijd de functie fa uit het package psych. Dit is een superhandig package voor o.a. psychometrie.
Package userfriendlyscience combineert een aantal van die functies in de functie fullFact (zie ?fullFact).
Wat ook heel handig kan zijn is de functie factorLoadingDiamondCIplot: die visualiseert de betrouwbaarheidsintervallen voor de factorladingen, bijvoorbeeld:

Dit is handig omdat besluiten over factorladingen op basis van puntschattingen (net als alle besluiten op basis van puntschattingen) geen goed idee zijn. Die puntschattingen verschillen namelijk van steekproef tot steekproef, en kunnen de volgende keer volslagen anders zijn. Als je steekproef groot genoeg is, kunnen ze van steekproef tot steekproef natuurlijk minder verschillen.
Hoeveel kunnen die puntschattingen van steekproef tot steekproef verschillen? Daarover heeft het betrouwbaarheidsinterval meer informatie.
De puntschatting is alleen geldig in je steekproef. Maar, je bent meestal niet geinteresseerd in je steekproef: die is alleen maar een middel om iets te zeggen over je populatie. Je betrouwbaarheidsinterval kun je gebruiken om een idee te krijgen van van wat plausibele waarden in de populatie zijn. Hoe groter de betrouwbaarheid, hoe plausibeler het is dat de populatiewaarde (bijvoorbeeld de populatie factorlading) in het betrouwbaarheidsinterval uit een gegeven steekproef ligt.
Dus: als je bijvoorbeeld 95% betrouwbaarheidsintervallen gebruikt (de 'default'), dan geldt dat als de betrouwbaarheidsintervallen van factorladingen overlappen, het niet vanzelfsprekend is om te concluderen dat die factorladingen in de populatie verschillen. Zoals je ziet geldt dat hier voor de meeste factorladingen: de evidentie dat die factorladingen verschillen is dus niet sterk.
Deze dataset ('Thurstone.33', een standaard dataset in R, zie het voorbeeld van deze functie) betreft dan ook maar 9 regels.
Dus, ik zou vooral de combinatie van fa en factorLoadingDiamondCIplot gebruiken, zoals in dit voorbeeld:
factorLoadingDiamondCIplot(fa(Thurstone.33, 2,
n.iter=100, n.obs=100),
alpha = .5);