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

    function mdev=InitStore(mdev,ind);
%INITSTORE initialises twostage models

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



if nargin==1
	ind= BMIndex(mdev);
end
if ind==0
	return
end

X= getdata(mdev,'FIT');
XG = X{end};

TS= mdev.TwoStage{ind};

Xcode= gcode(TS,double(XG));

L= model(mdev);
RF1= RFstart(L);
DatumType= get(L,'DatumType');
switch DatumType
    case {1,2}
        prf= children(mdev);
        Ydatum= double(prf(1).getdata('Y'));
    case 3
        pdatum= datumlink(mdev);
        Ydatum= double(pdatum.getdata('Y'));
    otherwise
        Ydatum= 0;
end

selrf= mdev.ResponseFeatures(ind,:);


Yrf= children(mdev,selrf+RF1,@getdata,'Y');
Yrf= [Yrf{:}];

TS= InitStore(TS,Xcode,Yrf,Ydatum);
ri= var(TS);
if isempty(ri);
    mdev.TwoStage{ind}= TS;
	TS= pevinit(mdev,ind);
end
mdev.TwoStage{ind}= TS;
if ind==2 && ismle(TS)
    mdev.MLE.Model= TS;
end

pointer(mdev);