x <- seq(0.1, 10, 0.2)
dlogis.exp.ext(x, 2.0, 5.0, 0.1)
plogis.exp.ext(x, 2.0, 5.0, 0.1)
qlogis.exp.ext(0.5, 2.0, 5.0, 0.1)
rlogis.exp.ext(10, 2.0, 5.0, 0.1)
hlogis.exp.ext(x, 2.0, 5.0, 0.1)
# Data
x <- stress31
# ML estimates
params = list(alpha=1.7919, beta=418.0473, lambda=0.1211)
#P–P (probability–probability) plot
pp.plot(x, params = params, pfun = plogis.exp.ext, fit.line=TRUE)
#Q-Q (quantile–quantile) plot
qq.plot(x, params = params, qfun = qlogis.exp.ext, fit.line=TRUE)
# Goodness-of-Fit(GoF) and Model Diagnostics
res <- gofic(x, params = params,
dfun = dlogis.exp.ext, pfun=plogis.exp.ext, plot=TRUE)
print.gofic(res)
Run the code above in your browser using DataLab