来自 R 中 Erlang 分布的样本

Question

的Erlang分布有两个参数：自然数ķ（形状参数）和实数拉姆达（速率参数）。你怎么能从Erlang 分布中随机抽取一个大小为n 的样本，使用 R 而不使用外来包？

Answer 1

这个问题的诀窍是要注意Erlang 分布是Gamma 分布的一个特例。来自 gamma 分布的采样在stats 包中实现。

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

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

n 是样本的大小。

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

rate 是 Erlang 分布的速率参数。