r - 从 R 中的 t 分布采样

Question

当我想指定自由度 (df)、均值和方差时，如何从 R 中的 t 分布中采样？例如，我如何从 8 df、均值 = 4、方差 = 16 的 t 分布中采样？我想我必须以某种方式使用非中心性参数，但我不确定如何。

Answer 1

你想要非中心t，还是只是一个缩放和移动（中心）的 t？要从具有位置m和比例参数的分布中抽样s，请执行

rt(N, df=df)*s + m

at 分布的方差是df/(df-2)为了得到一个特定的方差v：

rt(N, df=df)*sqrt(v * (df-2)/df) + m

Answer 2

如果您查看 的代码rt，您会看到来自非中心 t 的随机变量从以下位置返回：

rnorm(n, ncp)/sqrt(rchisq(n, df)/df)

因此，非中心参数rnorm作为分子分布的均值给出。