Diagnostic Plot for Basic Local Independence Models

Description

Plots BLIM residuals against fitted values.

Usage

## S3 method for class 'blim'
plot(x, xlab = "Predicted response probabilities",
     ylab = "Deviance residuals", ...)

Arguments

x	an object of class blim, typically the result of a call to blim.
xlab, ylab, …	graphical parameters passed to plot.

Details

The deviance residuals are plotted against the predicted response probabilities for each response pattern.

See Also

blim, residuals.blim.

Examples

library("pks")

## Compare MD and MDML estimation

data(DoignonFalmagne7)
blim1 <- blim(DoignonFalmagne7$K, DoignonFalmagne7$N.R, method="MD")
blim2 <- blim(DoignonFalmagne7$K, DoignonFalmagne7$N.R, method="MDML")

par(mfrow = 1:2)      # residuals versus fitted values
plot(blim1, main = "MD estimation",   ylim = c(-4, 4))
plot(blim2, main = "MDML estimation", ylim = c(-4, 4))