man/examples/example_31_norm_dmgs.R

In fitdistcp: Distribution Fitting with Calibrating Priors for Commonly Used Distributions

#
# example 1
x=fitdistcp::d30norm_example_data_v1
p=c(1:9)/10
q=qnorm_dmgs_cp(x,p)
xmin=min(q$ml_quantiles,q$cp_quantiles);
xmax=max(q$ml_quantiles,q$cp_quantiles);
plot(q$ml_quantiles,p,xlab="quantile estimates",xlim=c(xmin,xmax),
	sub="(from qnorm_dmgs_cp)",
	main="Normal_DMGS: quantile estimates");
points(q$cp_quantiles,p,col="red",lwd=2)