Hi @twiecki
I need to extend the context of my question.
I am working on the MMM solution and trying to implement the saturation function described in the paper “bayesian methods for media mix modeling with carryover and shape effects”.
The function has two parameters: slope and half saturation point K.
The function is directly applicable to the data on 0-1 interval. But will not work when data is not normalized or even when it goes through adstock transformation because it may cause adstock higher than 1. One of the solutions I have found is to internally renormalize parameter K according to the quantile of the input data but I am stuck with linspace and quantile
Same issue is with logistic function.
As am not the one who is using pymc3 for MMM I am sure there exist some other solution.