Bootstrap для данных опроса, возвращающих одинаковые значения для каждой итерации - PullRequest
Новая
       51

Bootstrap для данных опроса, возвращающих одинаковые значения для каждой итерации

0 голосов
/ 09 мая 2020

Я пытаюсь запустить bootstrap на взвешенной по опросу логистической c регрессии, чтобы вычислить доверительные интервалы для оценок. Однако, когда я запускаю bootstrap, он возвращает одни и те же значения для оценок на каждой итерации, поэтому я предполагаю, что он не берет случайную выборку.

Ниже мой код 

logit.bootstrapsurvey <- function(data, indicies) {
  d<-svydesign(~SDMVPSU,data=data, strata=~SDMVSTRA, weights=~WTMEC2YR,nest=TRUE)
 d<-subset(d,!is.na(pos))
  fit <- svyglm(pos~factor(age)+factor(race)+factor(poverty)+factor(car)+factor(sex)+
                  factor(mar)+factor(eduref), design = d, family = "binomial")

  return(coef(fit))
}

set.seed(999)
logit.boot <-  boot(data=data_18, statistic=logit.bootstrapsurvey, R=1000)

Я предполагаю, что что-то в функции, которую я написал, заставляет функцию загрузки запускаться с исходными данными, а не с новым набором данных. Однако я не могу понять, как это происходит, и действительно ли это проблема.

...