www.gusucode.com > mbctools 工具箱 matlab 源码程序 > mbctools/@mdev_local/localstats.m

    function [OK,Stats,InSig,labs,head]=localstats(mdev,SNo,Type,L,XL,Y)
%LOCALSTATS local statistics
%
% [Stats,InSig,labs,head]=localstats(mdev,SNo,Type,L,X,Y);

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


if nargin<=3
    % get local model and data
    [L,OK]=LocalModel(mdev,SNo);
    [XL,Y]= getdata(mdev);
    XL= XL(:,:,SNo);
    Y= Y(:,:,SNo);
else
    OK= mdev.FitOK(SNo);
end


if OK
    % summary stats
    TS= BestModel(mdev);

    % get global data
    [Xtp,Ytp]= getdata(mdevtestplan(mdev),'FIT');
    X= {XL Ytp(:,get(Xtp{2},'Name'),SNo)};

    ValRMSE= LocalValidationRMSE(mdev,SNo,L);
    Wc= mdev.GLSWeights{SNo};
    [labs,head,Stats,InSig]= localsummary(L,Type,TS,X,Y,Wc,ValRMSE);
else
    Stats=[];
    InSig=[];
    labs=[];
    head=[];
end