x <- seq(0.1, 1, 0.1)
dHL.nhe(x, 1.5, 0.8, 2)
pHL.nhe(x, 1.5, 0.8, 2)
qHL.nhe(0.5, 1.5, 0.8, 2)
rHL.nhe(10, 1.5, 0.8, 2)
hHL.nhe(x, 1.5, 0.8, 2)
#Data
x <- windshield
#ML Estimates
params = list(alpha =0.1649, beta=3.7152, lambda=0.5881)
#P–P (probability–probability) plot
pp.plot(x, params = params, pfun = pHL.nhe, fit.line=TRUE)
#Q-Q (quantile–quantile) plot
qq.plot(x, params = params, qfun = qHL.nhe, fit.line=TRUE)
# Goodness-of-Fit(GoF) and Model Diagnostics
out <- gofic(x, params = params,
dfun = dHL.nhe, pfun = pHL.nhe, plot=FALSE)
print.gofic(out)
Run the code above in your browser using DataLab