R里两个样本的bootstrap

热线电话：13121318867

啊啊啊啊啊吖

2018-11-27 阅读量: 726

R里两个样本的bootstrap

我的dataset有三列，第一列是treated，0表示是控制组，1表示是实验组，第二列是earn78，就是每个人的收入。然后我要看实验组和对照组的earn78值是不是不同。

但是我run出来的boot.mean.diff的t是NA...不知道是不是我function的codes（估计是t<-data

earn78[datatreated==1][i]和c<-data

earn78[datatreated==0][i]）写错了。

Bmeandiff <- function(data, i) {

t<-dataearn78[datatreated==1][i]

c<-data

earn78[datatreated==0][i]

return(mean(t)-mean(c))

}

# Set number of bootstrap replications

B=10000

boot.mean.diff <- boot(data=data, statistic=Bmeandiff, B)

0.0000

关注作者

发表评论

暂无数据

CDA考试动态

CDA报考指南

推荐帖子