matlab因子旋轉(zhuǎn)因子得分程序

matlab因子旋轉(zhuǎn)因子得分程序

ID:1486948

大?。?2.00 KB

頁數(shù):7頁

時間:2017-11-11

matlab因子旋轉(zhuǎn)因子得分程序_第1頁
matlab因子旋轉(zhuǎn)因子得分程序_第2頁
matlab因子旋轉(zhuǎn)因子得分程序_第3頁
matlab因子旋轉(zhuǎn)因子得分程序_第4頁
matlab因子旋轉(zhuǎn)因子得分程序_第5頁
資源描述:

《matlab因子旋轉(zhuǎn)因子得分程序》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫

1、因子旋轉(zhuǎn)因子得分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)=d

2、ata(:,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

3、','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''X2

4、''X3''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])boxon

5、holdoff%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%plotsecondfactoragainstt

6、hirdsubplot(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%MaximumLikelihoodFactorAnalysisaftervar

7、imaxrotationlambda=rotatefactors(ld,'Method','varimax');%rotatesthefactorloadingsmatrix%andestimatesfactorloadingsaftervarimaxvl=[lambda(:,1),lambda(:,2),lambda(:,3)];com=diag(vl*vl');%communalitiesarecalculatedpsi=diag(dat)-diag(vl*vl');%specificvariancesarecalculatedtbl=[vl,

8、com,psi]figure(2)subplot(2,2,1)holdon%plotfirstfactoragainsts

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動畫的文件,查看預(yù)覽時可能會顯示錯亂或異常,文件下載后無此問題,請放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫負責(zé)整理代發(fā)布。如果您對本文檔版權(quán)有爭議請及時聯(lián)系客服。
3. 下載前請仔細閱讀文檔內(nèi)容,確認文檔內(nèi)容符合您的需求后進行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時可能由于網(wǎng)絡(luò)波動等原因無法下載或下載錯誤,付費完成后未能成功下載的用戶請聯(lián)系客服處理。