search for: r2bayesx

Hi all I'm using the function bayesx to estimate a simple model, for example: library(R2BayesX) ## generate some data set.seed(111) n <- 200 ## regressor dat <- data.frame(x = runif(n, -3, 3)) ## response dat$y <- with(dat, 1.5 + sin(x) + rnorm(n, sd = 0.6)) ## estimate models with ## bayesx REML and MCMC b1 <- bayesx(y ~ sx(x), method = "REML", data = dat) I want...

Hi all, I have a problem with the library R2BayesX, when i try to use the command bayesx i get this error: dyld: Library not loaded: /usr/local/lib/libreadline.5.2.dylib Referenced from: /Library/Frameworks/R.framework/Versions/2.15/Resources/library/BayesXsrc/libs/i386/BayesX Reason: image not found I obtain this message also with the example...

Dear All, I have two questions regarding the use of the R2BayesX package for Bayesian analysis. First, is it possible to generate predictions based on the fitted model? According to Gelman and Hill (2007, pp. 361-363), there are at least two ways to do this in BUGS: (1) generate additional data points with the dependent variable coded as missing (and all the in...

Hi, I wonder if anyone can help me with this issue. I am using R2BayesX. It seems that the model can maximally contain 20 interactions. When the number of interaction terms exceed 20, the code stops working. Here is a piece of toy code. rm(list=ls()) library(BayesX) library(R2BayesX) #data generating model f2<-function(x1,x2,x3,x4) { y<-2*sin(pi*x1)*1.5+exp(2...

Hi all, I have a set of 54 files that I need to convert from ASCII grid format to .shp files to .bnd files for BayesX. I have the following R code to operate on those files: library(maptools) library(Grid2Polygons) library(BayesX) library(BayesXsrc) library(R2BayesX) readfunct <- function(x) { u <- readAsciiGrid(x) } modfilesmore <- paste0("MaxFloodDepth_", 1:54, ".txt") modeldepthsmore <- lapply(modfilesmore, readfunct) maxdepth.plys <- lapply(modeldepthsmore, Grid2Polygons(modeldepthsmore, level = FALSE)) layers <- p...