我正在尝试生成 1000 组 130 个随机数,这些随机数符合以下数字的偏态正态分布:
-10.4, -9.3, -6.8, -4.8, -5.7, 5.8, -4.5, -3.4, -2, 0.3, -0.4, -4.1, -6.9, -5.9, -2.5, -2, -2.8, -3.2, -4.4, -2, -1.4, 0.9, -1, -4.1, -11.7, 0.1
这些数字的平均值为 -3.99,标准差为 3.17,偏度为 -0.71,峰度为 0.22。
为了获得 1000 组 130 个随机数,我尝试了以下方法:
install.packages("sn")
library(sn)
p <- rmsn(n = 130,
xi = rep(-3.99, 1000),
Omega = diag(1000),
alpha = rep(-0.71, 1000),
tau = -0.71)
我得到 1000 个由 130 个随机数组成的向量,平均值为 -3.99。但是,它们没有偏斜 -0.71,而且我不知道如何将标准差设置为 3.17 或将峰度设置为 0.22。
任何帮助将非常感激!