R语言中ATE遇到的问题

jynikita

1157

收藏 2019-07-18

悬赏 5 个论坛币未解决

求助大神们，我做了一个随机对照研究，其中有2组为干预组，1组为没有干预的对照组，最后结果是ITT没有意义，又因为其中一组的干预组依从率比较低，所以用了Per-Protocol，结果是阳性，但是2个审稿人均说PP不严谨，均提出需要用CACE/LATE来分析，由于没有接触过任何这方面的方法，我查了很多文献，发现我的研究是one-sided non-compliance，就是组间没有Crossover，在对照组的人不会接受任何干预，而干预1组和干预2组会出现没有接受干预的情况，但不会出现交叉，所以这个研究只会出现Complier和Never-taker，按照文献说的我应该是使用ATE(Average treamtent effect)、ATT（Average treament effect on the treated）/TOT（Treatment effect on the treated）进行分析。

以我的干预组与对照组的ATE分析为例
我把X（随机分组）编码为干预组=1，对照组=0，Y（结果）编码为阳性=1，阴性=0，Ti（实际接受干预）编码为接受=1，不接受=0
然后在R中下载ATE程序包，然后：
> library(ATE)
> #binary treatment

> Y<-table(TVSS$Y)

> treat<-table(TVSS$Ti)

> X<-table(TVSS$X)

> fit1<-ATE(Y,treat,X)
然后就提示
Error in ATE(Y, treat, X) :
The treatment levels must be labelled 0,1,2,...

ATT也是
> fit2<-ATE(Y,treat,X,ATT=TRUE)
Error in ATE(Y, treat, X, ATT = TRUE) :
The treatment levels must be labelled 0,1,2,...

试过很多次也是这样，但是介绍里明明是这样说的，Ti是二分类变量：

非常迷惘，请各位大神指正啊！！谢谢了！！！

扫码加我拉你入群

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

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

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群