Compute p-values from full fit.
kronos_anova(fit, time)
fit
A lm model fit.
time
A string. Should be the column name containing the time values.