guass如何生成符合二项分布的序列

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

高斯如何生成符合二项分布的序列啊？找不到具体命令，只有负二项分布。。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：GUASS 二项分布 负二项分布 gauss 二项分布

new;cls;
p=0.4;
n=100;
r=100; @100行@
c=3; @生成3列100行二项式分布随机变量@
w=reshape(sumc((rndu(r*c,1).>(0|pbinom(seqa(0,1,n),p,n))')')-1,r,c);
print w;
proc pbinom(x,p,n); @累积概率分布@
local xx,d;
      xx = trunc(x).*(x.>=0); @ trunc取整或者删除小数@
      if maxc(maxc(xx)) ge 0;
         d = cumsumc(dbinom(seqa(0,1,maxc(maxc(xx))+1),p,n));
         d = reshape(d[xx+1],rows(x),cols(x));
      else;
      d = zeros(rows(x),cols(x));
      endif;
retp(d.*(x.>=0));
endp;
proc dbinom(x,p,n); @概率密度@
local xx,d;
    if p==0; d = (x.==0);
    elseif p==1; d = (x.==n);
    else;
    xx = x.*(x.==trunc(x)).*(x.>=0).*(x.<=n);
    d = exp(lnfact(n)-lnfact(xx)-lnfact(n-xx)+xx*ln(p)+(n-xx)*ln(1-p));
  endif; @ lnfact计算阶乘函数的自然对数@
retp(d.*(x.==trunc(x)).*(x.>=0).*(x.<=n));
endp;