如何用已知参数绘制截断正态分布

Question

我想生成一个截断的正态分布，其中任何一个R中都有已知的参数。注意，我不是在寻找伪随机数生成器。

假设I有一个均值为5，标准差为1的正态分布。我能画出一个截断正态分布的值吗，在点1和点10截断？

Answer 1

下面的函数为任何分布生成（伪）随机数，只要您有d'和'p'函数可用。在R的说法中，这意味着您必须具有dnorm和pnorm，它们是基本R的一部分。

#
# random numbers for truncated distribution
#
rtrunc <- function(n, distr, lower = -Inf, upper = Inf, ...){
    makefun <- function(prefix, FUN, ...){
        txt <- paste(prefix, FUN, "(x, ...)", sep = "")
        function(x, ...) eval(parse(text = txt))
    }
    if(length(n) > 1) n <- length(n)
    pfun <- makefun("p", distr, ...)
    qfun <- makefun("q", distr, ...)
    lo <- pfun(lower, ...)
    up <- pfun(upper, ...)
    u <- runif(n, lo, up)
    qfun(u, ...)
}

# Example:
x <- rtrunc(1, "norm", lower = 0, mean = 2, sd = 5)