View source: R/check_normality.R
check_normality | R Documentation |
Check model for (non-)normality of residuals.
check_normality(x, ...)
## S3 method for class 'merMod'
check_normality(x, effects = c("fixed", "random"), ...)
x |
A model object. |
... |
Currently not used. |
effects |
Should normality for residuals ( |
check_normality()
calls stats::shapiro.test
and checks the
standardized residuals (or studentized residuals for mixed models) for
normal distribution. Note that this formal test almost always yields
significant results for the distribution of residuals and visual inspection
(e.g. Q-Q plots) are preferable. For generalized linear models, no formal
statistical test is carried out. Rather, there's only a plot()
method for
GLMs. This plot shows a half-normal Q-Q plot of the absolute value of the
standardized deviance residuals is shown (in line with changes in
plot.lm()
for R 4.3+).
The p-value of the test statistics. A p-value < 0.05 indicates a significant deviation from normal distribution.
For mixed-effects models, studentized residuals, and not
standardized residuals, are used for the test. There is also a
plot()
-method
implemented in the see-package.
see::plot.see_check_normality()
for options to customize the plot.
m <<- lm(mpg ~ wt + cyl + gear + disp, data = mtcars)
check_normality(m)
# plot results
x <- check_normality(m)
plot(x)
# QQ-plot
plot(check_normality(m), type = "qq")
# PP-plot
plot(check_normality(m), type = "pp")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.