为什么rantbl()对cards里的分数概率不认识？

2430

收藏 2013-05-02

程序如下：data e;
input p1 p2;
do i=1 to 16;
a=rantbl(100,p1,p2);
end;
cards;
1/2 1/2
;
run;

日志里显示数据行“有对“p1”无效的数据”，不读。
但是直接用 a=rantbl(100,1/2,1/2);就可以。
因为我的概率有40多项，所以想用cards输入，请问高手这个怎么解决呢？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

墨岚

2013-5-2 16:21:18

问题不在rantbl,而是datalines/cards不认识这种分数表示。。那该怎么表示啊？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jingju11

2013-5-3 08:11:27

墨岚发表于 2013-5-2 16:21
问题不在rantbl,而是datalines/cards不认识这种分数表示。。那该怎么表示啊？

One more reason for the problem, is the function won't take characters as its arguments although you can read the fractions from cards. In fact, a more elegant way to do that is to load/write the values to a macro variable. That will give you a lot more flexibility on programing. Jingju

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

墨岚

2013-5-3 10:17:25

jingju11 发表于 2013-5-3 08:11
One more reason for the problem, is the function won't take characters as its arguments although y ...

Thanks! The idea that input fractions as characters solves my problem indeed('_temorary_' is very useful!). But I recognize code like "a =rantbl(123,.1,1/2,1/3,1/15);" can be read as well .
By the way, can rantbl() generate several radom numbers, instead of one number at a time? I need different numbers in one repeat...

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jingju11

2013-5-3 10:31:04

Yeah, it takes fractions (1/2)but not the character ('1/2').
I don't think think so. as you know, sas function returns one value at a time. If you want to have multiple columns, just insert one more loop, like:

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

墨岚

2013-5-3 10:43:52

jingju11 发表于 2013-5-3 10:31
Yeah, it takes fractions (1/2)but not the character ('1/2').
I don't think think so. as you know, ...

That's what I did.. Unfortunately, these three numbers are in fact three samples but not three numbers in one sample which is preferred. Still thank u for the confirmation!

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

jingju11

2013-5-3 10:56:12

random numbers not distinguished by samples. i tend to understand your samples as different datasets. using output in different positions will give you different forms of data...

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

墨岚

2013-5-3 11:29:13

jingju11 发表于 2013-5-3 10:56
random numbers not distinguished by samples. i tend to understand your samples as different datasets ...

比如现在有一个20项的多项分布，按照这个分布取得一批随机数不是我的最终目的，这批随机数实际上是我想做出的另一个分布的一个样本，因为这个样本是符合前面的多项分布的。每次抽样就是按照这个分布抽取固定数目的随机数比如10个，因为模型的原因这10个数不能重复，类似于一次不放回的抽样，所以我想问可否一次抽取多个不重复的随机数..rantbl的sas帮助里也没有看到这方面的说明.

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yongyitian

2013-5-3 12:09:15

There is a tricky way to avoid your input fraction problem. It may also give you more options to programm with array or loop.

Or just type a few more letters and numbers, it could save your time working out how the macro or array, or loops could make the program looks pretty.
data e;
input n1 d1 n2 d2 n3 d3 n4 d4 n5 d5;
   p1 = n1/d1;
   p2 = n2/d2;
   p3 = n3/d3;
   p4 = n4/d4;
   p5 = n5/d5;
   do i=1 to 16;
a=rantbl(100,p1,p2,p3,p4,p5);
output;
end;
datalines;
1 2 1 3 1 4 1 5 1 6
;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

墨岚

2013-5-3 12:26:05

yongyitian 发表于 2013-5-3 12:09
There is a tricky way to avoid your input fraction problem. It may also give you more options to pro ...

Tricky indeed.

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yongyitian

2013-5-3 12:44:19

墨岚发表于 2013-5-3 11:29
比如现在有一个20项的多项分布，按照这个分布取得一批随机数不是我的最终目的，这批随机数实际上是我想做 ...

use different seed to get a different random number.

do seed = 100 to 105
   a = rantbl(seed, p1,p2,p3,p4,p5);
      output;
   end;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jingju11

2013-5-11 08:26:01

墨岚发表于 2013-5-3 11:29
比如现在有一个20项的多项分布，按照这个分布取得一批随机数不是我的最终目的，这批随机数实际上是我想做 ...

非常复杂。没有理解。我想，table分布抓出来的是球的颜色编号。而多项分布是不同颜色球的数量组合的分布。我想对于任何一维的随机数，一次抓10个和每次一个抓十次没有任何区别。这应该也是随机数的核心之所在。我的一个小程序关于multinomial distribution.

京剧

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

墨岚

2013-5-12 20:36:07

谢谢！随机数的确应该那样产生。看来我的确把rantbl用错了。你给的博客很有用~

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

墨岚

2013-5-12 20:37:55

jingju11 发表于 2013-5-11 08:26
非常复杂。没有理解。我想，table分布抓出来的是球的颜色编号。而多项分布是不同颜色球的数量组合的分布。 ...

谢谢！随机数的确应该那样产生。看来我的确把rantbl用错了。你给的博客很有用~

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

墨岚

2013-5-12 23:58:08

jingju11 发表于 2013-5-11 08:26
非常复杂。没有理解。我想，table分布抓出来的是球的颜色编号。而多项分布是不同颜色球的数量组合的分布。 ...

刚才仔细看了一下你博客的那篇文章，对多项分布的解释和实现程序我都是认同的。这样我的情形可能就不是多项分布了，因为按你例子中的情形来说，我那n个球是一次性拿出的，没有补回。然后，，这n个球是相关的而不应该用产生n个随机数来模拟了。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群