VGAM (version 0.7-1)

qtplot.gumbel: Quantile Plot for Gumbel Regression

Description

Plots quantiles associated with a Gumbel model.

Usage

qtplot.gumbel(object, plot.it = TRUE,
    y.arg = TRUE, spline.fit = FALSE, label = TRUE,
    R = object@misc$R, percentiles = object@misc$percentiles,
    add.arg = FALSE, mpv = object@misc$mpv,
    xlab = NULL, ylab = "", main = "",
    pch = par()$pch, pcol.arg = par()$col,
    llty.arg = par()$lty, lcol.arg = par()$col, llwd.arg = par()$lwd,
    tcol.arg = par()$col, tadj = 1, ...)

Value

• The object with a list called qtplot in the post slot of object. (If plot.it=FALSE then just the list is returned.) The list contains components
• fitted.values
  The percentiles of the response, possibly including the MPV.
• percentiles
  The percentiles (small vector of values between 0 and 100.

Details

There should be a single covariate such as time. The quantiles specified by percentiles are plotted.

See Also

gumbel.

Examples

data(venice)
y = as.matrix(venice[,paste("r",1:10,sep="")])
fit1 = vgam(y ~ s(year, df=3), gumbel(R=365, mpv=TRUE),
            data=venice, trace=TRUE, na.action=na.pass)
fitted(fit1)[1:4,]

par(mfrow=c(1,1), bty="l", xpd=TRUE, las=1)
qtplot(fit1, mpv=TRUE, lcol=c(1,2,5), tcol=c(1,2,5), lwd=2,
       pcol="blue", tadj=0.4)

qtplot(fit1, perc=97, mpv=FALSE, lcol=3, tcol=3,
       lwd=2, tadj=0.4, add=TRUE) -> i
i@post$qtplot$fitted[1:4,]