{"published":"2018-09-22 13:17:23","description":"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\u2026","author_name":"ryamada22","html":"<iframe src=\"https://hatenablog-parts.com/embed?url=https%3A%2F%2Fryamada22.hatenablog.jp%2Fentry%2F20180922%2F1537589843\" title=\"\u30ce\u30f3\u30d1\u30e9\u30e1\u30c8\u30ea\u30c3\u30af\u30fb\u30d9\u30a4\u30ba\u5b9f\u8df5\u7de8 - ryamada\u306e\u907a\u4f1d\u5b66\u30fb\u907a\u4f1d\u7d71\u8a08\u5b66\u30e1\u30e2\" class=\"embed-card embed-blogcard\" scrolling=\"no\" frameborder=\"0\" style=\"display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;\"></iframe>","version":"1.0","height":"190","blog_url":"https://ryamada22.hatenablog.jp/","title":"\u30ce\u30f3\u30d1\u30e9\u30e1\u30c8\u30ea\u30c3\u30af\u30fb\u30d9\u30a4\u30ba\u5b9f\u8df5\u7de8","image_url":"https://cdn-ak.f.st-hatena.com/images/fotolife/r/ryamada22/20180922/20180922131657.jpg","url":"https://ryamada22.hatenablog.jp/entry/20180922/1537589843","provider_name":"Hatena Blog","categories":["\u30ce\u30f3\u30d1\u30e9\u30e1\u30c8\u30ea\u30c3\u30af","\u30d9\u30a4\u30ba","R","python"],"blog_title":"ryamada\u306e\u907a\u4f1d\u5b66\u30fb\u907a\u4f1d\u7d71\u8a08\u5b66\u30e1\u30e2","width":"100%","author_url":"https://blog.hatena.ne.jp/ryamada22/","provider_url":"https://hatena.blog","type":"rich"}