>>462
stanで検証

functions{
real jisaku_log(real y){
real temp;
temp = 2/3.0*y+2/3.0;
return log(temp);
}
}
data{
}
parameters{
real<lower=0,upper=1> p;
}
transformed parameters{
real q;
q = p/(p+1);
}
model{
p ~ jisaku();
}


mean se_mean sd 2.5% 25% 50% 75% 97.5% n_eff Rhat
p 0.5626 0.0074 0.2836 0.0366 0.3364 0.5874 0.8107 0.9839 1464 1.0006
q 0.3364 0.0035 0.1325 0.0353 0.2517 0.3701 0.4477 0.4959 1418 1.0007
lp__ -2.0063 0.0336 0.9343 -4.5941 -2.2177 -1.6547 -1.4262 -1.3608 773 1.0030