[ndata,txt1,MixedData1]=xlsread(samplecode_path);
XMAT=ndata;
NALT=4;
n=floor(size(XMAT,1)/NALT);
a=size(XMAT,1)-n*NALT;
if a==0
XMAT(:,end+1)=repmat([1:NALT]',n,1);
else
b=repmat([1:NALT]',n,1);
B=[1:a]';
XMAT(:,end+1)=[b;B];
end
nn=size(XMAT,2);
for i=1:NALT;
XMAT=[XMAT (XMAT(:,nn)==i)];%生成虚拟变量。
end
m1nsa=ndata(:,5);
lm1=log(m1nsa);
gm1=diff(lm1);
ppi=ndata(:,2);
dlppi=diff(log(ppi));
D1=XMAT(:,8);
D2=XMAT(:,9);
D3=XMAT(:,10);
d1=D1(2:end,:);d2= D2(2:end,:);d3=D3(2:end,:);
%%%%%%%%%%%%%%%%%以下命令请帮助修改:
X=[d1 d2 d3];
T=size(d1,1);
xdata=cell(T,1);
for indx=1:T
xdata{indx}=repmat(X(indx,:),2,1);%indx*eye(3);%
end
Spec = vgxset('n',2,'nAR',1,'bsolve',true(3,1),'Constant',true) %'bsolve',true(3,1)
[EstSpec,EstStdErrors] = vgxvarx(Spec,Y,xdata);
vgxdisp(EstSpec,EstStdErrors)
|