Two plots for the non-randomized PIT are currently available for checking the distributional assumption of the fitted iZIP model: the PIT histogram, and the uniform Q-Q plot for PIT.
histizipPIT( object, bins = 10, line = TRUE, colLine = "red", colHist = "royal blue", lwdLine = 2, main = NULL, ... ) qqizipPIT( object, bins = 10, col1 = "red", col2 = "black", lty1 = 1, lty2 = 2, type = "l", main = NULL, ... )
object | an object class "izip", obtained from a call to |
---|---|
bins | numeric; the number of bins shown in the PIT histogram or the PIT Q-Q plot. |
line | logical; if |
colLine | numeric or character: the colour of the line for comparison in PIT histogram. |
colHist | numeric or character; the colour of the histogram for PIT. |
lwdLine | numeric; the line widths for the comparison line in PIT histogram. |
main | character string; a main title for the plot. |
... | other arguments passed to plot.default and plot.ts. |
col1 | numeric or character; the colour of the sample uniform Q-Q plot in PIT. |
col2 | numeric or character; the colour of the theoretical uniform Q-Q plot in PIT. |
lty1 | integer or character string: the line types for the sample uniform Q-Q plot in PIT, see par(lty = .). |
lty2 | an integer or character string: the line types for the theoretical uniform Q-Q plot in PIT, see par(lty = .). |
type | 1-character string; the type of plot for the sample uniform Q-Q plot in PIT. |
The histogram and the Q-Q plot are used to compare the fitted profile with a standard uniform distribution. If they match relatively well, it means the iZIP distribution is appropriate for the data.
The gg_histizipPIT
and gg_qqizipPIT
functions
would provide the same two plots but in ggplot format.
Czado, C., Gneiting, T. and Held, L. (2009). Predictive model assessment for count data. Biometrics, 65, 1254--1261.
Dunsmuir, W.T.M. and Scott, D.J. (2015). The glarma
Package for Observation-Driven
Time Series Regression of Counts. Journal of Statistical Software,
67, 1--36.
## For examples see example(plot.izip)