MCMC診断をLatexテーブルに変換するには？

Question

codaパッケージを使用して、私のMCMCの要約統計量を計算しています。しかし、印刷された要約をラテックステーブルに変換するオプションはないようです。私はstargazerを試して、summary.mcmcの結果をデータフレームに強制しました。どちらの試みも失敗しました。

library(coda) 
mock_mcmc <- mcmc(rnorm(1000)) 
summary(mock_mcmc) 

summary.mcmcが

1. Empirical mean and standard deviation for each variable, 
    plus standard error of the mean: 

      Mean    SD  Naive SE Time-series SE 
     0.03180  0.98715  0.03122  0.03368 

2. Quantiles for each variable: 

    2.5%  25%  50%  75% 97.5% 
-1.89794 -0.65289 0.02952 0.67396 1.97158 

プリントアウトされますどのようにLaTeXファイルにテーブルを結果I出力を行います。

ここで再現可能な例ですか？要約統計量を手で計算することは可能だと私は理解していますが、分かりにくいコーダの便利な機能があるかどうかは不思議です。

Answer 1

コーダでそれを行う方法があるかどうかわかりません。以下のメソッドは、オブジェクトをデータフレームに変換してからxtableを再実行することによって、mcmcオブジェクトのxtableのメソッドを定義します。

Answer 2

重要なのは、出力から関連する量を抽出することです。適切な量​​を見つけるには、まさにあなたがこの表にしたいですかどうR.

library(coda) 
library(xtable) 

mock_mcmc = mcmc(rnorm(1000)) 
s = summary(mock_mcmc) 

stats_table = xtable(as.data.frame(t(s$statistics))) 
quant_table = xtable(as.data.frame(t(s$quantiles))) 

print(stats_table, file="stats_table.tex") 
print(quant_table, file="quant_table.tex")