Since almost any user will have multiple samples, and since the fitting model for one sample is just a special case of the multiple-sample one, I would suggest removing the sigfit.stan model and using sigfit_multi.stan always (perhaps renaming it to "sigfit_fitting.stan"), as well as removing the "multi" argument in the fit_signatures() function. That would also make the structure of the stanfit object coming from signature fitting consistent regardless of number of samples (always tri-dimensional).
Implemented, but needs testing!
The new model for fit_signatures() – with hierchical=FALSE – which is the previoius sigfit_multi.stan, is now called sigfit_fit.stan.
Since almost any user will have multiple samples, and since the fitting model for one sample is just a special case of the multiple-sample one, I would suggest removing the sigfit.stan model and using sigfit_multi.stan always (perhaps renaming it to "sigfit_fitting.stan"), as well as removing the "multi" argument in the fit_signatures() function. That would also make the structure of the stanfit object coming from signature fitting consistent regardless of number of samples (always tri-dimensional).