Как мне написать свою собственную функцию, которая создаст образец bootstrap - PullRequest
Новая
       8

Как мне написать свою собственную функцию, которая создаст образец bootstrap

1 голос
/ 26 апреля 2020

Я хочу написать функцию, которая дает мне доверительный интервал bootstrap для любого вектора (x), который будет выполняться l oop (n) раз. Я хочу оценить средний индекс массы тела для человека с сердечно-сосудистыми заболеваниями, используя имеющиеся у меня данные выборки (cvdonly $ cbmi).

Это то, что у меня есть, но оно не работает, потому что производимый им КИ очень отличается от среднего значения моего набора данных выборки: 

bs<-function(x,n){
  m_vec<-rep(NaN, n)
  for (i in 1:n){
    cbminorm<-rnorm(x)
    m_vec[i]<-mean(cbminorm)
  }
  ll<-round(quantile(m_vec, 0.025, na.rm = TRUE),2)
  ul<-round(quantile(m_vec, 0.975, na.rm = TRUE),2)
  paste("95% bootstrap inverval for the mean:", ll, "to", ul)
}

bs(cvdonly$cbmi, 100)
...