記事 での「MCMC」の検索結果 86件
MCMC再訪(1)
秋の某統計高座のテキスト作成のため、もう一度MCMCを最初からやってみる。
例題1:
正規分布することがわかっている ある母集団から、X = (5.89, 5.69, 4.71, 6...
タグ: MCMC winbugs
Zero-inflated Poisson
ZIP (Zero-inflated Poisson) モデルのベイズ推定。Martin et al. (2005)を参考にしてみる。R2WinBUGSを使用。
タグ: winbugs zero-inflated poisson MCMC
負の二項分布のベイズ推定
負の二項分布のパラメータをJAGSで推定するテスト。
データは
n <- 100
p <- 0.8
r <- 10
num <- rnbinom(n, r, ..
タグ: jags 負の二項分布 MCMC
MCMC (4)
R News 6/1はMCMC特集。MCMCpackとBUGSの使い方を解説。
タグ: MCMC
MCMC (3)
生態学会おわる。自由集会や統計合宿もあったので、MCMCの3回目。帰りの列車の中から、混合モデルをいろいろいじってみる。が、うまくいかない。
タグ: MCMC
多重共線性とか潜在変数とか共分散構造分析とか
多重共線性のことを考えることがあって、モデル選択では解決できなかったが、潜在変数を導入してみれば、と、そうすると共分散構造分析(SEM)になるなぁ、とそういえば、
豊田(2008)
ではS..
タグ: R SEM MCMC
R: glmmAKによる累積ロジットモデルのベイズ推定
glmmAKで累積ロジットモデルのベイズ推定をしてみる。累積ロジットモデルのベイズ推定(2)で使ったデータを使用。時間によるランダム効果は無視した。
# inverse logit
invlo..
タグ: 累積ロジットモデル glmmAK MCMC R
MCMC (2)
MCMCの2回目。線形回帰と、一般化線形モデル(相当のこと)をMCMCでやってみる。
タグ: MCMC
MCMCの勉強(1)
今さら感はあるが、MCMC (Markov Chain Monte Carlo; マルコフ連鎖モンテカルロ)を使えるようになろうと、まずは簡単な例から試してみた。
タグ: MCMC
累積ロジットモデルのベイズ推定
以前vglmで解析したデータをベイズ推定で解析してみた。nプロットについて、ある植物の被度(c1)を5段階で計測するというような場合を想定。被度は傾斜(slope)が急になるほど減るとする。
..
タグ: MCMC 累積ロジットモデル
事後確率分布に従う乱数
森林学会のポスター発表で話していて、MCMC推定した事後確率分布に従うような乱数を得られないか、という話になった。で、MCMCサンプルから値を得てくればどうか、というような結論になったのだが、うまくい..
タグ: MCMC R