将变量表示为 R 中逻辑回归中的多项式

Question

我的 Material 中有年龄作为协变量。连续变量。年龄在18-70岁之间。

我正在研究逻辑回归，并决定将年龄表示为多项式。

在数据中我有 4021 个观察结果，所以我只是复制了一段给你看看它们是什么样子的:

头(两个)

   gender       passinggrade age    prog
1    man          FALSE      69     FRIST
2    man             NA      70     FRIST
3 woman             NA       65     FRIST
4 woman           TRUE       68      FRIST
5 woman             NA       65     NMFIK
6    man          FALSE      70     FRIST

我的模型；

mod.fit<-glm(passinggrade ~prog+gender+age,family=binomial,data=both)

summary(mod.fit)

所以我想知道的是，我应该如何将年龄视为多项式？不知道我是否需要更改我的 R 代码中的某些内容？没有在 R 中做任何事情来“使年龄成为多项式”，我的问题很简单；你是怎么做到的？

Answer 1

您可以通过几种不同的方式执行此操作:

glm(passinggrade ~ prog + gender + poly(age, 3), ...

# Less preferred...
glm(passinggrade ~ prog + gender + age + I(age^2) + I(age^3), ...
glm(passinggrade ~ prog + gender + cbind(age, age^2, age^3), ...

参见 this post获取更多信息和讨论。