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

Question

我需要表征数据的偏斜正态分布的位置、形状和规模，然后使用这些参数从具有相同参数的偏斜正常分布中随机抽取值。在过去，我使用R中的sn包来完成这项工作。因此，例如，如果我有一个数据向量v，它遵循偏斜正态分布，我会使用sn.em函数来估计位置、形状和规模：

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

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

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

sn包中不再存在这两个函数如何使用sn包中的不同功能或不同的包来完成此操作？

Answer 1

在"sn"包中使用sn.mmple()函数。例如，您可以使用

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分布中绘制示例，请在同一个包中使用rssn()函数。例如，

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