Bootstrap方法在R语言中的运用

最新推荐文章于 2023-08-08 23:51:06 发布

zeni910

最新推荐文章于 2023-08-08 23:51:06 发布

阅读量1.2w

点赞数 4

分类专栏：统计学 r 文章标签： bootstrap r语言统计学

本文链接：https://blog.csdn.net/yxfw00/article/details/111599254

版权

Bootstrap本意是拎着靴带让自己跳起来，在统计学中是一种重采样的方法，通常在样本量小的时候，可以从中不断再次抽样。

1. Bootstrap 简单应用

Bootstrap bias偏差 &variance 方差

library(bootstrap)
bslogc<-function(x,B){
   
  n<-length(x)
  thetastar<-replicate(B,{
   xstar<-sample(x,n,replace = T)
                       median(xstar)} )
  thetastar
}
 
x <- diabetes$logCpeptide   
hist(x, freq = FALSE, main = "Histogram of diabetes data")
median(x)
bscpep<-bslogc(x,1000)
hist(bscpep,freq = FALSE, breaks = 7,main = "Histogram of bootstrap diabetes medians")
var(bscpep)
bias<-mean(bscpep)-median(x)

Bootstrap correlation

samplecorr <- function(data, B) {
   
  n <- nrow(data)
  res <- numeric(B)
  for(i in 1:B) {
   
    ind <- sample(n, n, replace = TRUE)
    bs_data <-  data[ind, ]
    res[i] <-<