假设我有这样两个自变量x1 x2 都是0,1编码的,  现在我想通过logistic模型通过x1 x2与反应变量y建立关系，从而生成一组基于二项分布 B(N, P) 的反应变量（也是0,1编码的）
关系根据 logit(p)=b0 + b1*x1+ b2*x2 现在b0 b1 b2都可已知x1 x2为上面提到的已知0,1编码数据，样本量可随意。
logit(p) = ln(P/1-P)  P为阳性结果发生概率
例如：x1    x2      y
            0      0            
            0      1
            1      0
            0      1
            .....  ....
            ...    .....                                                                     
请教各路大神这样一个过程在R中如何实现呢？