parterror.m

来自「基于粗糙集的知识约简算法和采用回归进行的知识约简表示」· M 代码 · 共 15 行

function partE=partError(X,Y)
%计算划分X的恢复错误
X=rowzero0(X);
Y=Y/100;
LOG_Y=log(Y);
%B
B=inv(X'*X)*(X')*LOG_Y;
%RSE:restore error
Y_C=X*B;
[c,r]=size(Y);
partE=0;
for i=1:c
    partE=partE+((Y(i,1)-Y_C(i,1))/Y(i,1))^2;
end
partE=partE/c;