Bootstrap在R

Question

我正在从名为nes2的数据中运行R中的非参数引导程序。这是我的代码：

co.store <- NULL
for(i in 1:500){ #we do this 500 times
  dat.bs <- nes2[sample(c(1:nrow(nes2)),nrow(nes2),replace=T),] #sampling that I did before. New data set
  logit.bs <- glm(formula=bs.model,data=dat.bs,family=binomial(link="logit")) #run logit with new data set
  co <- coefficients(logit.bs) #get coeff
  co.store <- rbind(co.store,co) #store coeff, rbind instead of c in order to have a stack of numbers. 

}

问题是我的新样本的系数完全相同。 我已经尝试了几件事，改变了整个模型中的公式（这是在公式而不是bs.model我试过放 nes2$turnout~nes2$age+nes2$agesqrd+hs.edu+college.edu+nes2$white+nes2$income这是整个模型）但它不会起作用

有什么想法吗？

谢谢！