similar to: getting started on Bayesian analysis

Hi Does anyone know of any R coding/functions for MCMC approaches? I am currently using BUGS but I wonder if the bazaar has produced anything? I think I am pushing BUGS to it's limit and possibly past it at the moment. John -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Dear R-helpers, Hi! All. I'm doing a project which needs MCMC simulation. I wonder whether there exists related packages in R. The only one I know is a MCMCpack package. What I want to do is implementing gibbs sampling and Metropolis-Hastings Algorithm to get the posterior of hierarchical bayesian models. Thanks in advance. Jun

Is there a way of running BUGS on OS X (from R)? I only see Windows versions on their website. If not, what are the alternatives for Bayesian analysis? Thanks, Tamas

The document R-exts contains the following example of using mathematics in R documentation. \deqn{p(x) = {\lambda^x\ \frac{e^{-\lambda}}{x!}} {p(x) = lambda^x exp(-lambda)/x!} There is a syntax error in there, but that's not my point. The problem is that using "R CMD Rd2dvi" I find that putting the alternate forms of the equation on top of each other doesn't work.

Hi all, I would like to start to use R's MCMC abilities to compute answers in Bayesian statistics. I don't have any specific problems in mind yet, but I would like to be able to compute/sample posterior probabilities for low-dimensional custom models, as well as handle "standard" Bayesian cases like linear regression and hierarchical models. R clearly has a lot of abilities in

I have been running a Gibbs Sampler to estimate levels of efficiency in the Louisiana Shrimp Industry. I created a matrix (samp) where I stored the results of each iteration for 86 variables. I run 10,000 iterations. So, the matrix samp is 10,000 x 86. I want to use the gelman-rubin test to check for convergence. To do that, I need at least two chains. If I run second chain with different starting

Dear R-users (and developers), I am looking for an efficient framework to carry out parameter estimations based on MCMC (optionally with specified priors). My goal is as follow: * take ANY R-function returning a likelihood-value (this function may itself call external programmes or other code!) * run a sampler that covers the multidimensional parameter space (thus creating a posterior

Hi, The department of economics at our university (Budapest) is planning a course on numerical methods in economics. They are trying to decide which software to use for that, and I would like to advocate R. The other alternative is Matlab. I have found comparisons in terms of computational time for matrix algebra, but I don't think that is relevant: the bottleneck for economists is usually

Dear all, I am writing an R code to fit a Bayesian mixed logit (BML) via MCMC / MH algorithms following Train (2009, ch. 12). Unfortunately, after many draws the covariance matrix of the correlated random parameters tend to become a matrix with almost perfect correlation, so I think there is a bug in the code I wrote but I do not seem to be able to find it.. dull I know. Has anybody written a

Hi all. Are there any R functions around that do quick logistic regression with a Gaussian prior distribution on the coefficients? I just want posterior mode, not MCMC. (I'm using it as a step within an iterative imputation algorithm.) This isn't hard to do: each step of a glm iteration simply linearizes the derivative of the log-likelihood, and, at this point, essentially no

Dear all, I am trying to use the logistic regression with MCMClogit (package: MCMCpack/Coda) and I want to put a beta prior on the parameters, but it's giving me error message (please see output below) no matter what shape 1 or 2 I use. It works perfect with the cauchy or normal priors. Do you know if there is a catch there somewhere? Thanks logpriorfun <- function(beta,shape1,shape2){

Hello, Is there a function in R by which one can run AR model with Bayesian Information Criteria (BIC)? To my knowledge, functions ar and ar.ols could select the order only by AIC. Thanks, Miao [[alternative HTML version deleted]]

Hi Friends, I'm trying to model the consumer decisions (Click-Through Rate and Conversion) in Search Engine Advertising using a hierarchical Bayesian binary logit. The input data is the weekly CTRs and Avg. Position for each search keyword. CTR is modeled as (for each keyword i and week j): Pij = exp(C + Bi x Positionij + A1 x Lengthi + A2 x Brandi + A3 x ProductSpecifici) / [1 + exp(C +

I have been checking my package coda with the current R-devel using "R CMD check". As usual, this has uncovered a large number of errors and inconsistencies in my documentation, for which I am very grateful. The only problem I have is with code/documentation mismatches when I have written a method for a generic function, e.g. * checking for code/documentation mismatches ... WARNING

Dear R developers, these days I'm working on some R code for fitting completely generic Bayesian Hierarchical Models in R, a la OpenBUGS and JAGS. A key feature of OpenBUGS and JAGS is that they automatically build an appropriate MCMC sampler from a generic model, specified as a directed acyclic graph (DAG). The spirit of my (would-be) implementation is instead more focused on experimentation

Hi, people. I was a bit intrigued by the message quoted below. Indeed, if pt() is given a matrix, it returns a matrix. Should this feature be documented? ?pt speaks about "a vector of quantiles", and says nothing about the type of what it returns. The same might presumably apply to other distribution-related functions. ----- Forwarded message from Martyn Plummer <plummer at

This message is in MIME format --_=XFMail.1.1.p0.Linux:971128122615:3052=_ Content-Type: text/plain; charset=us-ascii I have rounded up three buglets in R-0.50-a4. Two of them I can fix and a patch is supplied below. I hope this is useful for the current source (if these haven't been fixed already :) 1) cov cov() fails when it's argument is a matrix with one column and with column names

Dear R users: I am using the package coda (the last verison in CRAN) to analyse the output from a MCMC Bayesian analysis. And I get unconsitented results. I have export the chain using the read.table function and after I have transformed this data frame to an mcmc object using the mcmc function. I am interested in three variables, when I use the function effectiveSize I have these figures:

[Sent to R-help by mistake. Sorry] I have been checking my package coda with the current R-devel using "R CMD check". As usual, this has uncovered a large number of errors and inconsistencies in my documentation, for which I am very grateful. The only problem I have is with code/documentation mismatches when I have written a method for a generic function, e.g. * checking for

In the following code I have created the posterior density for a Bayesian survival model with four parameters. However, when I try to use the adapt function to perform integration in four dimensions (on my old version of R I get an error message saying that I have applied a non-function, although the function does work when I type kernel2(param0, theta0), or on the newer version of R the computer