начальная загрузка вмененного вывода cem - PullRequest
Новая
       5

начальная загрузка вмененного вывода cem

2 голосов
/ 20 сентября 2010

Я делаю грубое точное сопоставление с вмененным набором данных. Для проверки я делаю начальную загрузку результатов из 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)

В основном все работает гладко, за исключением того, что я не могу получить загрузочный CI bca.

...