我对R来说比较新,我试图使用Shaffer 2004创建并在R帮助文件系列中提供的更新版本的逻辑曝光链接功能代码来分析嵌套成功数据。 / p>
以下是我使用的代码:
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)
每次我尝试运行它时,都会收到以下错误:
&#34;错误:找不到有效的起始值:请指定一些&#34;
我尝试添加一个start = c(1,0)参数,但这没有效果。任何帮助将不胜感激!
答案 0 :(得分:1)
很晚才能复活这个问题,但我遇到了同样的问题并且刚刚找到答案。如果您的曝光日数据中有任何零值,则会出现上述有关无有效起始值的错误。将零更改为0.5或删除这些行可以解决问题: - )