用R中的每个不同的df和ncp实现t个随机变量？

Question

我正在尝试使用 rt() 生成 5 个随机 t 变量，这 5 个变量中的每一个都有一个特定的 df（分别从 1 到 5 ) 和特定的ncp（分别为seq(0, 1, l = 5)）。因此，5 个随机 t 变量，每个变量都有不同的 df 和不同的 ncp。

为了实现上述目标，我尝试了以下，但没有成功。什么是实现我上面描述的高效 R 代码？

vec.rt = Vectorize(function(n, df, ncp) rt(n, df, ncp), c("n", "df", "ncp"))

vec.rt(n = 5, df = 1:5, ncp = seq(0, 1, l = 5))

或者

mapply(FUN = rt, n = 5 , df = 1:5, ncp = seq(0, 1, l = 5))

通知：

rt(n = 5, df = 1:5, ncp = seq(0, 1, l = 5))

R 给出以下警告：

Warning message: In if (is.na(ncp)) { : the condition has length > 1 and only the first element will be used

Answer 1

重新表述您的问题有助于找到答案：您需要来自 5 个随机变量的样本长度为 1 (n = 1)，每个变量都有不同的参数。

mapply(FUN = rt, n = 1 , df = 1:5, ncp = seq(0, 1, l = 5))