www.gusucode.com > mbcmodels 工具箱 matlab 源码程序 > mbcmodels/@xreglinear/mse.m

    function s2=mse(m);
%MSE

%  Copyright 2000-2004 The MathWorks, Inc. and Ford Global Technologies, Inc.



y= m.Store.y;
Q= m.Store.Q;

y= ytrans(m,y(isfinite(y)));

yq= Q'*y;

sse= y'*y-yq'*yq;

df= length(y)-sum(~m.TermsOut);
s2 = sse/df;