x <- seq(5, 10, 0.2)
dlind.ginv.exp(x, 5.0, 1.5, 0.5)
plind.ginv.exp(x, 5.0, 1.5, 0.5)
qlind.ginv.exp(0.5, 5.0, 1.5, 0.5)
rlind.ginv.exp(10, 5.0, 1.5, 0.5)
hlind.ginv.exp(x, 5.0, 1.5, 0.5)
# Data
x <- conductors
# ML estimates
params = list(alpha=97.0105, lambda=29.9324, theta=0.9028)
#P–P (probability–probability) plot
pp.plot(x, params = params, pfun = plind.ginv.exp, fit.line=TRUE)
#Q-Q (quantile–quantile) plot
qq.plot(x, params = params, qfun = qlind.ginv.exp, fit.line=TRUE)
# Goodness-of-Fit(GoF) and Model Diagnostics
out <- gofic(x, params = params,
dfun = dlind.ginv.exp, pfun=plind.ginv.exp, plot=FALSE)
print.gofic(out)
Run the code above in your browser using DataLab