既知のパラメータを持つ切り捨て正規分布をプロットする方法R

Question

Rのいずれかに既知のパラメータを持つ切り捨て正規分布を生成したいと思います。私は擬似乱数ジェネレータを求めていません。

平均5と標準偏差1の正規分布を仮定します。ポイント1と10で切り捨てられた切り捨て正規分布の値をプロットできますか？

Answer 1

この問題の解決策を見つけました。

これは以前にjuanが提案したtruncnormパッケージのdtruncnorm関数を使用して行うことができます。私は、平均5、標準偏差1、限界1と10の問題の例を使ってこれを実証します。

まず、分布からプロットする点の数のベクトルを作成する必要があります、それは等間隔になります。ですから、等間隔の分布から10点をプロットしたい場合たとえば、あなたのベクトル（「VEC」）は次のようになります。

vec=seq(from=1,by=1,length.out = 10) 

上記の我々は1ずつ増加、1から始まる10点をプロットすることを保証します、値に10

まで私たちは、その後dtruncnorm（）関数内でこれらを入れて、「テスト」変数に保存し、それをプロットします

test=dtruncnorm(vec,a=1,b=10,mean=5,sd=1) 
plot(test) 

見ての通り、これは密度をプロット離散点で。あなたは簡単にしようとすると、ポイントの数を増やすことによってプロットが連続作ることができます：

vec=seq(from=1,by=0.1,length.out = 100) 
test=dtruncnorm(vec,a=1,b=10,mean=5,sd=1) 
plot(test) 

ご質問

Answer 2

次の関数は、 '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)