[求助]我编的bootstrap程序，高手帮我看看有什么问题！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

建立var（2）模型并在此基础上做bootstrap因果关系检验，程序写出来拿到R上运行，还是有些问题。
我把程序发给高手们看看，问题出在什么地方。
另外，如果已经计算出了残差，是否可以直接赋值为a，运行r以下的程序。
诚盼指点！
x1=scan(1)  
x2=scan(2)             
x=ts(cbind(x1,x2)
l=var(x,p=2,type="const")   显示有问题
l1=coef(l)$x1
l2=coef(l)$x2
a=residusals(l)

m1=matrix(0,26,1000)
for(i,in 1:1000){m1[,i]=sample(a[,1],26)}   a[,1]显示有问题
m2=matrix(0,26,1000)
for(i in 1:1000){m2[,i]=sample(a[,2],26)}
x11=matrix(0,26,1000)
x21=matrix(0,26,1000)
for(i in 1:1000){for(j in 1:26){x11[j+2,i]=l1[5]+l1[1]x1[j+1,i]+l1[3]x1[j,i]+l1[2]x2[j+1,i]+l1[4]x2[j,i]+m1[j,i]}}
for(i in 1:1000){for(j in 1:26){x22[j+2,i]=l2[5]+l2[1]x1[j+1,i]+l2[3]x1[j,i]+l2[2]x2[j+1,i]+l2[4]x2[j,i]+m2[j,i]}}

得到X11 X21

[此贴子已经被作者于2009-3-24 8:35:53编辑过]

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Bootstrap Bootstra boots Trap boot 程序 高手 Bootstrap

无解，鉴定完毕！！！

哪位高手帮助解决了，请吃大餐！：）

做VAR(2)是用var函数吗？你用的是哪个package?stats有个var函数是求方差的，如果误用了var，那后面的结果肯定都是错了！另外你的scan()使用好像也是错的！

谢谢，R中的vars package 可以做VAR模型。至于scan调用数据，发上来时我简化了，这一点不存在问题。还请帮忙看看下面的程序！

请求高手帮忙！