R语言，bootstrap截断回归，面板数据，解决的朋友有额外现金报酬，急求，感谢

首先这是非面板，先建立返回回归向量的函数，再做回归，这个可以用
betas<-function(formula,data,indices){  
  d<-data[indices,]  
  fit<-lm(formula,data=d)  
  return(fit$coef)  
}  
results<-boot(data=paneldata,statistic=betas,R=500,formula=c~q+pf+lf)
print(results)
boot.ci(results,type="bca",index=2)

然后按照同样的方法扩展到面板为什么就不行了呢？
betas<-function(formula,data,indices){  
  d<-data[indices,]  
  fit<-plm(formula,data=d,model="within" ,index=c("id","year"))  
  return(fit$coef)  
}  
results<-boot(data=paneldata,statistic=betas,R=500,formula=c~q+pf+lf)  
print(results)  
boot.ci(results,type="bca",index=2)  

出的结果还是非面板的结果，应该是建立返回回归向量的公式有问题，报错如下
Error in pdim.default(index[[1]], index[[2]]) :
  duplicate couples (id-time) In addition: Warning messages:
1: In pdata.frame(data, index) :
  duplicate couples (id-time) in resulting pdata.frame
to find out which, use e.g. table(index(your_pdataframe), useNA = "ifany")
2: In is.pbalanced.default(index[[1]], index[[2]]) :
  duplicate couples (id-time)

3: In is.pbalanced.default(index[[1]], index[[2]]) :
  duplicate couples (id-time)

以及boot.ci(results,type="bca",index=2)  ，index=2是什么意思？到这里获得了置信区间，之后怎么办呢？怎么得到其他需要的参数？Z统计值，P值，显著性什么的？

可加847865855QQ私聊，完全解决问题现金100红包，求大神，在线等==急求

http://blog.csdn.net/yujunbeta/article/details/9255965

请问楼主解决了吗？我也需要做这个。

阿捷赫 发表于 2017-7-4 11:10
请问楼主解决了吗？我也需要做这个。

并没有，我已经放弃了R语言转投了MATLAB的怀抱，可能国内用R的太少了，每次遇到问题完全找不到答案太烦了。。。

lingdoujia 发表于 2017-9-30 13:44
并没有，我已经放弃了R语言转投了MATLAB的怀抱，可能国内用R的太少了，每次遇到问题完全找不到答案太烦了 ...

请问matlab咋弄的 我就想要一个能运行的程序 拜托 谢谢