我正在尝试在调查加权逻辑回归上运行引导程序,以计算估计值的置信区间。但是,当我运行引导程序时,它在每次迭代时都会返回相同的估算值,因此我假设它不是在抽取随机样本。
下面是我的代码
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)
我的假设是,我编写的函数中的某些内容导致引导函数在原始数据而非新数据集上运行。但是我无法弄清楚它是如何做到的,或者这是否确实是问题所在。