restore_e.m

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

function [RSE]=restore_E(filename,attrnum)
%求restore error
%X
X_R=txttomt01(filename,attrnum);
[X,XATTR]=rowzero(X_R);
[c1,r1]=size(X);
ITEMNUM=c1;
%Y
Y=txttomtY(filename);
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);
RSE_ALL=0;
for i=1:c
    RSE_ALL=RSE_ALL+((Y(i,1)-Y_C(i,1))/Y(i,1))^2
end
RSE=RSE_ALL/ITEMNUM;