裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

混合分布

2010年07月15日 | ブログラミング
ある平均値と標準偏差を持つ2つの正規分布の混合分布の乱数を返す関数を定義する
gen <-function(m1=0, sd1=1, m2=1, sd2=1, p=0.5) # p は,1番目の分布に従うサンプルの割合
{
    return(ifelse(runif(1) < p, rnorm(1, m1, sd1), rnorm(1, m2, sd2)))

}
n <- 100000
x <- numeric(n)
for (i in 1:n) x[i] <- gen(m2=3, sd2=0.5, p=0.9)
hist(x,breaks=50, freq=FALSE, main="")

コメント    この記事についてブログを書く
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする
« まあ,893(やくざ)ですな | トップ | 2直線の交点座標を求めるプロ... »
最新の画像もっと見る

コメントを投稿

ブログラミング」カテゴリの最新記事