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