plotPosteriorFit now returns theta value