matlab因子旋转因子得分程序

因子旋转因子得分clearall;closeall;clc;%loaddatadata=xlsread('stockdata');%transformdataxt=data;xt(:,[1,3,10,11])=data(:,[1,3,10,11])./(10.'9);xt(:,[2,4])=data(:,[2,4])./(10."10);xt(:,5)=data(:,5)./(10."5);xt(:,[6,12,13])=data(:,[6,12,13])./(10.'11);xt(:,7)=data(:,7)./(10."2);xt(:,8)=data(:,8).*10;xt(:,9)=data(:,9)./10;fori=1:13zz(:,i)=(xt(:,i)-mean(xt(:,i)))/(sqrt(var(xt(:,i))));end%correlationmatrixdat=corr(zz);%MaximumLikelihoodFactorAnalysiswithoutvarimaxrotation%factanalperformsmaximum-likelihoodfactoranalysis[lambda,psi,T,stats,F]=factoran(zz,3,'xtype','data','rotate','none');%estimatedfactorloadings%theestimatedfactorloadingsmatrixld=lambda;com=diag(ld*ld');%communalitiesarecalculatedpsi=diag(dat)-diag(ld*ld');%specificvariancesarecalculatedtbl=[lambda(:,1),lambda(:,2),lambda(:,3),com,psi]figure(1)subplot(2,2,1)s=['X1''X5''X6''X7''X8''X9''X10''X11''X12''X13''X14'];%plotfirstfactoragainstsecondholdonplot(lambda(:,1),lambda(:,2),'w')ylim([-0.7,0.7])title('Factors21-theta21')ylabel('y')xlabel('x')fori=1:13text(lambda(i,1),lambda(i,2),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%plotfirstfactoragainstthirdsubplot(2,2,3)holdonplot(lambda(:,1),lambda(:,3),'w')xlabel('x')ylabel('y')title('Factors31-theta31')ylim([-0.5,0.5])fori=1:13text(lambda(i,1),lambda(i,3),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%plotsecondfactoragainstthirdsubplot(2,2,2)holdonplot(lambda(:,2),lambda(:,3),'w')xlabel('x')ylabel('y')title('Factors32-theta32')xlim([-0.7,0.7])ylim([-0.5,0.5])fori=1:13text(lambda(i,2),lambda(i,3),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%MaximumLikelihoodFactorAnalysisaftervarimaxrotationlambda=rotatefactors(ld,'Method','varimax');%rotatesthefactorloadingsmatrix%andestimatesfactorloadingsaftervarimaxvl=[lambda(:,1),lambda(:,2),lambda(:,3)];com=diag(vl*vl');%communalitiesarecalculatedpsi=diag(dat)-diag(vl*vl');%specificvariancesarecalculatedtbl=[vl,com,psi]figure(2)subplot(2,2,1)holdon%plotfirstfactoragainstsecondplot(lambda(:,1),lambda(:,2),'w')xlabel('x')ylabel('y')title('Factors21-theta21')fori=1:13text(lambda(i,1),lambda(i,2),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%plotfirstfactoragainstthirdsubplot(2,2,3)plot(lambda(:,1),lambda(:,3),'w')xlabel('x')ylabel('y')title('Factors31-theta31')fori=1:13text(lambda(i,1),lambda(i,3),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%plotsecondfactoragainstthirdsubplot(2,2,2)holdonplot(lambda(:,2),lambda(:,3),'w')xlabel('x')ylabel('y')title('Factors32-theta32')fori=1:13text(lambda(i,2),lambda(i,3),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%PrincipalComponentMethodaftervarimaxrotation%spectraldecomposition[eigvec,eigval]=eigs(dat);E=ones(size(eigvec(:,1:3)))*eigval(1:3,1:3);%theestimatedfactorloadingsmatrixQ=sqrt(E).*eigvec(:,1:3);lambda=rotatefactors(Q,'Method','varimax');%rotatesthefactorloadingsmatrix%andestimatesfactorloadingsaftervarimaxld=[lambda(:,1),lambda(:,2),lambda(:,3)];com=diag(ld*ld');%communalitiesarecalculatedpsi=diag(dat)-diag(ld*ld');%specificvariancesarecalculatedtbl=[ld,com,psi]figure(3)%plotfirstfactoragainstsecondsubplot(2,2,1)holdonplot(lambda(:,1),lambda(:,2),'w')xlabel('x')ylabel('y')title('Factors21-theta21')fori=1:13text(lambda(i,1),lambda(i,2),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%plotfirstfactoragainstthirdsubplot(2,2,3)plot(lambda(:,1),lambda(:,3),'w')xlabel('x')ylabel('y')title('Factors31-theta31')fori=1:13text(lambda(i,1),lambda(i,3),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%plotsecondfactoragainstthirdsubplot(2,2,2)holdonplot(lambda(:,2),lambda(:,3),'w')xlabel('x')ylabel('y')title('Factors32-theta32')fori=1:13text(lambda(i,2),lambda(i,3),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%PrincipalFactorMethodaftervarimaxrotation%inverseofthecorrelationmatrixdatf=inv(dat);%preliminaryestimateofpsipsiini=diag(diag(1./f));psi=psiini;[eigvec,eigval]=eigs(dat-psi);EE=ones(size(eigvec(:,1:3)))*eigval(1:3,1:3);QQ=sqrt(EE).*eigvec(:,1:3);psiold=psi;psi=diag(1-sum((QQ.*QQ)'));z=psi-psiold;convergence=z;lambda=rotatefactors(QQ,'Method','varimax');%rotatesthefactorloadingsmatrix%andestimatesfactorloadingsaftervarimaxld=[lambda(:,1),lambda(:,2),lambda(:,3)];com=diag(ld*ld');%communalitiesarecalculatedpsi=diag(dat)-diag(ld*ld');%specificvariancesarecalculatedtbl=[ld,com,psi]figure(4)%plotfirstfactoragainstsecondsubplot(2,2,1)holdonplot(lambda(:,1),lambda(:,2),'w')xlabel('x')ylabel('y')title('Factors21-theta21')fori=1:13text(lambda(i,1),lambda(i,2),s(i,1:3));end;line([-1,1],[0,0])line([0,0],[-1,1])boxonholdoff%plotfirstfactoragainstthirdsubplot(2,2,3)plot(lambda(:,1),lambda(:,3),'w')xlabel('x')ylabel('y')title('Factors31-theta31')fori=1:13text(lambda(i,1),lambda(i,3),s(i,1:3));end;line([-1,1],[0,0])line([0