Bayesian Computation with R: Second Edition (Use R!) 2値反応回帰モデルの一種であるプロビットモデルの回帰係数の同時事後分布からのサンプリング 回帰係数の事前分布は精度行列がゼロ行列の場合とそうでない場合の両方に対応する。 後者の場合は対数周…

Bayesian Computation with R: Second Edition (Use R!)離散確率分布に対して最高確率区間を計算する。 最高確率区間：こちらから与える確率を含む最短区間 確率分布が単峰でない場合は，2つ以上の区間に分かれる。 確率分布が左右対称でない場合は，区間の…

Bayesian Computation with R: Second Edition (Use R!)（ヒストグラムのように）均等な区間ごとに確率が定義された確率分布の密度を計算する使用法 histprior(p,midpts,prob) 引数 p:連続型確率変数において，対応する密度を計算する確率変数の値をまとめた…

Bayesian Computation with R: Second Edition (Use R!)百分位点と対応する事前累積確率に関する2組の情報に基づき，それにマッチするベータ分布のパラメータを求める。使用法 beta.select(quantile1, quantile2) 引数 quantile1 ある100分位点x1（例 50%点…

Bayesian Computation with R: Second Edition (Use R!)比率に関する離散分布を用いた場合の将来の二項分布による実験の成功数に関する予測分布の計算使用法 pdiscp(p, probs, n, s) 引数 p:比率の候補となる値で構成されるベクトル probs:各比率の確率で構…

Bayesian Computation with R: Second Edition (Use R!)離散事前分布を用いた場合の比率に関する事後分布の計算使用法 pdisc(p, prior, data) 引数 p:（推測する）比率の候補となる値で構成されるベクトル prior:候補となる各比率の事前確率で構成されるベク…

Bayesian Computation with R: Second Edition (Use R!) データ分布 y~N(mu,sig2)事前分布 p(mu,sig2)proportional to (sig2)^(-1)とする。このときの(mu,sig2)の事後分布（カーネル）の対数に関するある(mu_0,sig2_0)における値。 使用法 normchi2post(thet…

Bayesian Computation with R: Second Edition (Use R!) Zellnerのg-priorbeta (given sigma2) ~N(b0,c0*sigma2*(X'X)^(-1)) p(sigma2)proportional to (sigma2)^(-1)の場合の，回帰モデルの事後分布のカーネルの対数を計算。注 この関数の出力が関数laplace…

Bayesian Computation with R: Second Edition (Use R!)作者: Jim Albert出版社/メーカー: Springer発売日: 2009/05/15メディア: ペーパーバック購入: 1人 クリック: 12回この商品を含むブログ (20件) を見る無情報事前分布を用いた線形回帰モデルに関する（…

Bayesian Computation with R: Second Edition (Use R!)作者: Jim Albert出版社/メーカー: Springer発売日: 2009/05/15メディア: ペーパーバック購入: 1人 クリック: 12回この商品を含むブログ (20件) を見る無情報事前分布を用いた線形回帰モデルの （被説…

Bayesian Computation with R: Second Edition (Use R!)無情報事前分布を用いた線形回帰モデルの（被説明変数の）期待値に関する事後分布からのサンプリング使用法 blinregexpected(X1,theta.sample) 引数 X1：行列（行が説明変数の組に対応） 推定に用いな…

Bayesian Computation with R: Second Edition (Use R!) 正規分布のパラメータ（mu,sigma2）の無情報事前分布を前提にして そのパラメータの事後分布からシミュレーションした上で その事後分布に基づき，関心のある統計量の事後予測分布を求める。使用法 no…

Bayesian Computation with R: Second Edition (Use R!) データ分布y~N(Xbeta,sigma2*I)無情報事前分布p(bea,sigma2)proportional to (sigma2)^(-1)または Zellnerのg-priorbeta (given sigma2) ~N(b0,c0*sigma2*(X'X)^(-1)) p(sigma2)proportional to (sigm…

Bayesian Computation with R: Second Edition (Use R!)データ分布 y~N(mu,sigma2) 事前分布 p(mu,sigma2)proportional to (sigma2)^(-1) の場合の(mu,sigma2)の事後分布からのサンプリング。 この場合，解析的な結果を得ることができるので， サンプリング…

Bayesian Computation with R: Second Edition (Use R!)指定した整数値ごとに行列の行を集計して新しい行列を作成する。 使用法 regroup(data,g) 引数 data:行列 g:1からdataの行数nの間の正整数値 gごとに行列dataの行を集計した新しい行列 行列dataの行数…

Bayesian Computation with R: Second Edition (Use R!)グリッド上で定義された一般的2変量密度関数より乱数を発生させる使用法 simcontour(logf,limits,data,m) 引数 logf:適用する対数密度を定義した関数 limits:グラフを描く領域c(xlo,xhi,ylo,yhi)を指定…

Bayesian Computation with R: Second Edition (Use R!) 2変量密度関数の等高線の描画一般的な2変量密度関数用の等高線グラフ 実際には，対数密度関数を使ってグラフを描く。 密度関数の定数部分は除いても描ける。 モード（最頻値）での密度関数の高さを基…

Bayesian Computation with R: Second Edition (Use R!) 2変量正規分布の密度関数の自然対数の計算 定数部分を除くカーネルの自然対数 使用法 lbinorm(xy,par)引数 xy: 密度関数の対数を評価する点(x,y) par:平均ベクトルmと分散共分散行列vをリスト形式にま…

Bayesian Computation with R: Second Edition (Use R!)切断された確率分布からのサンプリング累積分布関数と（分位数を出力する）逆累積分布関数（分位数） が利用可能な場合に適用できる。 Rの場合，例えば正規分布なら 累積分布関数はpnorm 分位数を出力…

Bayesian Computation with R: Second Edition (Use R!)逆ガンマ分布IG(a,b)からのサンプリング ガンマ分布G(a,b)に従う確率変数Xの逆数Y=1/Xは （定義上）逆ガンマ分布に従うという事実を直接用いる。使用法 rigamma(n, a, b) 引数 n: 発生させる乱数の数 a…

Bayesian Computation with R: Second Edition (Use R!) 多変量t分布からのサンプリング 使用法 rmt(n = 1, mean = rep(0, d), S, df = Inf) #正値定符号行列Sを必ず指定しなければならない。 #自由度dfはデフォールトでは無限大Inf # dfを指定しなければ，…

Bayesian Computation with R: Second Edition (Use R!)-多変量正規分布からのサンプリング使用法 rmnorm(n = 1, mean = rep(0, d), varcov) #分散共分散行列varcovは必ず入力。 引数 n サンプリングする回数 mean 平均ベクトル varcov 分散共分散行列（正値…