| 查看: 894 | 回复: 1 | |||
刀光木虫 (著名写手)
铁虫
|
[求助]
SIMPLS算法
|
|
clear all clc load('data.mat'); %预先编写数据文件data.mat,并保存到当前工作路径下 X=data(:,1:5); y=data(:,6:8); % 1. 以SIMPLS算法进行 N=5; [XL,yl,XS,YS,beta,PCTVAR,MSE,stats] = plsregress(X,y,N); % 按照重组之后的成分对原始方差的解释力度 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') % 计算拟合优度，以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] % 权重分布图 figure(4) plot(1:N,stats.W,'o-'); legend({'c1','c2','c3','c4','c5','c6','c7'},'Location','NW') xlabel('Predictor'); ylabel('Weight'); 有哪个大神研究过这个程序，请问plsregress函数里面的各参数是指什么？% 权重分布图里面的c1-c7是指什么，% 权重分布图里面是指什么权重？感激不尽！！！ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有81人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
感觉是下周放榜了
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
w20073201
木虫 (正式写手)
小蜜蜂
- 应助: 5 (幼儿园)
- 金币: 2857.4
- 散金: 31
- 帖子: 372
- 在线: 355.6小时
- 虫号: 793821
- 注册: 2009-06-13
- 性别: MM
- 专业: 中药学其他科学问题
2楼2014-08-20 11:13:47
