Is it fair/legitimate to use least square or optimization estimated parameter values to initialize the chain?

[monnahc]

FYI there is another similar discussion that might be of interest to you:

Overdispersed initial values -- general questions

I know the general advice is to initialize Stan (or any MCMC algorithm) from widely dispersed values relative to the posterior. My understanding is this helps catch issues like multimodality or chains not mixing well, which might not occur if started from the same point. But if the idea is to find local minima, why not use the optimizer? Start from the same dispersed values and optimize and check to see they all find the same mode (assuming for now it exists). It seems much more computationally…