二番目に簡単な rstan コード - Qiita
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); alpha ~ normal(0, 100); beta ~ normal(0, 100); s ~ inv_gamma(0.001, 0.001); } ' datastan <- list(N=N, x=x, y=y) fit <- s
Stan で相関係数を推定する - Qiita
library(rstan) N <- 1000 x <- rnorm(N, mean = 50, sd = 10) y <- 10 + 0.8 * x + rnorm(N, mean =0, sd = 7) data <- data.frame(x, y) stancode <- ' data{ int<lower=0> N; vector[2] y[N]; } parameters { vector[2] mu; real<lower=0> var_x; real<lower=0> var_y; real cov; } transformed parameters { matrix[2,2] sigma; sigma[1,1] <- var_x; sigma[2,1] <- cov; sigma[1,2] <- cov; sigma[2,2] <- var_y; } model{ fo
Rstanのインストール - 驚異のアニヲタ社会復帰の予備
巷で話題のRstanをインストールだけしておいた。 OSはubuntu 12.04 LTS. なんか困ったらhttp://heartruptcy.blog.fc2.com/かAnalyze IT.を見れば問題ない。 # 依存関係 install.packages("inline") install.packages("Rcpp") install.packages("coda") # あると便利 # Rstan options(repos = c(getOption("repos"), rstan = "http://wiki.rstan-repo.googlecode.com/git/")) install.packages('rstan', type = 'source') Eight Schoolsの例をとりあえずやってみる。計算結果のstanfitクラスオブジェクトのデフォルトプロ
1
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
処理を実行中です
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く
Stan で increment_log_prob() - Qiita
stanで階層ベイズ（空間構造を入れる） | @DataSci
