廿TT

譬如水怙牛過窓櫺 頭角四蹄都過了 因甚麼尾巴過不得

Stan

RStan でポアソン過程の停止時刻を推定する

問題設定 n = 20 人の患者さんが繰り返し訪問する病院の窓口を T = 100 時間観察した。患者 i はそれぞれある時期 に達すると病院への訪問をやめる(理由は不明:病気が治ったか、遠くに引っ越したか、通院に飽きたか、死んだか)。この訪問がレート λ (全…

RStanで自殺者数推移の変化点を調べる

日本の自殺者数は平成10年に急増して3万人を超えて以来高い水準で推移していましたが、平成24年に久しぶりに3万人を下回り、今のところ減少傾向です。このデータに変曲点のある回帰直線、を当てはめて変化点をうまく検出できるか試してみます。コードは [RS…

[RStan]変曲点のある単回帰分析

豊田秀樹『マルコフ連鎖モンテカルロ法』5章8節の例題。マルコフ連鎖モンテカルロ法 (統計ライブラリー)作者: 豊田秀樹出版社/メーカー: 朝倉書店発売日: 2008/05メディア: 単行本購入: 11人 クリック: 168回この商品を含むブログ (13件) を見る変曲点のある…

久保(2012)ベイズモデル入門のための例題

概要 データ解析のための統計モデリング入門――一般化線形モデル・階層ベイズモデル・MCMC (確率と情報の科学)作者: 久保拓弥出版社/メーカー: 岩波書店発売日: 2012/05/19メディア: 単行本購入: 16人 クリック: 163回この商品を含むブログ (23件) を見る生態…

ワイブル分布のパラメータ推定(Stan vs survreg)

※初公開時は Stan のコードが間違ってたせいでパラメータがうまく求まっていなかった。(修正:7/28)ワイブル分布のパラメータ推定は意外とむずかしくって、打ち切りが多いときとかパラメータの数が多いときはシェイプパラメータが変なところに飛んでっちゃ…