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 |
plot.it |
Logical. Plot it? If |
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 |
percentiles |
See |
add.arg |
Logical. Add the plot to an existing plot? |
mpv |
See |
xlab |
Caption for the x-axis. See |
ylab |
Caption for the y-axis. See |
main |
Title of the plot. See |
pch |
Plotting character. See |
pcol.arg |
Color of the points.
See the |
llty.arg |
Line type. Line type.
See the |
lcol.arg |
Color of the lines.
See the |
llwd.arg |
Line width.
See the |
tcol.arg |
Color of the text
(if |
tadj |
Text justification.
See the |
... |
Arguments passed into the |
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
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) head(fitted(fit1)) ## 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 head(i@post$qtplot$fitted) ## End(Not run)