Ошибка допустимых начальных значений для логистической регрессии GLM

Я относительно новичок в R и пытаюсь анализировать данные об успешном вложении, используя обновленную версию кода функции связи логистической экспозиции, созданную Shaffer 2004 и предоставленную в семействе файлов справки R.

Вот код, который я использую:

       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), но это не помогло. Любая помощь будет принята с благодарностью!