I am trying to set up the function using
np.random.standard t and np.random.binomial.
X|(Y = 1) is a t − distribution with 1 degree of freedom with mean µ1 X|(Y =