I am new to stan and bayesian analysis and looking at modelling some spectroscopic data (UV-Vis), do you think for this kind of ‘signal’ data with discrete counts over a range of energies and a number of peaks (some which can be indistinguishable), a gaussian multivariate mixture model is appropriate?

Are there any resources dedicated to getting started with MV-GMM’s in Stan in any case?

I want to be able to get an estimate of parameters, and possibly be able to identify the number of peaks within the data (this is often indistinguishable), so that the absorbance can be identified. The parameters include the FWHM, amplitudes and any peak positions. Eventually this will also include a ‘background’ estimate for example.

it is my fault sorry I misunderstood your question. Y i.e. the excitation here can be modelled as gaussian function to account for the peak shape as a function of the wavelength(x). so e(v) ~ A*exp(-(x-xo)/fwhm)^2).
where xo= centre/peak position which is real and positive. x=wavelength values, again real and positive fwhm=peak width, real and positive and Amplitude is also real and positive.

I have tried fitting this in pystan using a general non-linear approach and although the fit is visually good with convergence. The parameter values have high levels of deviation. So I was looking at whether another model could better describe the data. But I think I may possibly be out of my depth here