悬赏 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命令是什么意思。求高人指点,帮忙解读该程序