- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 0 个
- 通用积分
- 0
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 189 点
- 帖子
- 15
- 精华
- 0
- 在线时间
- 15 小时
- 注册时间
- 2009-4-19
- 最后登录
- 2014-5-5
高中生
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
 - 0 个
- 通用积分
- 0
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 189 点
- 帖子
- 15
- 精华
- 0
- 在线时间
- 15 小时
- 注册时间
- 2009-4-19
- 最后登录
- 2014-5-5
|
5论坛币
|
%检验交叉有效性停止,取m个成分
% F0=r1t1+...+rmtm,且th=Eh-1*wh=E0*wh^*
% F0=r1E0w1^*+...+rmE0wm^* = E0(sum(rh*wh^*))
% y(估计)=a0+a1*x1+ … +ap*xp
y=F0;%测得的真实数据
yh1=r1*t1;%求出y在取t1成分回归建模后的拟合值(所有样本)
yh2=r1*t1+r2*t2;
yh3=r1*t1+r2*t2+r3*t3;
yh4=r1*t1+r2*t2+r3*t3+r4*t4;%此处如果是多个因变量,应用r'代替r
SS0=sum(y.^2);
press1=getpress(E0,F0,1)
Q1=1-press1/SS0;
SS1=sum((y-yh1).^2);
press2=getpress(E0,F0,2)
Q2=1-press2/SS1;
SS2=sum((y-yh2).^2);
press3=getpress(E0,F0,3)
Q3=1-press3/SS2;
SS3=sum((y-yh3).^2);
press4=getpress(E0,F0,4)
Q4=1-press4/SS3;
SS4=sum((y-yh4).^2);
press5=getpress(E0,F0,5)
Q5=1-press5/SS4;
RMSEC=sqrt(press5/11)
这是pls1程序中用交叉有效性确定提取成分个数的一段程序。不清楚getpress命令是什么意思。求高人指点,帮忙解读该程序 |
|