我的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)
三个资料Q群下载不了也转发不了,先放这里Fine_tuning.zipLangChain.zipdata_clear.rar