Question

我试图从二元结果的对数二项式回归中获得 RR。有两个分类变量：治疗和组，两个数字变量：年龄和 BMI。但我收到一个错误 错误：找不到有效的起始值：请指定一些。请问如何解决此错误？

N <- 50
data.1 <- data.frame(Outcome=sample(c(0, 0, 1), N, rep=T), Age=runif(N, 8, 58),
                     BMI=rnorm(N, 25, 6), Group=rep(c(0, 1), length.out=N),
                     treatment=rep(c('1', '2', '3'), length.out=N))
data.1$Group <- as.factor(data.1$Group)
coefini <- exp(coef(glm(Outcome ~ Group + treatment + Age + BMI, data=data.1, 
                        family=binomial(link="logit"))))
fit2 <- glm(Outcome ~ Group + treatment + Age + BMI, data=data.1, 
            family=binomial(link="log"), start=coefini)

Answer 1

似乎是因为逻辑回归的系数不适用于对数二项式回归。用 coefini=coef(glm(Outcome~Group+treatment+Age+BMI, data=data.1,family =binomial(link = "log") )) 替换第三行就可以了。（删除 exp 并将链接更改为 log。）

具有多个类别预测变量和数字预测变量的二元结果的对数二项式回归

1 个答案: