1. hoxo_m

    No comment

    hoxo_m
Changes in body
Source | HTML | Preview
@@ -1,47 +1,51 @@
Kosugitti 先生の[世界一簡単なrstanコード](http://kosugitti.net/archives/4823)に触発されて、
二番目に簡単な rstan コードを書いてみました。
```r
library(rstan)
N <- 1000
x <- rnorm(N, mean = 50, sd = 10)
y <- 10 + 0.8 * x + rnorm(N, mean =0, sd = 7)
stancode <- '
data{
int<lower=0> N;
real x[N];
real y[N];
}
parameters {
real alpha;
real beta;
real<lower=0> s;
}
model{
for(i in 1:N)
y[i] ~ normal(alpha + beta * x[i], s);
s ~ inv_gamma(0.001, 0.001);
}
'
datastan <- list(N=N, x=x, y=y)
fit <- stan(model_code=stancode, data=datastan, iter=1000, chain=4)
traceplot(fit, ask=T)
print(fit, digit=2)
```
![Rplot01.png](https://qiita-image-store.s3.amazonaws.com/0/2570/dfec861b-f5d9-8ab0-e714-b7db8480cb0b.png)
```:結果
mean se_mean sd 2.5% 25% 50% 75% 97.5% n_eff Rhat
alpha 11.60 0.06 1.11 9.29 10.90 11.60 12.33 13.73 396 1.01
beta 0.77 0.00 0.02 0.73 0.76 0.77 0.79 0.82 404 1.01
s 7.01 0.01 0.15 6.71 6.90 7.01 7.10 7.32 438 1.00
lp__ -2447.58 0.06 1.26 -2450.88 -2448.09 -2447.25 -2446.70 -2446.22 451 1.00
```
alpha の推定がうまくいってないような気がするんですが。。。
ご査収ください。
+
+関連記事:
+[Stan で相関係数を推定する](http://qiita.com/hoxo_m/items/0f1b05681f5d6c4b560a)
+[Stan で increment_log_prob()](http://qiita.com/hoxo_m/items/56a63faf29e1af4db27d)