99999久久久久久亚洲,欧美人与禽猛交狂配,高清日韩av在线影院,一个人在线高清免费观看,啦啦啦在线视频免费观看www

熱線電話：13121318867

登錄
2018-11-27 閱讀量: 817
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
0
關(guān)注作者
收藏
評論(0)

發(fā)表評論

暫無數(shù)據(jù)
推薦帖子