暑いですね。夏ですね。 今年も広島ベイズ塾は夏合宿を行いました。そこで発表した，Stanコードの書き方中級編についてのスライドをアップしました。 内容としては，回帰分析から一般化線形モデル，欠損値のあるモデ … Continue reading →

誰得なのかよくわからない記事を書きます。※最後に追記あります 独立成分分析を調べていたら，ラプラス分布と正規分布の両方を含んだ「一般化正規分布」というのが出てきたので，調べてみました。 一般化正規分布は，平均μと，尺度パ … Continue reading →

  Twitterで質問があったので，書いておきます。 Stanでは，functions{}ブロックでユーザーが関数を定義できます。 普通の関数は，Rと同じように作ればいいのでそれほど難しくはないというか，違和 … Continue reading →

  先日，Stanの初心者講習の資料をアップしました。思った以上に反響があって驚いています。 さて，知り合いから媒介分析をStanでやるには？という質問があったので，コードを書いておきます。レベル的には初心者レ … Continue reading →

  先日関学で行われた，Harukara.Rで行われた，Stan超初心者講習会のスライドをスライドシェアにアップしました。 プログラミングを知らない人向けに書いているので，その辺りがわかっている人からすると，や … Continue reading →

  すでに関数自体は作っていたglmmstanですが，それをパッケージ化しました。 関数名と同じ，glmmstanパッケージです。   インストール glmmstanパッケージは，CRANに登録してい … Continue reading →

glmmstanを作ってから，MCMCで分析するのが楽しくて楽しくて仕方がない清水です。 で，glmmstanでWAICも計算できるのですが，どうもわからないことが出てきました。で，いろいろ考えた結果，「わからん」という … Continue reading →

  今回はRの話です。 社会心理学会の方法論セミナーでもGLMMをとりあげましたが，階層ベイズの話も久保先生のトークの中にありました。 GLMMでは，変量効果の数が増えると最尤法だと推定が難しくなったりするので … Continue reading →