Dear Dieter, That sounds like a good solution, thanks.
Shige On Mon, Apr 12, 2010 at 2:27 AM, Dieter Menne <dieter.me...@menne-biomed.de> wrote: > > > Shige Song wrote: >> >> >> What is the preferred way to get Bayesian analysis results (such as >> those from MCMCpacki, MCMCglmm, and DPpackage) into LaTeX table >> automatically? >> > > It depends on what you want as output. Let's assume the summary of an > MCMCpack function: > First direct it to a variable, and check the available data frames. Then > output these with latex. > > If you do this very often: latex is a generic, so you could write a simple > function that creates customized output for a class "summary.mcmc". For > example, it could output both statistics and quantiles together, and do > reasonable rounding based on the standard errors. > > Dieter > > > library(MCMCpack) > library(Hmisc) > x<-rep(1:10,5) > y<-rnorm(50,mean=x) > qreg <- summary(MCMCquantreg(y~x)) > str(qreg) > #List of 6 > # $ statistics: num [1:3, 1:4] 0.3479 0.9029 0.3413 0.2602 0.0391 ... > # ..- attr(*, "dimnames")=List of 2 > # .. ..$ : chr [1:3] "(Intercept)" "x" "sigma" > # .. ..$ : chr [1:4] "Mean" "SD" "Naive SE" "Time-series SE" > # $ quantiles : num [1:3, 1:5] -0.139 0.825 0.257 0.168 0.878 ... > # ..- attr(*, "dimnames")=List of 2 > # .. ..$ : chr [1:3] "(Intercept)" "x" "sigma" > # .. ..$ : chr [1:5] "2.5%" "25%" "50%" "75%" ... > ... > # - attr(*, "class")= chr "summary.mcmc" > > latex(qreg$statistics,file="") > latex(qreg$quantiles,file="") > > -- > View this message in context: > http://n4.nabble.com/MCMC-results-into-LaTeX-tp1836393p1836723.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.