x <- seq(1.0, 5, 0.25)
dhc.inv.NHE(x, 2.0, 0.5, 2.5)
phc.inv.NHE(x, 2.0, 0.5, 2.5)
qhc.inv.NHE(0.5, 2.0, 0.5, 2.5)
rhc.inv.NHE(10, 2.0, 0.5, 2.5)
hhc.inv.NHE(x, 2.0, 0.5, 2.5)
# Data
x <- relief
# ML estimates
params = list(beta=79.7799, lambda=0.1129, theta=154.1769)
#P–P (probability–probability) plot
pp.plot(x, params = params, pfun = phc.inv.NHE, fit.line=TRUE)
#Q-Q (quantile–quantile) plot
qq.plot(x, params = params, qfun = qhc.inv.NHE, fit.line=TRUE)
# Goodness-of-Fit(GoF) and Model Diagnostics
res <- gofic(x, params = params,
dfun = dhc.inv.NHE, pfun=phc.inv.NHE, plot=FALSE)
print.gofic(res)
Run the code above in your browser using DataLab