load cities
stdr=std(ratings);
sr=ratings./repmat(stdr,329,1);
[ Loadings, specificVar, T, stats]=factoran(sr,4,'rotate','varimax');
%subplot(121)
plot(Loadings (:,1),Loadings (:,2),'b.' )
text( Loadings (:,1),Loadings (:,2),num2str((1:10)'))
xlabel('czynnik 1')
ylabel('czynnik 2')
grid
xlim([-1 1])
ylim([-1 1])
axis square