绘制R中的偏态正态分布

Question

如果给出案例数，平均值，标准差，中位数和MAD，我如何在R中绘制偏态正态分布。

一个例子是我有1'196个案例，平均成本是6'389，标准差是5'158，中位数是4'930和MAD 1'366。而且我们知道收费的案例总是花费一些钱，所以成本必须始终是积极的。

我能找到的这个问题的最佳答案来自https://math.stackexchange.com/a/17995/54064，并建议使用sn包。但是我无法找到如何将它用于我的具体用例。

Answer 1

我在fGarch包中取得了一些成功。

require("fGarch")
hist(rsnorm(1000, mean = 0, sd = 1, xi = 15))

mmm <- replicate(300, {
  x <- rsnorm(1196, mean = 6389, sd = 5158, xi = 15)
  c(mean = mean(x), sd = sd(x))
})

> mean(mmm[1, ])
[1] 6404.312
> mean(mmm[2, ])
[1] 5169.572