This function generates a data set for testing purposes. There is a single instrument that takes on values of 1 and 2, and the distribution of the values for the binary instrument is uniform. The MTRs are m0 ~ 1 and m1 ~ 1. All unobservables u are integrated out.
gendist3(subN = 5, p1 = 0.4, p2 = 0.6)
integer, default set to 5. This is the number of individuals possessing each value of the instrument. So the total number of observations is subN * 2.
the probability of treatment for those with the instrument Z = 1.
the probability of treatment for those with the instrument Z = 2.
a data.frame.