function yy_after=correction(yy_before,sd,b,posi_prop,nega_prop,probdiffthreshold,K)
% yy_before:改正前预测值
% sd:预测不确定性
% yy_after:改正后预测值
% b:yy_before bar的中值
% posi_prop:正误差比例
% nega_prop:负误差比例

posi_prob=interp1(b,posi_prop,yy_before,'linear','extrap'); % 正改正概率，线性插值保证正反概率相加恒为一
nega_prob=interp1(b,nega_prop,yy_before,'linear','extrap'); % 负改正概率
ae=K*sd; % 改正值

yy_after=yy_before;
for i=1:length(yy_before)
    probdiff=abs(posi_prob(i)-nega_prob(i));
    if posi_prob(i)>nega_prob(i) && probdiff>probdiffthreshold
        yy_after(i)=yy_before(i)+ae(i);
    elseif posi_prob(i)<nega_prob(i) && probdiff>probdiffthreshold
        yy_after(i)=yy_before(i)-ae(i);
    end
end