tests/mcmc.R

## BayesX MCMC testing
library("BayesXsrc")
mcmc <- run.bayesx("mcmc.prg", verbose = FALSE)
fx1 <- read.table("mcmc_f_x1_pspline.res", header = TRUE)
fx2 <- read.table("mcmc_f_x2_pspline.res", header = TRUE)
print(round(head(fx1), digits = 3))
print(round(head(fx2), digits = 3))

Try the BayesXsrc package in your browser

Any scripts or data that you put into this service are public.

BayesXsrc documentation built on Feb. 16, 2023, 6:56 p.m.