similar to: Two questions about R2BayesX package

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 in the bayesx help: ## generate some data

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*x2)/3+2 * sin(4 * pi * (x3

Hi, this is the first announcement of a *developer version* of the R package BRugs, which is based on OpenBUGS 2.1.0. OpenBUGS 2.1.0 has been released earlier this week (http://mathstat.helsinki.fi/openbugs/). BRugs contains a complete OpenBUGS installation as well as R code to access OpenBUGS functionality from R. Simply install the package, load it and say ?BRugs and ?BRugsFit

Dear All, I am trying to use the glmmBUGS package with OpenBUGS under Linux. Here is a short list of code I modified from the diseasemapping package document: ---------------------------------------- library(diseasemapping) data(popdata) data(casedata) model = getRates(casedata, popdata, ~age*sex, breaks=seq(0, 90, by=10) ) ontario = getSMR(popdata,model, casedata) spplot(ontario,

This is slightly off-topic, but there was a discussion about MCMC a couple of weeks ago. OpenBUGS can operate from R, at least in Windows (there are some problems in Linux at the moment). The interface with R is one are that there are plans to work on: it's a bit basic at the moment. Bob -- Bob O'Hara Dept. of Mathematics and Statistics P.O. Box 68 (Gustaf H??llstr??min katu 2b)

Dear R-users, I am trying to use the following code to reproduce results from Prof. Gelman's book, but have the listed error for R2WinBUGS version (the openbugs version is good). I am using R-2.6.1 on windows XP, and all the R packages are most current ones. schools.bug can be found at http://www.stat.columbia.edu/~gelman/bugsR/runningbugs.html . Can anyone help me to figure out what's

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 While I'm not an R expert, I have used R on Windows XP. Now I've moved to Ubuntu (Intrepid), and I'm trying to configure R to work with the Gelman and Hill _Data Analysis Using Regression and Multilevel/Hierarchical Models_. So far, it's not working. I start by following the instructions for installing arm and BRugs at

Greetings: I'm going to encourage some students to try Bayesian ideas for hierarchical models. I want to run the WinBUGS and R examples in Tony Lancaster's An Introduction to Modern Bayesian Econometrics. That features MS Windows and "bugs" from R2WinBUGS. Today, I want to ask how people are doing this in Linux? I have found a plethora of possibilities, some of which are not

Dear list, I'm trying to embed OpenBUGS in an R package for use of it on 64-bit Linux. In order to get the CLI working one has to compile C code contained in $OpenBUGS/Manuals/CBugs.html (copied to say CBugs.c) using gcc -m32 -o bugs CBugs.c -ldl I put the OpenBUGS distribution in the ./inst subdirectory of the package root. Where should I now put the CBugs.c and how

Dear All, I am trying to graph a simple scatter plot where the x axis is year and the y axis is a percentage (percentage of infant death). Instead of plotting the raw data, I want to plot summary statistics such as mean and median. Here is the problem: the value range of y is between 0 and 1, but since infant death is a rare event, the mean and median is very low (something like 5%), which shows

Dear All, After playing with lmer for couple of days, I have to say that I am amazed! I've been using quite some multilevel/mixed modeling packages, lme4 is a strong candidate for the overall winner, especially for multilevel generzlized linear models. Now go back to my two-level poisson model with cross-classified model. I've been testing various different model specificatios for the

Dear All, Are there R packages that can estimate survival model with long-term survivors? This is sometimes known as "cure" model or "split-population" model. Thanks. Shige [[alternative HTML version deleted]]

Actually, I re-read the post and think it needs clarification. We may both be right. If the question is "I am building a model and want to know if I should retain this random effect?" (or something like that) then the LRT should be used to compare the fitted model against another model. This would be accomplished via anova(). In other multilevel programs, the variance components are

Dear All, When I tried to run the following code (taken from the *spdplyr* package vignettes): library(spdplyr) library(maptools) data(wrld_simpl) worldcorner <- wrld_simpl %>% mutate(lon = coordinates(wrld_simpl)[,1], lat = coordinates(wrld_simpl)[,2]) %>% filter(lat < -20, lon > 60) %>% dplyr::select(NAME) I got the following error messages: Error in (function (cl,

Dear All, Can anyone give me some hints about how to set starting values for a lmer model? For complicated models like this, good starting values can help the numerical computation and make the model converge faster. Thanks! Shige [[alternative HTML version deleted]]

Dear Colleagues, In the following simple ggplot2 code: -------------------------------- m <- ggplot(d.fig, aes(time, prob)) m + stat_summary(fun.data = "median_hilow", conf.int = .95, geom = "smooth") + facet_wrap(~ Cohort, nrow=1) + coord_cartesian(ylim = c(0, .03)) -------------------------------- Is there a way to replace the y-axis label from "prob" to

Hi all, I used both OpenBugs and R function bugs{R2WinBUGS} to run a linear mixed effects model based on the same data set and initial values. I got the same summary statistics but different posterior samples. However, if I order these two sets of samples, one is generated from OpenBugs and the other is generated from R, they turn to be the same. And the samples from R do not have any

Hi, I'm working on the codes below however every time I run them when they get to OpenBUGS I keep getting the error message: array index is greater than array upper bound for hab. Any help would be greatly appreciated, Suzie Codes: ungulate <- read.csv(file.choose ()) #ungulate ungulate <- as.matrix(ungulate);colnames(ungulate)<-NULL;rownames(ungulate)<-NULL

Dear All, I am trying to compile&install the package "RGtk2" on my Ubuntu 10.04 box. I did not have problem with earlier versions, but with the new version, I got the following error message : ------------------------------------------------------------------------------------------------- * installing *source* package ?RGtk2? ... checking for pkg-config... /usr/bin/pkg-config