x <- seq(0.1, 2.0, 0.1)
dlogis.chen(x, 1.5, 1.5, 0.1)
plogis.chen(x, 1.5, 1.5, 0.1)
qlogis.chen(0.5, 1.5, 1.5, 0.1)
rlogis.chen(10, 2.0, 5.0, 0.1)
hlogis.chen(x, 1.5, 1.5, 0.1)
# Data
x <- bladder
# ML estimates
params = list(alpha=4.46424, beta=0.15506, lambda=0.24904)
#P–P (probability–probability) plot
pp.plot(x, params = params, pfun = plogis.chen, fit.line=TRUE)
#Q-Q (quantile–quantile) plot
qq.plot(x, params = params, qfun = qlogis.chen, fit.line=TRUE)
# Goodness-of-Fit(GoF) and Model Diagnostics
out <- gofic(x, params = params,
dfun = dlogis.chen, pfun=plogis.chen, plot=FALSE)
print.gofic(out)
Run the code above in your browser using DataLab