cosyong 给定一组观测值,比如 x=c(5,4,9,6,21,17,11,20,7,10,21,15,13,16,8) 如何用boot包去计算var(s^2)的bootstrap估计,其中s^2是观测值x的样本方差?
winthander x=c(5,4,9,6,21,17,11,20,7,10,21,15,13,16,8) var_sample <- function(x,i) var(x) library(boot) var_boot <- boot(x,var_sample,R=500)
cosyong 非常感谢!不过你给出的是var(x)的bootstrap估计,并非var(s^2)的bootstrap估计 你写的程序结果为: ORDINARY NONPARAMETRIC BOOTSTRAP Call: boot(data = x, statistic = var_sample, R = 500) Bootstrap Statistics : original bias std. error t1* 34.31429 -2.550057 7.621101 但直接计算var(s^2)的bootstrap估计,结果应该为50左右,但我不知道如何用boot包去计算