r - 自举估算的cem输出

Question

我正在对估算的数据集进行粗化的精确匹配。为了验证我对 cem 的结果进行了引导。但是我无法让boot.ci输出正常工作。我收到错误消息：

Error in bca.ci(boot.out, conf, index[1L], L = L, t = t.o, t0 = t0.o,  : 
  estimated adjustment 'a' is NA

物体在那里，但由于某种原因，我无法得到偏差校正的结果。

library(boot)
library(cem)
library(Amelia)

am.bo<-function(data,ind) {
  data<-data[ind,]
  a9<-amelia(data,ords=c("gck","gc.m","gc.e","gc.v"),id=c("as2","d30"))
  os<-a9$imputations[1:5]
  os2<-cem(treatment="as2",datalist=os,data=data,drop="bc")
  os3<-att(os2,gck~bc,data=os)
  print(os3,T)->g
  intercept<-g[1,1]
  bac<-g[1,2]
  vec<-c(intercept,bc)
  return(vec)
}

dlply(a4,.(as),function(x) boot(x,am.bo,R=30)) -> bootas_gcadm
lapply(bootas_gcadm,boot.ci)

基本上一切都运行顺利，除了我无法获得bca引导 CI。