做列联表分析，可以检测分类变量的相关性

我想做的是数据模拟，而不是分析。谢谢

proc iml;
N=50;T=4;
bt1=0.5;bt0=0;
seed=j(T,N,0);
mu=j(T,N,0);
cov=
{1.00  .70    .50    .30,
  .70  1.00    .60   .40,
  .50    .60  1.00   .50,
  .30    .40    .50 1.00  };
sig=t(root(cov));
Z=sig*rannor(seed)+mu;
X=ranuni(seed)*2-1;
Y=t(Z<=(bt0+bt1*X));
create y from y;
append from y;
quit;
这个帖子里，我解释的比较详细：https://bbs.pinggu.org/thread-2121268-1-1.html
做列联表分析
proc freq data=y;
table col1*col2/expected chisq nocol norow nopercent;
run;
可以看到生成的二分类变量有明显的相关性。

非常感谢，ziyenano的解答。那我怎样能够确定我的二分类变量的分布情况，比如我想让二分类变量服从二项分布B(n,0.2)。