I need to generate random samples from the Positive Stable (PS) density, and also compute its quantile (inverse cdf). For your convenience, I attach the corresponding paper (see bottom right of Page 9), which has the PS density and the Laplace transform. Then, I need to use it inside Proc NLMIXED as a random effect (which only allows Normal). Maybe I can rely on the Probability Integral Transformation (PIT) technique (not sure)? For that, I need the inverse cdf of the PS first. Is this do-able by Proc IML? If that is possible, how can I connect IML with NLMIXED?
PS is a somewhat crazy density. R has random sample draws and quantile functions of a particular parameterization of PS in the package 'stabledist', but there is nothing I see in SAS! Thoughts?
The authors of the paper state that "The method can be implemented using SAS IML," so I suggest you contact the authors to see if they will share some of their code. I suspect what they did was to use the NLP functions in SAS/IML to fit the parameters of this problem (maybe an MLE formulation?). They mention that they also called PROC PHREG as part of the analysis. They might have done that from inside their IML program, since you can call SAS procedures from IML
I doubt that they called an IML function from PROC NLMIXED. And the RANDOM statement only supports normally distributed effects, so using NLMIXED doesn't sound like a fruitful method.
Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.
If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website.