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

    function [L,OK,Stats] = fit(L,x,y)
%FIT general fit method for local and global models
%
% [L,OK,Stats] = fit(L,x,y)

%  Copyright 2006-2007 The MathWorks, Inc. 

if ~isa(x,'sweepset') || ~isa(y,'sweepset')
    % make a sweepset if necessary
    ss.varNames = [get(L,'symbols'),{varname(L),'logno'}];
    ss.data = [double(x),double(y),ones(size(y))];
    Data = struct2sweepset(sweepset,ss);
    % Sweepset has only one test
    Data = DefineSweepSet(Data,'logno',0.5);
    x = Data(:,1:nfactors(L));
    y = Data(:,end-1);
end
[L,Bhat,Wchat,OK]= fitmodel(L,x,y,allparameters(L),{[]});

if OK
    % set model
    L = update(L,Bhat);
    L = EvalDelG(L);
    L= pevinit(L,x,y);
    Stats  = FitSummary(L,x,y,Wchat);
else 
    Stats = [];
end