| 查看: 2711 | 回复: 0 | ||
刀光木虫 (著名写手)
铁虫
|
[求助]
PLS和SIMPLS算法
|
|
clear all clc load('data.mat'); %预先编写数据文件data.mat,并保存到当前工作路径下 X=data(:,1:5); y=data(:,6:8); % 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'); 请问[XL,yl,XS,YS,beta,PCTVAR,MSE,stats] = plsregress(X,y,N);中各参数是指什么？legend({'c1','c2','c3','c4','c5','c6','c7'},'Location','NW')权重分布图是指什么权重？如果要加一个预测模型要怎么加？ |
回复此楼» 猜你喜欢
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有7人回复
-
申请2026年博士
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
点击这里搜索更多相关资源