目的描述
出于模型的需要，我们的团队选择做一次因子分析，通常这部分在队伍中是会有同学专门负责这块的，至于为什么笔者就不在这里多说了。
解决思路
在MATLAB中封装了有关因子分析的方法--factoran，读者可以通过help命令来查看如何调用这个方法。
需要读者注意的是，factoran方法并不会做数据规范化，所以读者需要自己来做这个操作。
代码展示
%数据单位化
TX_F_D=zscore(TX_D)
%5,求和
SUM_F_D=TX_F_D+AZ_F_D+CA_F_D+NM_F_D
%6，因子分析
[SUM_lambda,SUM_psi,SUM_T,SUM_stats,SUM_F]=factoran(SUM_F_D,6)
%求取贡献率
SUM_Contribute=Factor_Contribute(SUM_lambda,index)
%求各项因子的得分
ALL_F=Factor_F(TX_F_D,AZ_F_D,CA_F_D,NM_F_D,SUM_lambda,SUM_Contribute)
%画图
subplot(2,2,1)
plot(YEARS,ALL_F{1}(:,1),'r-',YEARS,ALL_F{2}(:,1),'g--',YEARS,ALL_F{3}(:,1),'b:',YEARS,ALL_F{4}(:,1))
xlabel('YEARS')
ylabel('F1')
legend('TX','AZ','CA','NM','Location','SouthEast')
%画单个的图
figure
plot(YEARS,ALL_F{1}(:,1),'r-',YEARS,ALL_F{2}(:,1),'g--',YEARS,ALL_F{3}(:,1),'b:',YEARS,ALL_F{4}(:,1))
xlabel('YEARS')
ylabel('F1')
legend('TX','AZ','CA','NM','Location','SouthEast')
figure
plot(YEARS,ALL_F{1}(:,2),'r-',YEARS,ALL_F{2}(:,2),'g--',YEARS,ALL_F{3}(:,2),'b:',YEARS,ALL_F{4}(:,3))
xlabel('YEARS')
ylabel('F2')
legend('TX','AZ','CA','NM','Location','SouthEast')
figure
plot(YEARS,ALL_F{1}(:,3),'r-',YEARS,ALL_F{2}(:,3),'g--',YEARS,ALL_F{3}(:,3),'b:',YEARS,ALL_F{4}(:,3))
xlabel('YEARS')
ylabel('F3')
legend('TX','AZ','CA','NM','Location','SouthEast')
figure
plot(YEARS,ALL_F{1}(:,4),'r-',YEARS,ALL_F{2}(:,4),'g--',YEARS,ALL_F{3}(:,4),'b:',YEARS,ALL_F{4}(:,4))
xlabel('YEARS')
ylabel('F')
legend('TX','AZ','CA','NM','Location','SouthEast')
结果展示
由于结果有多种多样的，直接给出MATLAB的工作空间，有兴趣的读者可以自行下载。
链接 密码：zffr