Question

我使用coda包来计算我的MCMC的摘要统计信息。但是，似乎没有选项将打印的摘要转换为Latex表。我已尝试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文件中？我了解可以手动计算汇总统计数据，但我很好奇是否有一个我不了解的尾声的便捷功能。

Answer 1

我不知道是否有办法在尾声中做到这一点。下面的方法通过将对象转换为数据帧然后重新运行xtable来定义mcmc对象的xtable方法。

library(coda)
library(dplyr)
library(magrittr)
library(xtable)

xtable.summary.mcmc = function(x, ...)
  x %>%
  use_series(statistics) %>%
  c(x %>%
      use_series(quantiles) ) %>%
  as.list %>%
  dplyr::as_data_frame() %>%
  xtable(...)

1000 %>% 
  rnorm %>% 
  mcmc %>% 
  summary %>%
  xtable

如何将MCMC诊断转换为Latex表？

1 个答案: