parazhu 发表于 2013-2-19 08:35 
不好意思,我看错了。我猜测你的datsam是你原始的数据,那么datsam是一个dataframe,而indice是抽样之后的位 ...
谢谢parazhu的回复!我按照你的提示改动了,并且试了几次后发现大括号里的几个选项要分行写才能运行,我原来是同一行输入的。如下的程序运行就可以了:
>-library(boot)
>-s1 <- function(datsam, indices){
+d<-datsam[indices,]
+fitlr <- glm(case ~ packyr * radon+age, family=binomial, data=d)
+s <- (exp(fitlr$coef[2]+fitlr$coef[3]+fitlr$coef[5])-1)/(exp(fitlr$coef[2])-1+exp(fitlr$coef[3])-1)}