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

Question

我的年龄是我材料中的协变量。连续变量。年龄在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。