2017-06-01から1日間の記事一覧

ディリクレ過程のメモ

Stick Breaking Process R package dirichletprocessがあるが、あえて自作すると: # Dirichlet process alpha <- 5 n <- 1000000 v <- rbeta(n,1,alpha) hist(v) ps <- rep(0,n) logv <- log(v) log1v <- log(1-v) cumsumlog1ps <- cumsum(log1v) ps[1] <- …