混合ディストリビューションからサンプルを生成する必要があります
40% のサンプルは Gaussian (平均 = 2、標準偏差 = 8) から取得
20% のサンプルは Cauchy (場所 = 25、スケール = 2) からのものです。
40% のサンプルは Gaussian (平均 = 10、sd=6) から取得
これを行うために、私は次の関数を書きました:
dmix <- function(x){
prob <- (0.4 * dnorm(x,mean=2,sd=8)) + (0.2 * dcauchy(x,location=25,scale=2)) + (0.4 * dnorm(x,mean=10,sd=6))
return (prob)
}
そして、次のようにテストしました:
foo = seq(-5,5,by = 0.01)
vector = NULL
for (i in 1:1000){
vector[i] <- dmix(foo[i])
}
hist(vector)
このようなヒストグラムを取得しています(これは間違っていることがわかっています)-
私は何を間違っていますか?誰でもいくつかの指針を教えてもらえますか?