Challenge!!
The shortest Bayesian estimation in the world(?)
世界一短いベイズ推定への挑戦。
The shortest(2017/05/19 16:14)
2017/05/19 16:14 Daiki Hojo(@dastatis)
“Thanks to all contributors. Improper? Never mind!”
library(rstan)
stan_hist(stan(,,"parameters{real p;}"))
The shortest proper sampling
2017/05/19 16:10 Kentaro Matsuura(@berobero11) “It does not need the model{} block, in the recent version!”
library(rstan)
stan_hist(stan(,,"parameters{simplex[2] p;}"))
2017/05/19 14:20 Mr.Unadon(@MrUnadon) “Thanks to all contributors below. Now that this is the shortest stan code, I think.”
library(rstan)
stan_hist(stan(,,"parameters{simplex[2] p;}model{}"))
Excelent! graphics code!
2017/05/19 12:00 Daiki Hojo(@dastatis)
“In the case of focusing the graphic expression, this may be the shortest.”
library(rstan)
stan_hist(stan(,,"parameters{real p;}model{p~gamma(9,1);}"))
Excelent! sampling from the data!
2017/05/19 09:20 Hiroshi Shimizu(@simizu706) “In the case of estimation from data, this code will be the shortest.” (reported only stan code)
library(rstan)
#stan model code
"parameters{real p;}model{9 ~ poisson_log(p);}"
2017/05/19 09:00 Hiroshi Shimizu(@simizu706)
“This will be the shortest stan code. Warning? Never mind.” (reported only stan code)
library(rstan)
#stan model code
"parameters{real p;}model{p~beta(9,9);}"
2017/05/19 00:20 Mr.Unadon(@MrUnadon)
library(rstan)
hist(extract(stan(,,"parameters{real p;}model{p~gamma(9,1);}"))$p)
2017/05/19 00:00 Thanks to Mr.kos59125(He gave me an excellent advice), and revised by…
Mr.Unadon(@MrUnadon)library(rstan)
hist(extract(stan(,,"parameters{simplex[2]p;}model{1~bernoulli(p[1]);}"))$p[,1])
Excelent! Theoretical expression!
2017/05/18 22:40 Reported and Advised by…
kos59125(@kos59125)library(rstan)
hist(extract(stan(model_code="parameters{simplex[2]p;}model{1~bernoulli(p[1]);}"))$p[,1])
2017/05/18 21:40 First Report
Mr.Unadon(@MrUnadon)library(rstan)
hist(extract(stan(model_code="parameters{real<lower=0,upper=1>p;}model{int N;N=1;N~bernoulli(p);}"))$p)