x=c(5,4,9,6,21,17,11,20,7,10,21,15,13,16,8)
n=length(x)
boot=10000
s2=NA*seq(1:boot)
for(i in 1:boot)
{
xsamp=sample(x,n,replace = T)
s2[i]=sum((xsamp-mean(xsamp))^2)/(n-1)
}
var(s2)
#[1] 57.49942

2# epoh

s2=NA*seq(1:boot)     epoh老师，您好！这一行是什么意思？NA并没有赋值，NA乘以seq序列是做什么用的？
非常感谢！

不好意思，最近出去玩了，才致许久没有回帖，对不起哈。毫无疑问，最佳答案是epoh大师。看来bootstrap方法我还没有领会精神，得要加油了。谨再谢！

但是epoh大师也是直接手工计算的，并没有用boot package来计算

后来我在cos统计之都上找到了答案，现在贴出来供大家分享和讨论：
> library(boot)#当然要先下载安装boot包
> set.seed(200)
> x=c(5,4,9,6,21,17,11,20,7,10,21,15,13,16,8)
> x.boot <- boot(x,function(x,i)var(x[i]),R=999)
> x.boot

ORDINARY NONPARAMETRIC BOOTSTRAP


Call:
boot(data = x, statistic = function(x, i) var(x[i]), R = 999)


Bootstrap Statistics :
    original    bias    std. error
t1* 34.31429 -2.203737     7.32227
大家看到了std.error 下面对应的7.32227吗
var(x)的bootstrap估计就是7.32227的平方 53.61564