Skip to contents

Formula: DFBETAS_ij = (b_j - b_(i)j) / (s_(i) sqrt(C_jj)), C = (X'X)^-1

Usage

Dfbetas(X, y, intercept = TRUE)

Arguments

X

See Usage.

y

See Usage.

intercept

See Usage.

Value

List with dfbetas (n by p), cutoff, leverage, beta, sigma_i, n, p.

References

Belsley, Kuh and Welsch (1980), Regression Diagnostics: Identifying Influential Data and Sources of Collinearity, Wiley. The definitions and the cutoffs used here (2/sqrt(n) for DFBETAS, 2 sqrt(p/n) for DFFITS, scaling by the delete-one root mean square s_(i)) are as documented by the SAS and R reference implementations, which cite BKW for them.

Examples

V <- c(1, 2, 3, 4, 5, 6, 7, 8)
Dfbetas(V, V)
#> $dfbetas
#>             [,1]        [,2]
#> [1,] -0.18370455  0.15461042
#> [2,]  0.00000000  0.00000000
#> [3,]  0.15042309 -0.09494992
#> [4,]  0.18497895 -0.06227319
#> [5,]  0.07399158  0.06227319
#> [6,] -0.03926958  0.19830164
#> [7,] -0.46864708  0.98606323
#> [8,] -1.16533692  1.96155433
#> 
#> $cutoff
#> [1] 0.7071068
#> 
#> $leverage
#> [1] 0.4166667 0.2738095 0.1785714 0.1309524 0.1309524 0.1785714 0.2738095
#> [8] 0.4166667
#> 
#> $beta
#> [1] 6.097555e-16 1.000000e+00
#> 
#> $sigma_i
#> [1] 1.329622e-15 1.335963e-15 1.317870e-15 1.266196e-15 1.266196e-15
#> [6] 1.262033e-15 9.569522e-16 8.384101e-16
#> 
#> $n
#> [1] 8
#> 
#> $p
#> [1] 2
#> 
#> $method
#> [1] "DFBETAS (Belsley-Kuh-Welsch)"
#> 
#> attr(,"class")
#> [1] "morie_rich_result" "list"