x <- seq(2, 5, 0.25)
dexpo.inv.chen(x, 0.5, 2.5, 1.5)
pexpo.inv.chen(x, 0.5, 2.5, 1.5)
qexpo.inv.chen(0.5, 0.5, 2.5, 1.5)
rexpo.inv.chen(10, 0.5, 2.5, 1.5)
hexpo.inv.chen(x, 0.5, 2.5, 1.5)
# Data
x <- headneck44
# ML estimates
params = list(alpha=0.3947, lambda=15.5330, theta=8.1726)
#P–P (probability–probability) plot
pp.plot(x, params = params, pfun = pexpo.inv.chen, fit.line=TRUE)
#Q-Q (quantile–quantile) plot
qq.plot(x, params = params, qfun = qexpo.inv.chen, fit.line=TRUE)
# Goodness-of-Fit(GoF) and Model Diagnostics
# Display plot and print numerical summary
gofic(x, params = params,
dfun = dexpo.inv.chen, pfun=pexpo.inv.chen, plot=TRUE, verbose = TRUE)
Run the code above in your browser using DataLab