问题是找出一个存在方差最小的门限值. 我写出来代码,不过中间有错误," ??? Error using ==> le
Matrix dimensions must agree." 不知道如何改..现在就卡在这里..希望有高手能帮我看下,点拨一下,谢谢了!!
数据在附件里
load dat
n=96;
q=sort(q);
GA= q(floor(n*0.1): floor(n*0.9));
S=zeros(1,length(GA));
for i = 1:length(GA);
gama= GA(i);
XL=X(q<=gama,: );
YL=y(q<=gama);
B1=(XL'*XL)\(XL*YL);
eL=YL-XL*B1;
XR=X(q>gama,: );
YR=y(q>gama);
B2=(XR'*XR)\(XR*YR);
eR=YR-XR*B2;
S(i)=sum(eL)+sum(eR);
end
[C,I]=min(S);
Gamamin=GA(I);