eba (version 1.5-2)

plot.eba: Diagnostic Plot for EBA Models

Description

Plots EBA residuals against fitted values.

Usage

## S3 method for class 'eba':
plot(x, xlab = "Predicted choice probabilities",
  ylab = "Deviance residuals", \dots)

Arguments

x

an object of class eba, typically the result of a call to eba

xlab, ylab, ...

graphical parameters passed to plot.

Details

The deviance residuals are plotted against the predicted choice probabilities for the upper triangle of the paired-comparison matrix.

See Also

eba, residuals.eba.

Examples

## Compare two choice models

data(celebrities)                     # absolute choice frequencies
btl1 <- eba(celebrities)              # fit BTL model
A <- list(c(1,10), c(2,10), c(3,10), c(4,11), c(5,11), c(6,11),
          c(7,12), c(8,12), c(9,12))  # the structure of aspects
eba1 <- eba(celebrities, A)           # fit a preference tree (EBA)
anova(btl1, eba1)                     # model comparison based on likelihoods

par(mfrow = c(1,2))                   # residuals versus fitted values
plot(btl1, main = "BTL", ylim = c(-4, 4.5))  # BTL doesn't fit well
plot(eba1, main = "EBA", ylim = c(-4, 4.5))  # EBA fits better