 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 844 | 回復(fù): 1 | ||
刀光木蟲(chóng) (著名寫(xiě)手)
鐵蟲(chóng)
|
[求助]
SIMPLS算法
|
|
clear all clc load('data.mat'); %預(yù)先編寫(xiě)數(shù)據(jù)文件data.mat,并保存到當(dāng)前工作路徑下 X=data(:,1:5); y=data(:,6:8); % 1. 以SIMPLS算法進(jìn)行 N=5; [XL,yl,XS,YS,beta,PCTVAR,MSE,stats] = plsregress(X,y,N); % 按照重組之后的成分對(duì)原始方差的解釋力度 figure(1) plot(1:N,cumsum(100*PCTVAR(2,  ),'-bo');xlabel('Number of PLS components'); ylabel('Percent Variance Explained in y'); % 殘差圖 figure(2) yfit = [ones(size(X,1),1) X]*beta; residuals = y-yfit; corrcoef(y,yfit); stem(residuals) xlabel('Observation'); ylabel('Residual'); % 擬合圖（y yfit） figure(3) plot(y,yfit,'o') % 計(jì)算擬合優(yōu)度，以R^2表示 y1=y(:,1); y2=y(:,2); y3=y(:,3); yfit1=yfit(:,1); yfit2=yfit(:,2); yfit3=yfit(:,3); TSS1 = sum((y1-mean(y1)).^2); RSS1 = sum((y1-yfit1).^2); Rsquared1 = 1-RSS1/TSS1 TSS2 = sum((y2-mean(y2)).^2); RSS2 = sum((y2-yfit2).^2); Rsquared2 = 1-RSS2/TSS2 TSS3 = sum((y3-mean(y3)).^2); RSS3 = sum((y3-yfit3).^2); Rsquared3 = 1-RSS3/TSS3 myRSS=[Rsquared1,Rsquared2,Rsquared3] % 權(quán)重分布圖 figure(4) plot(1:N,stats.W,'o-'); legend({'c1','c2','c3','c4','c5','c6','c7'},'Location','NW') xlabel('Predictor'); ylabel('Weight'); 有哪個(gè)大神研究過(guò)這個(gè)程序，請(qǐng)問(wèn)plsregress函數(shù)里面的各參數(shù)是指什么？% 權(quán)重分布圖里面的c1-c7是指什么，% 權(quán)重分布圖里面是指什么權(quán)重？感激不盡�。。� |
回復(fù)此樓» 猜你喜歡
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有7人回復(fù)
-
265求調(diào)劑
已經(jīng)有5人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有3人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
w20073201
木蟲(chóng) (正式寫(xiě)手)
小蜜蜂
- 應(yīng)助: 5 (幼兒園)
- 金幣: 2861.4
- 散金: 31
- 帖子: 372
- 在線: 354.4小時(shí)
- 蟲(chóng)號(hào): 793821
- 注冊(cè): 2009-06-13
- 性別: MM
- 專(zhuān)業(yè): 中藥學(xué)其他科學(xué)問(wèn)題
2樓2014-08-20 11:13:47
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 刀光 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
