similar to: MCMC revisited and BUGS ...

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

I see what Dean means about how Digium/Asterisk might have struck a balance between "the cathedral and the bazaar" antipodes of the SW development world. Nicholas Carr's "The Ignorance of Crowds" finally states his "politics" when it says "When you move from the bazaar to the cathedral, it?s best to leave your democratic ideals behind." But treating

Hello, I am about to move all of my modelling work into R, and I have been investigating the present state of MCMC and Bayesian methods in R. Following a thread on the mailing list in 2000, I have looked at mcmcpack and Hydra. Three years down the line, is there anything new in this area? I have used both MCSim and WinBUGS in the past. The first one seems promising, but is too focused towards

I am trying to get summary statistics from WinBUGS/JAGS output in the form of mcmc.list objects, using the summary() function. However, I get odd warning messages: Warning messages: 1: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : algorithm did not converge 2: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : algorithm did

I'm trying to read CODA/mcmc files (see the coda package), as generated by jags/WinBUGS/OpenBUGS, into a big.matrix. I can't load the whole mcmc object produced by read.coda() into memory since I'm using a laptop for this analysis (currently I'm unfunded). Right now I'm doing it by creating the filebacked.big.matrix, reading a chunk of data at a time from the chain

On Sat, 28 Dec 2002, Adon Irani wrote: | woah do0d . that's a crackp0t theory if i ever heard one ., | | i'll offer an alternate theory . KiSS - keEp it simple sir ; why should | the icecast2 developers be guilt-tripped into repeating icecast's | performance ? Icecast2 is icecast2. I'm talking about it's lack of innovation (ogg streaming! woopdefuckingdoo) as well as

Greetings all, I'm a newcomer to Bayesian stats, and I'm trying to calculate the Deviance Information Criterion "by hand" from some MCMC output. However, having consulted several sources, I am left confused as to the exact terms to use. The most common formula can be written as DIC = 2*Mean(Deviance over the whole sampled posterior distribution) - Deviance(Mean

Dear Sir/Madam, I am currently doing a project which need to use Bayesian Model. I need to use WINBUGS to do MCMC calculation, and therefore I need to between R and WINBUGS. I tried to submitted install.packages("R2WinBUGS") from R, but the system says, "Warning: unable to access index for repository...". I am not sure what happened with it. Could anyone tell me

Dear list members, Best greetings and apologies for cross-posting. There are available places for the course: "Bayesian Data Analysis with R and WinBUGS", please find the description of the course below. If you have any question don't hesitate to contact me. Best regards, Pablo ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *Three days course in: Bayesian Data

Dear list members, Apologies for cross-posting. Please, find below the information of two statistical courses with R: 1) Statistical Analysis with R 2) Bayesian Data Analysis with R and WinBUGS If you have any question don't hesitate to contact me. Best regards, Pablo ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ *Two days course in: Statistical Analysis with R *Where:

Dear all I used R2WinBUGS package's bugs() function to generate MCMC results. Then I tried to save the simulation draws in R, using read.bugs() function. Here is a simple test: ###################### library(coda) library(R2WinBUGS) #fake some data to test beta0=1 beta1=1.5 beta2=-1 beta3=2 N=200 x1=rnorm(N, mean=0,sd=1) x2=rnorm(N, mean=0,sd=1) x3=rnorm(N, mean=0,sd=1) lambda2= exp(beta0+

Apologies for cross-posting. Places are available in the following course: "Introduction to Bayesian Data Analysis using WinBUGS and R" Best Regards, Pablo ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Course: Introduction to Bayesian Data Analysis using WinBUGS and R Where: Linux Hotel, Essen-Horst, Germany

Dear R Community, I am happy to introduce the latest version 1.3-0 of the 'dclone' R package. The package provides low level functions for implementing maximum likelihood estimating procedures for complex models using data cloning and Bayesian Markov chain Monte Carlo methods with support for JAGS, WinBUGS and OpenBUGS. Data cloning is a global optimization approach and a variant of

Dear R Community, I am happy to introduce the latest version 1.3-0 of the 'dclone' R package. The package provides low level functions for implementing maximum likelihood estimating procedures for complex models using data cloning and Bayesian Markov chain Monte Carlo methods with support for JAGS, WinBUGS and OpenBUGS. Data cloning is a global optimization approach and a variant of

Dear list members, Apologies for cross-posting. Please, find below the information of an introductory Bayesian Data Analysis course with WinBUGS and R. If you have any question don't hesitate to contact me. Best regards, Pablo ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Course: Introduction to Bayesian Data Analysis using WinBUGS and R Where: Linux Hotel, Essen-Horst,

I am trying to use the function 'filter.NA=TRUE' in Geneland. The function appears to be set on TRUE by default, as it appears as TRUE in the 'parameter.txt' file output and hence I do not need to enter the function per se (as it is an 'Unused argument otherwise') . Hence all my missing data (individuals that I have not yet scored at that specific loci) are scored as

Hi all, I would like to compute ergodic mean using MCMC output from WinBUGS. I tried using CODA package, but it seems that it is not implemented yet. Could anyone help me to compute this? Attached to this email are my output and index files. Kind regards, Raquel -- Raquel Rangel de Meireles Guimar?es Doutoranda em Demografia raquel at cedeplar.ufmg.br

A naive question from a non-statistician: I'm looking into running a Bayesian analysis of a model with high dimensionality. It's not a standard model (the likelihood requires a lot of code to implement), and I'm using a Linux machine. Was wondering if someone has any thoughts on what the advantages of OpenBugs are as opposed to just R (or should I be thinking WinBUGS under Wine?)?

Dear List members, I am using R to generate MCMC time series plots. This is the code I use; however, everytime an error message will come out saying could not find function "plot.mcmc." I have coda package and Lattice package installed. Does anyone know what may get wrong here? Thanks so much for your suggestions. Hongli Li *library (coda) dat1=read.csv("Itemtime.csv")

Hello, I have quite a tough problem, which might be able to be solved by MCMC. I am fairly new to MCMC (in the learning process) - so apologize if the answer is totally obvious, and any hints, links etc are greatly appreciated. I'll illustrate the problem in a version cut-down to the essentials - the real problem is ways more complex. Suppose I have a Markovian series of poisson