``` M <- lm(height ~ sex * age * ntiles(age, n=10, format = "interval"), data = NCHS) gmodel(M) # doesn't work ```