perl - perl中的随机二项分布数据集

Question

我尝试在 perl 中完成我在 R 中取得的成功，但很难与我的下游需求相结合。

在 R 中我做了以下操作

library("MASS")
d <- rnegbin(100000, mu = 250, theta = 2)
hist(d, breaks=1000, xlim=c(0,1000))

生成我需要的漂亮图表，峰值在 180-200 之间，尾部在右边。

有人可以帮我使用 Math::Random

编写 perl 等效代码吗

我试过了，但没有得到正确的形状

use Math::Random qw(random_negative_binomial);

# random_negative_binomial($n, $ne, $p)
# When called in an array context, returns an array of $n outcomes 
# generated from the negative binomial distribution with number of 
# events $ne and probability of an event in each trial $p. 
# When called in a scalar context, generates and returns only one 
# such outcome as a scalar, regardless of the value of $n.
# Argument restrictions: $ne is rounded using int(), the result must be positive. 
# $p must be between 0 and 1 exclusive.

# I tried different variable values but never got the right shape
my @dist = random_negative_binomial($n, $ne, $p);

我需要什么值来模拟 R 结果？我需要 X 上相同的值范围和相同的一般形状

感谢任何帮助，我没有找到那个包的插图示例

史蒂芬

Answer 1

我对统计了解不多，但由于没有其他人挺身而出:我会使用 Perl 数据语言 PDL(我将其用于其他事情)并获取 PDL::Stats::Distr 模块。您可以在此处找到与您的示例有些相似的示例 http://pdl-stats.sourceforge.net/Distr.htm .该模块包括 pmf_binomial(质量函数)和 mme_binomial(分布)。您还需要 PGPLOT 模块。

您将需要一些随机数据:

    $data = pdl  1..100000;    ## generate linear 1 - 100000
    $data = $data->random;     ## make them random between 0..1