| qtplot.gumbel {VGAM} | R Documentation |
Plots quantiles associated with a Gumbel model.
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, ...)
object |
A VGAM extremes model of the
Gumbel type, produced by modelling functions such as vglm
and vgam with a family function either
"gumbel" or "egumbel". |
plot.it |
Logical. Plot it? If FALSE no plot will be done.
|
y.arg |
Logical. Add the raw data on to the plot? |
spline.fit |
Logical. Use a spline fit through the fitted percentiles? This can be useful if there are large gaps between some values along the covariate. |
label |
Logical. Label the percentiles? |
R |
See gumbel. |
percentiles |
See gumbel. |
add.arg |
Logical. Add the plot to an existing plot? |
mpv |
See gumbel. |
xlab |
Caption for the x-axis. See par. |
ylab |
Caption for the y-axis. See par. |
main |
Title of the plot. See title. |
pch |
Plotting character. See par. |
pcol.arg |
Color of the points.
See the col argument of par. |
llty.arg |
Line type. Line type.
See the lty argument of par. |
lcol.arg |
Color of the lines.
See the col argument of par.
|
llwd.arg |
Line width.
See the lwd argument of par.
|
tcol.arg |
Color of the text
(if label is TRUE).
See the col argument of par.
|
tadj |
Text justification.
See the adj argument of par.
|
... |
Arguments passed into the plot function
when setting up the entire plot. Useful arguments here include
sub and las.
|
There should be a single covariate such as time.
The quantiles specified by percentiles are plotted.
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. |
Unlike gumbel, one cannot have
percentiles=NULL.
Thomas W. Yee
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,]
## Not run:
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,]
## End(Not run)