请问产生p1后计算p的目的是什么？ 

  fij=1+1*treat+rj;p1=exp(fij)/(1+(exp(fij)));

  eij=rannor(0);ez=sqrt((1-p1)*p1)*eij;p=p1+ez;

  ez是在反映过度离散的情况吗，而且在模拟产生y的时候，p是没有被用到的

  y=ranbin(0,1,p1);

我是直接按照书上多水平logistic模型的公式编的， ez是水平1的误差

yij=（pai）ij+ez=exp(fij)/(1+(exp(fij)))+ez;

fij=1+1*treat+rj;

ez=sqrt((1-p1)*p1)*eij

其中，ez服从二项分布

由于自变量和因变量为非线性关系，所以第一水平误差项不用加，只添加第二水平误差项即可。