Hello,

I am hoping to get some help on how to program a particular model in pymc. It is a two stage logistic regression hierarchical cascade where:

Event A has a probability of occurring and if it does, then there is some probability of Event B happening. If Event A does not happen then there is no possibility of Event B happening.

I have samples where B has occurred, where A has occurred but not B and where A has not occurred. One thing to note is that this is not a balanced dataset

Some ideas I had:

- The second level, B, probability is the multiplication of A’s logistic probability and B’s logistic probability
- The second level’s training data should only be those events which A was observed

Any help would be greatly appreciated. Thank you.