Categorical model with continuous dependent variable

Questions
4

If s only has binary value 0 and 1, you should model it as a Bernoulli:
obs = pm.Bernoulli('obs', p=p1, observed=s)