x <- seq(0.1, 1, 0.1)
dinv.expo.exp.pois(x, 1.5, 0.8, 2)
pinv.expo.exp.pois(x, 1.5, 0.8, 2)
qinv.expo.exp.pois(0.5, 1.5, 0.8, 2)
rinv.expo.exp.pois(10, 1.5, 0.8, 2)
hinv.expo.exp.pois(x, 1.5, 0.8, 2)
#Data
x <- conductors
#ML Estimates
params = list(alpha =40.5895, beta=22.7519, lambda=2.9979)
#P–P (probability–probability) plot
pp.plot(x, params = params, pfun = pinv.expo.exp.pois, fit.line=TRUE)
#Q-Q (quantile–quantile) plot
qq.plot(x, params = params, qfun = qinv.expo.exp.pois, fit.line=TRUE)
# Goodness-of-Fit(GoF) and Model Diagnostics
res <- gofic(x, params = params, dfun = dinv.expo.exp.pois,
pfun = pinv.expo.exp.pois, plot=FALSE)
print.gofic(res)
Run the code above in your browser using DataLab