Threshold heteroskedastic models