Create a bernoulli distribution for the zi variable: Plot regression with python

Question

A set of the code in Python where I want to assign bernoulli values to the zi variable