def create_z_row(power_set_i, p):
    """
    Create a row of Z for a given subset of the power set of p

    @param power_set_i: the ith row of the power set
    @param p: the number of covariates

    @return: a vector of 0's and 1's (1 corresponds to index in s)
    """
    return [int(x in power_set_i) for x in range(p)]