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);