请教使用rantbl的一个问题

shiyiming

rantbl(seed,prob,....)这个格式能够用类似rantbl(seed,p)的形式替换么？（这里p是一个向量，元素是这个离散分布的概率）

我试验了我写的那个格式好像不能

那么怎么能把向量p放到rantbl里使用呢？

谢谢啦

shiyiming

[code:yk66m3kt]
data _null_;
        array _p{*} p1 - p5;
                array _s{5} _temporary_;
                do j=1 to dim(_p); _p[j]=j; end;               
                s=sum(of _p[*]);
                do j=1 to dim(_p); _p[j]=_p[j]/s; put _p[j]= 5.3  '   '   @;  _s[j]=0; end;
                put //;
                do i=1 to 1000;
                    x=rantbl(87554, of _p[*]);
                    _s[x]+1;
                end;
                do j=1 to dim(_p);
                    put _s[j]= 'out of 1000 ' /;
                end;
run;
[/code:yk66m3kt]

shiyiming

谢谢你的回复

这个使用array的方法，我在书上刚好找到。

不过我有个问题，就是我是用iml做的一些do loop。如果还用你写的在data下使用的方法就要把循环分开，把矩阵再放到data set里了。这样就有点太麻烦了，是不是在iml这个procedure下就没有方法能够把这个rantbl函数也用向量的方式表示离散的概率呢？

再次表示感谢

呵呵