SAS macro to generate randomized quantile residuals from PROC GLIMMIX?

[quentinread]

I am a big user of the R package DHARMa (https://cran.r-project.org/web/packages/DHARMa/vignettes/DHARMa.html) which produces scaled quantile residuals from fitted GLMMs and does a variety of diagnostic tests on them such as uniformity, overdispersion, and autocorrelation. 

I know that it is possible to write a macro in SAS to generate quantile residuals via simulation from a fitted GLIMMIX model and run the diagnostic tests. But I was curious if anyone in the community already has a macro that can do that, that they are willing to share.

Thanks in advance for your help!