fit <- lm(weight ~ Time + Diet, data = ChickWeight)
print_anova(anova(fit))
print_anova(car::Anova(fit, type = "III"))
# Underlying numeric p-values are untouched:
a <- anova(fit)
print_anova(a)
a[["Pr(>F)"]] # full-precision doubles
Run the code above in your browser using DataLab