View source: R/ols-pure-error-anova.R
ols_pure_error_anova | R Documentation |
Assess how much of the error in prediction is due to lack of model fit.
ols_pure_error_anova(model, ...)
model |
An object of class |
... |
Other parameters. |
The residual sum of squares resulting from a regression can be decomposed into 2 components:
Due to lack of fit
Due to random variation
If most of the error is due to lack of fit and not just random error, the model should be discarded and a new model must be built.
ols_pure_error_anova
returns an object of class
"ols_pure_error_anova"
. An object of class "ols_pure_error_anova"
is a
list containing the following components:
lackoffit |
lack of fit sum of squares |
pure_error |
pure error sum of squares |
rss |
regression sum of squares |
ess |
error sum of squares |
total |
total sum of squares |
rms |
regression mean square |
ems |
error mean square |
lms |
lack of fit mean square |
pms |
pure error mean square |
rf |
f statistic |
lf |
lack of fit f statistic |
pr |
p-value of f statistic |
pl |
p-value pf lack of fit f statistic |
mpred |
|
df_rss |
regression sum of squares degrees of freedom |
df_ess |
error sum of squares degrees of freedom |
df_lof |
lack of fit degrees of freedom |
df_error |
pure error degrees of freedom |
final |
data.frame; contains computed values used for the lack of fit f test |
resp |
character vector; name of |
preds |
character vector; name of |
The lack of fit F test works only with simple linear regression. Moreover, it is important that the data contains repeat observations i.e. replicates for at least one of the values of the predictor x. This test generally only applies to datasets with plenty of replicates.
Kutner, MH, Nachtscheim CJ, Neter J and Li W., 2004, Applied Linear Statistical Models (5th edition). Chicago, IL., McGraw Hill/Irwin.
model <- lm(mpg ~ disp, data = mtcars)
ols_pure_error_anova(model)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.