fitThresh function

Fits multi-state threshold model using the discrete approximation