Association between Modic type and outcomes was explored with a mixed model and then by two-way analysis of variance ( ANOVA ) at each time point with Modic and treatment groups ( surgical , n = 126 ; nonsurgical , n = 117 ) as fixed factors , adjusted for disc degeneration , age , sex , smoking , and duration of radicular pain .