What is the exact non-linear model that you want to fit to the betas? One option is to use -stim_times_IM to generate a set of betas, one per event, and then fit a model to that. That is easier than setting up a new model for 3dNLfim. If there are too many events, particularly with other conditions, then IM is not a great option.
But in any case, can you describe the model you would like to fit to?