Added model fit conditional on expected p-value for 10x speedup.