我想跑一个var的过程┌ x ┐ ┌ 1 ┐ ┌ 0.5 0.1 ┐┌ y(-1) ┐ ┌ 0 0 ┐┌ y(-2) ┐ ┌ e1 ┐
│ │=│ │+ │ │ │ │+│ ││ │+│ │
└ y ┘ └ 0 ┘ └ 0.4 0.5 ┘└ x(-1) ┘ └ 0.25 0 ┘└ x(-2) ┘ └ e2 ┘
┌ e1 ┐ ┌ 1 0 ┐
│ │~NID(0,I) ,I=│ │
└ e2 ┘ └ 0 1 ┘
所以写了以下程式
create u 105
smpl 1 2
genr x=3
genr y=3
smpl 1 105
rndseed 1200
rndseed 1000
genr e1=nrnd
genr e2=nrnd
genr [email=trend=@trend(0]trend=@trend(0[/email])
smpl 3 105
genr x=1+0.5*x(-1)+0.1*y(-1)+0*x(-2)+0*y(-2)+e1
genr y=0+0.4*x(-1)+0.5*y(-1)+0.25*x(-2)+0*y(-2)+e2
但是却无法动作 所以小弟我又换了以下程式
create u 105smpl 1 2
genr x=3genr y=3
smpl 1 105
rndseed 1300
genr e1=nrnd
genr e2=nrnd
smpl 3 105
vector(2) e
e.fill e1, e2,
matrix(2,2) I
I.fill(by=r) 1,0, 0,1,
matrix(2,1) a0
a0.fill(by=r) 1, 0,
vector(2) z
z.fill x, y,
vector(2) z1
z1.fill x(-1), y(-1),
vector(2) z2
z2.fill x(-2), y(-2),
matrix(2,2) a1
a1.fill(by=r) 0.5,0.1, 0.4,0.5,
matrix(2,2) a2
a2.fill(by=r) 0,0, 0.25,0,
z=a0+a1*z1+a2*z2+I*e
仍然不行 请问程式哪边出错?还请高人解答 谢谢