密度 f(x) = 0.6φ(x)+0.4ψ(x) † φを平均-1,分散1の正規分布に従う確率変数,ψを平均2,分散1の正規分布に従う確率変数として f(x) = 0.6φ(x)+0.4ψ(x) となる密度関数 truedensity() を定義する. truedensity <- function (x) { 0.6/sqrt(2*pi)*exp(-(x+1)^2/2) + 0.4/sqrt(2*pi)*exp(-(x-2)^2/2) } > curve(truedensity, xlim=c(-6,6), ylim=c(0,0.3), col=2) f(x) に従う乱数を生成する関数 † 次に, f(x) に従う乱数を生成する関数を定義する. generator_tmp <- function(n) { data1 <- rnorm(n)-1 # φ(x) に従う乱数 data2 <