r - 两片正态(或 split 正态)分布的密度

Question

两件式正态分布是否有密度函数:



在 CRAN 上？以为我会在编写代码之前检查一下。我已经检查了分发任务 View 。它没有在那里列出。我查看了几个可能的包，但无济于事。

更新:我添加了 dsplitnorm , psplitnorm , qsplitnorm和 rsplitnorm函数到 fanplot包裹。

Answer 1

如果您选择构建您自己的发行版本，您可能会对 感兴趣。地区 .它(以及相关的软件包 distrEx , distrSim , distrTEst , distrTeach 和 Doc)已写入 distr为从现有发行版构建新发行版提供统一的界面。 (我在 distrDoc 软件包附带的精彩插图的帮助下构建了这个示例，可以通过键入 vignette("distr") 获得。)

这实现了 split normal distribution ，这可能不是您所追求的。使用 地区 工具集，但是，调整它以适应您的确切需求应该不会太难。

library(distr)

## Construct the distribution object.
## Here, it's a split normal distribution with mode=0, and lower- and
## upper-half standard deviations of 1 and 2, respectively.
splitNorm <- UnivarMixingDistribution(Truncate(Norm(0,2), upper=0), 
                                      Truncate(Norm(0,1), lower=0), 
                                      mixCoeff=c(0.5, 0.5))
## Construct its density function ...
dsplitNorm <- d(splitNorm)
## ... and a function for sampling random variates from it
rsplitNorm <- r(splitNorm)

## Compare the density it returns to that from rnorm()
dsplitNorm(-1)
# [1] 0.1760327    
dnorm(-1, sd=2)
# [1] 0.1760327    

## Sample and plot a million random variates from the distribution
x <- rsplitNorm(1e6)         
hist(x, breaks=100, col="grey")

## Plot the distribution's continuous density
plot(splitNorm, to.draw.arg="d")