ryamada22 https://blog.hatena.ne.jp/ryamada22/ ryamadaの遺伝学・遺伝統計学メモ https://ryamada22.hatenablog.jp/ ノンパラメトリック ベイズ R python Non-parametric bayesian clustering Data set simulation n <- 1000 d <- 4 k <- 5 s <- sample(1:k,n,replace=TRUE) m <- matrix(rnorm(d*k,sd=6),k,d) X <- matrix(0,n,d) for(i in 1:k){ tmp <- which(s == i) r <- matrix(rnorm(d * length(tmp))) X[tmp,] <- r for(j in 1:d){ X[tmp,j] <- X[tmp,j] + m[i,j] } } plo… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fryamada22.hatenablog.jp%2Fentry%2F20180922%2F1537589843" title="ノンパラメトリック・ベイズ実践編 - ryamadaの遺伝学・遺伝統計学メモ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.f.st-hatena.com/images/fotolife/r/ryamada22/20180922/20180922131657.jpg Hatena Blog https://hatena.blog 2018-09-22 13:17:23 rich https://ryamada22.hatenablog.jp/entry/20180922/1537589843 1.0 100%