Question

我正在尝试按年龄对我的 glm 模型进行分层。

我的主要问题是：我可以创建一个名为 BYDECADE 的新列并按如下方式运行代码，而不是在 glm 中执行此子集设置吗？

 fit <- glm(COMPLICATION~BYDECADE, 
                 data = data.clean, 
                 family = binomial)

AGE BYDECADE
43     40
45     40
24     20
21     20
23     20
37     30
34     30
36     30

以下有效，但显然我没有捕捉到所有 50 岁以上的人的十年

fit <- glm(COMPLICATION~AGE, 
             data = subset(data.clean, AGE > 50), 
             family = binomial)

但以下对我不起作用，我不知道如何去做。

fit <- glm(COMPLICATION~AGE, 
             data = subset(data.clean, AGE > 50 & AGE <= 60), 
             family = binomial)

Answer 1

这里有很多关于拟合分层 (G) LM 的问题。这是一种方法。

## convert AGE back to numeric:
data.clean <- transform(data.clean, AGE=as.factor(as.character(AGE)))
fits <- lme4::lmList(COMPLICATION~AGE|BYDECADE, 
                 data = data.clean, 
                 family = binomial)

分层逻辑回归

1 个答案: