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