R语言里两个样本的bootstrap码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

zahuy

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

但是我run出来的boot.mean.diff的t是NA...不知道是不是我function的codes（ 估计是t<-data$earn78[data$treated==1][i]和c<-data$earn78[data$treated==0][i]）写错了。



Bmeandiff <- function(data, i) {

  t<-data$earn78[data$treated==1][i]

  c<-data$earn78[data$treated==0][i]

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

# Set number of bootstrap replications
B=10000


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

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：实验组 是不是 控制组 对照组 每个人

没有呀

是两样本bootstrap吗

楼主还在吗，请教两样本bootstrap代码怎么写