Lack of fit of the regression model.

This function performs a hypothesis test for lack of fit in a simple regression model.

lack_fit(mod)

Arguments

mod	A lm object.

Value

lack_fit function returns a list with the statistica and p value.

Examples

x <- c(1.0, 1.0, 2.0, 3.3, 3.3, 4.0, 4.0, 4.0, 4.7, 5.0,
       5.6, 5.6, 5.6, 6.0, 6.0, 6.5, 6.9)
y <- c(10.84, 9.30, 16.35, 22.88, 24.35, 24.56, 25.86,
       29.16, 24.59, 22.25, 25.90, 27.20, 25.61, 25.45,
       26.56, 21.03, 21.46)
plot(x=x, y=y, pch=20)

mod <- lm(y ~ x)
lack_fit(mod)
#> $SSep
#> [1] 15.56303
#> 
#> $SSres
#> [1] 250.1338
#> 
#> $SSfa
#> [1] 234.5708
#> 
#> $Fo
#> [1] 13.18827
#> 
#> $p.value
#> [1] 0.001388715
#>