x <- seq(0.1, 2.0, 0.2)
dpois.gpz(x, 2.0, 0.5, 0.2)
ppois.gpz(x, 2.0, 0.5, 0.2)
qpois.gpz(0.5, 2.0, 0.5, 0.2)
rpois.gpz(10, 2.0, 0.5, 0.2)
hpois.gpz(x, 2.0, 0.5, 0.2)
# Data
x <- stress
# ML estimates
params = list(alpha=0.2211, beta=0.6540, lambda=6.5456)
#P–P (probability–probability) plot
pp.plot(x, params = params, pfun = ppois.gpz, fit.line=TRUE)
#Q-Q (quantile–quantile) plot
qq.plot(x, params = params, qfun = qpois.gpz, fit.line=TRUE)
# Goodness-of-Fit(GoF) and Model Diagnostics
out <- gofic(x, params = params,
dfun = dpois.gpz, pfun=ppois.gpz, plot=FALSE)
print.gofic(out)
Run the code above in your browser using DataLab