前回cmdstanpyを使って単回帰を解説したので、今回は重回帰分析用のコードをまとめてみる。 使用するデータ x1(説明変数) 1 2 3 4 5 x2(説明変数) 2 6 6 9 6 x3(目的変数) 1 3 4 7 9 過去の計算によると下記式が得られていた。 pythonコード 使用したpython…

stanで重回帰分析をしようとすると、こんな感じのモデルを書くことになる。 data { int<lower=0> N; //データ数 int<lower=0> M; //説明変数の数 real X[N, M]; real Y[N]; } parameters { vector[M] b; real<lower=0> sigma; } model { Y ~ normal (X*b, sigma); } このXbはちゃんと書</lower=0></lower=0></lower=0>…