Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/73.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181

Warning: file_get_contents(/data/phpspider/zhask/data//catemap/3/reactjs/21.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 自举逻辑模型-某些子集不收敛_R_Logistic Regression_Glm_Statistics Bootstrap - Fatal编程技术网

R 自举逻辑模型-某些子集不收敛

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

我想引导一个逻辑模型。整个数据集的模型收敛良好。但是,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 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))