x <- seq(0.1, 10, 0.2)
dmod.inv.gen.exp(x, 2.0, 0.5, 0.2)
pmod.inv.gen.exp(x, 2.0, 0.5, 0.2)
qmod.inv.gen.exp(0.5, 2.0, 0.5, 0.2)
rmod.inv.gen.exp(10, 2.0, 0.5, 0.2)
hmod.inv.gen.exp(x, 2.0, 0.5, 0.2)
# Data
x <- fibers69
# ML estimates
params = list(alpha=30.7790, beta=0.1942, lambda=14.8297)
#P–P (probability–probability) plot
pp.plot(x, params = params, pfun = pmod.inv.gen.exp, fit.line=TRUE)
#Q-Q (quantile–quantile) plot
qq.plot(x, params = params, qfun = qmod.inv.gen.exp, fit.line=TRUE)
# Goodness-of-Fit(GoF) and Model Diagnostics
out <- gofic(x, params = params,
dfun = dmod.inv.gen.exp, pfun=pmod.inv.gen.exp, plot=TRUE)
print.gofic(out)
Run the code above in your browser using DataLab