Geeglm y ~ x1 + I(x^2)

Question

En glm, si se especifica un modelo como

geeglm(qsmk ~ sex + race + age + I(edad*edad), 
             family=binomial, data=data,
             id=id, corstr="independence")

¿Qué indica la función I, I(edad^2) o I(edad*edad)?

Answer 1

I() significa "tal cual" (ver ?I): esto es necesario porque operadores como ^ y * tienen diferentes significados en un contexto de fórmula que cuando se están realizando cálculos regulares. I(age^2) y I(age*age) son equivalentes; ambos significan agregar un término "edad al cuadrado" al modelo. Si la fórmula contuviera términos no protegidos age^2 o age*age, denotarían la interacción de edad consigo misma (o, "el efecto principal de edad más la interacción de edad consigo misma"). La interacción de cualquier variable consigo misma se define como la variable original (es decir, age^2 es equivalente a edad en un contexto de fórmula; lo mismo sucede con edad*edad, porque edad+edad^2 -> edad+edad -> edad).