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

我需要描述数据的偏正态分布的位置、形状和规模，然后使用这些参数从具有相同参数的偏正态分布中随机抽取值。在过去，我使用R中的

sn

包来实现这一点。因此，例如，如果我有一个数据向量

v

，它遵循斜正态分布，我将使用

sn.em

函数来估计位置、形状和比例：

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

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

这两个函数在

sn

包中不再存在。如何使用

sn

包中的不同函数或不同的包一起执行此操作？

在“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 

---

可能功能已被替换/重命名：您是否检查了@user20650？检查参考手册时运气不佳。似乎在任何地方都没有记录。快速浏览一下手册，可能是

rsn

（第20页）和

selm

（第45页）@user20650使用

rsn

功能，尽管不清楚

selm

中的输出对应于位置（

xi

）、刻度（

omega

）和倾斜（

alpha

）运行

rsn

功能所需的参数。是的，不确定，我只看到了selm，因为我注意到了线性模型&mle，乍一看，它似乎在正确的球场上。我所要做的就是阅读手册，确保sn.em没有被弃用/替换/改进，如果不清楚，请发电子邮件给维护人员，询问删除该功能的原因。（因为如果要使用函数，通常可以使用旧源加载旧的包版本）

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