fit_lm <- lm(weight ~ Time + Diet, data = ChickWeight)
print_summary(fit_lm)
fit_lmer <- lme4::lmer(weight ~ Time + (1 | Chick), data = ChickWeight)
print_summary(fit_lmer)
# Underlying numeric p-values are untouched:
sm <- summary(fit_lm)
sm$coefficients[, "Pr(>|t|)"] # full-precision doubles
Run the code above in your browser using DataLab