Сравните две модели регрессии в R

Question

age25=subset(juul,juul[,"age"]>25.00)## create a subset of age greater than 25
modelgf=lm(age25[,"igf1"]~age25[,"age"])
age20=subset(juul,juul[,"age"]<20.00)
modelgf2=lm(age20[,"igf1"]~age20[,"age"])

Я пытался сравнить модели modelgf и modelgf2, используя anova(m1,m2).Тем не менее, я получаю предупреждение:

In anova.lmlist(object, ...) :
  models with response ‘"age20[, \"igf1\"]"’ removed because response differs from model 1

Есть ли другие способы сравнить эти две модели?

Answer 1

Вот, пожалуйста:

# Dummy for Age>25
juul[,"ageCat25"] <- juul[,"ageCat"] > 25.00
# Collinear dummy for Age<20
juul[,"ageCat20"] <- ifelse(!juul[,"ageCat25"] & juul[,"age"]<20.00, TRUE, juul[,"ageCat25"])
m1 <- lm(foo ~ ageCat25, juul)
m2 <- lm(foo ~ ageCat20, juul)
anova(m1,m2)

Интерпретация оставлена ​​на ФП.