<div dir="ltr">The way to do this is to use a differentiable logistic as Ben suggested. 
In addition you should use phases and gradually
<br>increase the steepness.  Finally as a technical matter you should use 
exp(z) in place of exp(-z) to avoid overflow when
<br>necessary.  See the attached changepoint.tpl
<br>
<br>To run this with multiple phase you would use the command line option 
-maxph 3
<br>
<br>Note use of value() function which will be needed in random effects model.
</div>