betareg (version 1.0-1)

plot.betareg: Plot Diagnostics for an betareg Object

Description

This function returns four plots: a plot of residuals against fitted values, a plot of standardized residuals against fitted values, a generalized leverage plot against fitted values and a plot of Cook's distances versus row labels.

Usage

plot.betareg(x, which = 1:4, caption = c("Deviance residuals vs indices of obs.", 
    "Standardized residuals vs indices of obs.", "Generalized leverage vs. Predicted values", "Cook's distance plot"), 
    panel = points, sub.caption = deparse(x$call), main = "", 
    ask = prod(par("mfcol")) < length(which) && dev.interactive(), 
    ..., id.n = 3, labels.id = names(residuals(x)), cex.id = 0.75)

Arguments

x

Fitted model by betareg.

which

If a subset of the plots is required, specify a subset of the numbers '1:4'.

caption

Captions to appear above the plots.

panel

Panel function. A useful alternative to 'points' is 'panel.smooth'.

sub.caption

common title-above figures if there are multiple; used as 'sub' (s.'title') otherwise.

main

title to each plot-in addition to the above 'caption'.

ask

logical; if `TRUE', the user is asked before each plot, see `par(ask=.)'.

...

optional arguments.

id.n

number of points to be labelled in each plot, starting with the most extreme.

labels.id

vector of labels, from which the labels for extreme points will be chosen. 'NULL' uses observation numbers.

cex.id

magnification of point labels.

References

FERRARI, S.L.P., CRIBARI-NETO, F. (2004). Beta regression for modeling rates and proportions. Journal of Applied Statistics, forthcoming.

See Also

betareg, br.fit, predict.betareg, residuals.betareg

Examples

data(pratergrouped)
attach(pratergrouped)
y <- V11/100
X <- cbind(V1,V2,V3,V4,V5,V6,V7,V8,V9,V10)
fit <- betareg(y ~ X)
par(mfrow=c(2,2))
plot(fit)