r - GLM 逻辑回归的有效起始值错误

标签 r glm

我对 R 比较陌生，我正在尝试使用由 Shaffer 2004 创建并在 R 帮助文件 ?family 中提供的逻辑曝光链接函数代码的更新版本来分析嵌套成功数据。

这是我正在使用的代码:

       logexp <- function(days = 1) 
       {
       linkfun <- function(mu) qlogis(mu^(1/days)) 
       linkinv <- function(eta) plogis(eta)^days 
       mu.eta <- function(eta) days * plogis(eta)^(days-1) * binomial()$mu_eta  
       valideta <- function(eta) TRUE
       ink <- paste0("logexp(", days, ")")
       structure(list(linkfun = linkfun, linkinv = linkinv,
             mu.eta = mu.eta, valideta = valideta, name = link),
        class = "link-glm")
       }
       nestdata=read.table(file.choose(), header=TRUE)   ` 
       model=glm(survive~trtmnt,family=bionomial(logexp(days=nestdata$expos)),
       data=nestdata)

每次我尝试运行它时，都会收到以下错误:

“错误:找不到有效的起始值:请指定一些”

我尝试添加 start=c(1,0) 参数，但这没有效果。任何帮助将不胜感激!

最佳答案

很晚才重新提出这个问题，但我遇到了同样的问题并且刚刚找到了答案。如果您的暴露天数数据中有任何零值，您会收到上述关于没有有效起始值的错误。将零更改为 0.5 或删除这些行可以解决问题:-)

关于r - GLM 逻辑回归的有效起始值错误，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/22875502/