2021-09-01から1ヶ月間の記事一覧

Metropolis Hastings サンプリングする

# 値の台 xx <- seq(from=-5,to=20,length=1000) # 適当な(密度関数) dd <- dnorm(xx) + dnorm(xx,1,2) + dnorm(xx,10,2) # 確率密度関数っぽい形 plot(xx,dd) # library(mcmc) # 確率密度に比例したあたいの対数を返す関数を書く h <- function(x){ log(dno…