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をリスト形式にま…