Plot the quantile-quantile graph for censored data.
Plot the quantile-quantile graph for censored data.
This function may be rather slow for many iterations/many data because it relies on numerical inversion of the mixture Cumulative Distribution Function. set.seed(150520) data(salinity) out <- MixNRMI1cens(xleft = salinityleft,xright=salinityright, extras = TRUE, Nit = 100) BNPdensity:::qq_plot_censored(out)
qq_plot_censored(fit, thinning_to =500)
Arguments
fit: The result of the fit, obtained through the function MixNRMI1 or MixNRMI2, MixMRMI1cens or MixMRMI2cens
thinning_to: How many iterations to compute the mean posterior quantiles