R 自举逻辑模型-某些子集不收敛
我想引导一个逻辑模型。整个数据集的模型收敛良好。但是,boot函数选择不再收敛的子集。我能做什么R 自举逻辑模型-某些子集不收敛,r,logistic-regression,glm,statistics-bootstrap,R,Logistic Regression,Glm,Statistics Bootstrap,我想引导一个逻辑模型。整个数据集的模型收敛良好。但是,boot函数选择不再收敛的子集。我能做什么 library(boot) set.seed(2) y <- c(rep(0,10),rep(1,10)) x <- c(rnorm(10,2,1),rnorm(10,6,1)) dat = data.frame(x, y) fit <- glm(y ~ x, quasibinomial(), data=dat) # Model with all data
library(boot)
set.seed(2)
y <- c(rep(0,10),rep(1,10))
x <- c(rnorm(10,2,1),rnorm(10,6,1))
dat = data.frame(x, y)
fit <- glm(y ~ x, quasibinomial(), data=dat) # Model with all data workes fine
bs <- function(data, indices) {
d <- data[indices,]
fitboot <- glm(y ~ x, family = quasibinomial(), data=d)
return(coef(fitboot))
}
results <- boot(data=dat, statistic=bs, R=10) # I get warnings
这似乎是由于选择了子集
有趣的是,这个子集是有效的:
fit <- glm(y ~ x, quasibinomial(), data=dat[1:13,])
fit也许可以尝试获取Wilcox包的源代码,该包有几个允许引导的函数。也许它还提供了其他有用的函数供将来使用。迭代重新加权最小二乘算法在默认的25次迭代中无法收敛。如果按如下所示增加迭代次数,警告将消失:
fitboot <- glm(y ~ x, family = quasibinomial(), data=d,
control=glm.control(maxit=50))
fitboot非常感谢您的回答,但我最终还是按照上面的建议做了。谢谢!!!这很有效!!!你帮了我很多。我花了一整天的时间在这上面,现在它起作用了。谢谢
fit <- glm(y ~ x, quasibinomial(), data=dat[1:14,])
fitboot <- glm(y ~ x, family = quasibinomial(), data=d,
control=glm.control(maxit=50))