来自 R 中 Erlang 分布的样本

Question

Erlang distribution 有两个参数：自然数k（形状参数）和实数lambda（速率参数）。如何在没有外来包的情况下使用 R 从 Erlang 分布中抽取大小为 n 的随机样本？

Answer 1

这个问题的诀窍是注意the Erlang distribution 是the Gamma distribution 的一个特例。从 gamma 分布采样在 stats package 中实现。

以下函数从标准 R 向量中的 Erlang(k, lambda) 分布返回大小为 n 的样本：

rgamma(n, shape=k, rate = lambda)

n 是样本的大小。

shape 是参数k。这是 Erlang 分布的形状参数（对于 Erlang，这必须是自然数 >=1）。

rate 是 Erlang 分布的速率参数。