www.gusucode.com > mbctools 工具箱 matlab 源码程序 > mbctools/@modeldev/diagnosticStats.m
function [data,factors,standardPlotStr,olIndex]= diagnosticStats(mdev) %DIAGNOSTICSTATS % Copyright 2000-2015 The MathWorks, Inc. and Ford Global Technologies, Inc. m= model(mdev); [X,Y,DataOK]= FitData(mdev); Y(~DataOK)= NaN; if status(mdev) % get new stats [data,factors,standardPlotStr]= diagnosticStats(m,X,Y); if mdev.ModelStage==1 % add monitor variables to scatter plot TP= mdevtestplan(mdev); mvars= getMonitor(TP); if ~isempty(mvars) && ~isempty(mvars.values) MDATA= getdata(TP,'ALLDATA'); MDATA= MDATA(:,mvars.values); if size(X,1)==size(X,3) MDATA= MDATA(DataOK,:); data = [data double(MDATA)]; factors= [factors(:); mvars.values(:)]'; end end end % get outlier indices olIndex= outliers(m,data,factors); if ~isRespFeat(mdev) % add validation residual plots standardPlotStr= [standardPlotStr 'Validation residuals']; end else data = [double(X(DataOK,:)),double(Y(DataOK))]; factors = [get(X,'name');get(Y,'name')]; olIndex = []; standardPlotStr= {}; end