21/02/20 10:04:00.43 bTTCnZQZ.net
>>148
simpleboot便利だな。
一々、sampleする関数を組まなくてもすむ。
library(simpleboot)
b.mean=one.boot(data=X,FUN=mean,R=5000,student = TRUE,M=5000)
boot.ci(b.mean)
X=c(40,45,46,47,49,52,52,56,69,72)
sim <- function(){
b=simpleboot::one.boot(X,mean,R=5000)
mean(b$t<55)
}
sim()
y=replicate(1e3,sim())
BEST::plotPost(y,col=4,showCurve = T)