使用R中的sn包从偏正态分布中绘制数据

Question

我需要表征数据的偏态正态分布的位置，形状和比例，然后使用这些参数从具有相同参数的偏态正态分布中随机抽取值。 过去，我是使用R中的sn包来完成此操作的。因此，例如，如果我有一个数据向量v ，它遵循偏正态分布，那么我将使用sn.em估计位置，形状和比例功能：

require(sn)
data(ais)
v <- ais$Fe
sn.em(,v)

然后，我将使用rsnorm函数从具有相同参数的偏态正常dsitrubtion中随机抽取：

rsnorm(100, shape = x, location = y, scale = z)

这两个功能在sn软件包中不再存在。 我该如何使用sn软件包中的不同功能或一起使用不同的功能来做到这一点？

Answer 1

在“ sn”包中使用sn.mple（）函数。 例如，您可以使用

cp.est <- sn.mple(y=v,opt.method = "nlminb")$cp 
dp.est <- cp2dp(cp.est,family="SN")
dp.est
       xi     omega     alpha 
20.244158 73.840301  9.142412 

要从SN分布中抽取样本，请在同一软件包中使用rsn（）函数。 例如，

rsn(n=100, xi=20.24, omega=73.84, alpha=9.14)