data a (type=corr);
_type_='corr';
input x1-x3;
cards;
1   .   .
.3  1   .
.5 .4   1
;
proc factor n=3;
run;

proc iml;
f={    0.78114        -0.46758         0.41375,
    0.70560         0.68785         0.17028,
    0.83518        -0.14380        -0.53084 };
data=rannor(j(5000,3,0));
z = data * t(f);     /*impose inner correlation*/
x1=z[,1]*.1 + 5; /*linear transformation*/
x2=z[,2]*.2 + 3;
x3=z[,3];
z=x1||x2||x3;
create dta1 from z [colname={x1 x2 x3}];
append from z;
quit;

proc corr data=dta1;
var x1 x2 x3;
run;

非常感谢你的回复，不过我想得到的是服从一个无序多分类的多项分布随机数，不知道你是否知道呢？

x=rand('TABLE',0.2,0.3,0.4,0.1);
产生数字1,2,3,4 with probability as above
Jingju

非常感谢！我用这个程序可以产生多项分布