Fitting multiple cones to 2D data, how to impose ordering constraint on parameters

One more related question, how you would you achieve the same effect for discrete variables? it seems that this transformation only works for continuous variables (I receive the some what vague error ValueError: Transformations for discrete distributions)