Residuals for EBA Models

Computes deviance and Pearson residuals for eba objects.

# S3 method for eba
residuals(object, type = c("deviance", "pearson"), ...)

Arguments

object	an object of class eba, typically the result of a call to eba
type	the type of residuals which should be returned; the alternatives are: "deviance" (default) and "pearson"
...	further arguments passed to or from other methods; none are used in this method.

Details

Residuals are computed from the upper triangle of the paired-comparison matrix. See residuals.glm for details.

Value

A vector of residuals having as many elements as pairs of stimuli.

See also

eba, residuals.glm, plot.eba.

Examples

data(celebrities)                # absolute choice frequencies
btl1 <- eba(celebrities)         # fit Bradley-Terry-Luce model
sum( resid(btl1)^2 )             # deviance G2
#> [1] 78.21721
deviance(btl1)
#>   -2logL 
#> 78.21721 
sum( resid(btl1, "pearson")^2 )  # Pearson X2
#> [1] 77.24702