24小时热门版块排行榜

>论坛更新日志 (3701)
>文献求助 (346)
>虫友互识 (288)
>导师招生 (274)
>博后之家 (139)
>硕博家园 (138)
>休闲灌水 (136)
>论文投稿 (109)
>考博 (108)
>基金申请 (81)
>教师之家 (48)
>招聘信息布告栏 (47)
>考研 (40)
>公派出国 (39)
>绿色求助(高悬赏) (32)
>找工作 (29)

返回列表

刀光

木虫 (著名写手)

铁虫

应助: 1 (幼儿园)
金币: 21087.6
散金: 261
帖子: 2642
在线: 120.4小时
虫号: 2369003
注册: 2013-03-22
性别: GG
专业: 计算机软件

[求助] 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是指什么，% 权重分布图里面是指什么权重？感激不尽！！！

回复此楼