Questions about BUGS language for defining Bayesian models that can be run by the WinBUGS software package. It is advised to supplement your question with a reproducible BUGS example (http://stackoverflow.com/q/5963269).

learn more… | top users | synonyms

0
votes
0answers
5 views

Wishart distribution in BUGS with prior on the scale matrix

I am attempting to sample from a wishart in WinBUGS where the scale matrix has priors on the entries. dwish() in WinBUGS won't allow this, so I am attempting to ...
0
votes
0answers
19 views

Confidence intervals of proportions for clustered data

I am trying to compute the proportion of "yes"'s in data that is clustered (participants within cities). If I ignore the cities, I get very nice and narrow confidence intervals. If I compute the ...
0
votes
1answer
31 views

Undefined real result error at WinBUGS

I am currently working on my thesis and interested in estimating a multilevel differential item functioning model and I using at WinBUGS. Until I had done model check-up, there are no errors. However, ...
0
votes
0answers
11 views

Measurements from two raters. Should I use multilevel/“random-effects” model?

I have several variables that were measured on patients, for most patients the variables of interest were measured by 2 independent "raters". Most of the variables are binary. I need to compute the ...
3
votes
0answers
40 views

Is my OpenBUGS / WinBUGS model well specified?

I've just started trying to use OpenBUGS for Bayesian analysis of stochastic volatility models. In particular, I'm trying to calculate stochastic covariance, similar to the DC-MSV model specified by ...
0
votes
0answers
15 views

Simulation m values from posterior distribution with WinBUGS

I try to simulate from posterior distribution with WinBUGS. My data came from Multinomial distribution, i.e. : $y_i~\text{Multi}(n;p_1,p_2,p_3)$. A common prior for multinomial is Dirichlet, i.e.: ...
0
votes
0answers
18 views

Simulation from Dirichlet distribution with WinBUGS

I have a question. Now I am learning WinBUGS, doing bayesian statistics. How, can I simulate a Dirichlet distribution (which is the posterior, for my model ...
1
vote
0answers
42 views

What is the form of link function in this BUGS multinomial regression example?

The alligators example from openbugs examples repository is the same example that comes with winbugs. Basically this is a multinomial logistic regression example in which the outcome variable has 5 ...
7
votes
2answers
225 views

Bayesian estimation of $N$ of a binomial distribution

This question is a technical follow-up of this question. I have trouble understanding and replicating the model presented in Raftery (1988): Inference for the binomial $N$ parameter: a hierarchical ...
2
votes
0answers
55 views

Dawid Skene model for BUGS

I'm a beginner in BUGS. I'm trying to code the Dawid Skene model in BUGS. The model is as follows: Currently, I'm adapting the code from Stan, here. Here is the code: ...
1
vote
0answers
20 views

Longitudinal data format in WinBUGS or OpenBUGS, different number of items per time point

I am working with longitudinal data, where there are a different number of "items" per time point. I am working with the MCMCdynamicIRT1d function in ...
1
vote
0answers
15 views

WinBUGS multinomial all parameters are random variables

everyone here. I want to use multinomial distribution, but p[] and N are all random. model is correct and data loaded , model compilied, but it can't still run , it says:value of binomial f[1] must be ...
0
votes
0answers
10 views

WinBUGS multinomial with unknown parameter and size

I want to use multinomial distribution, but i don't know the p[] and N. So p[] and ...
1
vote
0answers
35 views

Zero trick in WinBUGS for truncation

I am using the Zero-trick to create a prior for correlation parameter. Data y is a bi-variate normal with covariance matrix with same diagonal elements and ...
1
vote
0answers
25 views

How are deterministic/logical nodes represented mathematically in directed graphical models?

Various software for performing inference on graphical models support logical(a.k.a deterministic) nodes. PyMc, Winbugs, Smile are the ones I'm aware of. Based on the different inference methods ...
1
vote
0answers
46 views

Multigroup Longitudinal Item Response WinBUGS OpenBUGS

I was able to fit the Longitudinal IRT model in Winbugs for an ordinal response by extending the BUGS code I took from the paper by Curtis in JSS http://www.jstatsoft.org/v36/c01/paper/ However, I am ...
0
votes
0answers
22 views

Bayes Hierarchical model - sampling with additional constraints

I am trying to create a bayesian model for the following problem: One has a barrel that contains many white balls and a few red balls. One has an uninformed prior on the concentration of red balls ...
1
vote
0answers
13 views

WinBUGS sript: display correlation [closed]

Does anybody know if it is possible to display the correlation matrix using WinBUGS script? (with the GUI it can be found in th Inference menu -> correlations -> print) Thank you everybody!
2
votes
0answers
54 views

Dummy coding from a bayesian stand of view

If we want to use categorical variables in regression context we are allowed to use dummy codings such as: http://www.ats.ucla.edu/stat/r/library/contrast_coding.htm Is this also required in ...
2
votes
0answers
39 views

Truncated normal in WinBUGS

I want to use the following truncated normal distribution in WinBUGS to estimate parameters of SEM using Bayesian analysis. ...
1
vote
0answers
34 views

Inverse normal distribution in winbugs

I want to use the inverse normal distribution in winbugs but I don't find it there. Please help me to find a solution. I have here truncated normal distribution for y ...
3
votes
2answers
107 views

Joint prior distributions in WinBUGS

Suppose we have a hierarchical model summarised by the following: $y_{i} \sim N(\mu_{i}, \sigma^{2})$, for $i = 1, \ldots, n$; (For these purposes, assume $\sigma^{2}$ is known) where $\mu_{i} = ...
2
votes
0answers
41 views

How to make WinBUGS close when “Trap” occurs [closed]

I'm using R2WinBUGS to process several datasets in WinBUGS in "batch mode". Occasionally traps occur, mainly due to bad parameters (I guess), and I managed to avoid ...
2
votes
0answers
273 views

WinBUGS/OpenBUGS Gamma distribution parameters [closed]

According to the WinBUGS manual the gamma distribution is defined by: dgamma(r,mu) However, what is r? Is it the shape, scale ...
1
vote
1answer
89 views

Specification of logical node (with distribution?) in WINBUGS

For a piece of homework I have an assignment using WINBUGS which I must admit confuses me to say the least. Tangential to my question but I have a few stochastic nodes that are to be gamma ...
1
vote
1answer
101 views

JAGS burn-in phase takes ZERO time? [duplicate]

My breath was taken away by this. I tried to run my JAGS model with different number burn-in samples but it still takes the same amout of time!! ...
4
votes
0answers
265 views

CAR spatial models in JAGS

WinBUGS comes with the GeoBUGS add-on, which contains a number of predefined model structures that are suitable for modelling spatial data structures e.g. geostatical structures (spatial.exp), ...
0
votes
1answer
182 views

Undefined real result in a zero-inflated negative binomial

I have run a zero-inflated Poisson model in WinBUGS without problems, and now I am trying to run its equivalent negative binomial. However, I get an "undefined real result" trap message over and over. ...
1
vote
1answer
553 views

Estimation of the shape parameter of gamma distribution in Winbugs

I want to estimate the shape parameter of gamma distribution in Winbugs. I select gamma distribution as prior for shape parameter. Data set is generated using MATLAB as: ...
0
votes
0answers
123 views

How exactly is car.normal model specified?

Looking at the GeoBUGS manual version 12, there is a confusion how exactly the car.normal model is specified: 1) at page 12, I read car.normal and car.l1 [top] ...
3
votes
1answer
365 views

Define Priors for Dirichlet Distribution parameters in JAGS

I'm defining a Multinomial-Dirichlet model in JAGS and want to assign some hyperpriors to the parameters of the Dirichlet distribution. In the WinBugs manual I read that "the parameters of Dirichlet ...
1
vote
1answer
159 views

How can I estimate the precision of a normal using a Gibbs sampler?

I am trying to estimate the precision $\tau$ of a normal distribution with either WinBUGS or OpenBUGS: $c \sim \text{normal}(\mu,\tau)$ $\mu \rightarrow \lambda \cdot t^{-\beta}$ $\tau \sim ...
1
vote
0answers
125 views

The effective sample size from WinBUGS results

I am running a Bayesian regression model by WinGUBS via R2WinBUGS package in R. Everything looks fine except for one parameter: ...
1
vote
1answer
460 views

Trap 66 in WinBUGS in a hierarchical Bayesian modeling

I want to analyze a multilevel multidimensional model in WinBUGS. the model is as below (N=2362 students responding to K=45 items of a test, students are nested within J=116 schools): ...
1
vote
1answer
437 views

WinBUGS: Multiple definitions of a node

So this question is about the BUGS modeling language. So you either know it or have no clue. I'm a newbie to this so it's been driving me mad. I want to define a simple two-state hidden Markov model ...
2
votes
0answers
117 views

How to define a BUGS/JAGS or WinBUGS model with a likelihood node that includes two integrals?

I would like to implement a BUGS/WINBUGS or JAGS model that specifies the likelihood node by marginalizing over two parameters. Is this feasible in BUGS/WINBUGS or JAGS? This is the posterior ...
2
votes
3answers
152 views

Bayesian Vs MLE regression - getting different results

I've set up a Bayesian regression model in WinBUGS to determine values for the unknown parameters (b1, b2, b3, b4) and intercept value (b0) in a linear regression model. The code is as follows: ...
0
votes
1answer
64 views

Looking for step by step example of sampling from DAG in Bayesian model

I am looking for a tutorial type example that shows the step by step process sampling from a simple hierarchical model. For example, I am trying to study the distribution of ...
1
vote
0answers
93 views

How can I compute covariances between parameters in WinBUGS?

I am wondering if there is an easy way to compute the covariances between parameters in WinBUGS/OpenBUGS. It is easy to obtain the variances, but for subsequent analysis, I need the covariances. ...
2
votes
1answer
536 views

WinBUGs-TIme Series Code

I am looking to get started using time series (arima) (regression) (etc) evaluating log stock returns in WinBUGs. Does anyone have a simple code to get me jump started?
1
vote
1answer
94 views

Multiplying a matrix by a scalar which has a prior distribution in OpenBUGS

So I am having a problem specifying my model in OpenBUGS. A set of vectors in a linear regression model is given a multivariate normal prior with a constant mean vector and a constant precision matrix ...
3
votes
0answers
143 views

How to use an initial posterior for recursive / sequential updating in WinBUGS

I am using WinBUGS to estimate / update the parameters of a model. The model is: $$ \begin{aligned} D(T,B,a)&= B*(a_0+a_1T+a_2T^2+a_3T^3)+error(B,T,a) \\ error &= \mathcal N(0, ...
2
votes
1answer
312 views

Logistic Regression - Bayesian Approach - Assessing Classification Precision

I have recently begun to read about bayesian statistics and I am playing around with the R2WinBUGS package. I'm trying to fit a logistic regression to the spam data (that can be found on the webpage ...
1
vote
2answers
481 views

Thinning chains in BUGS/JAGS

Hi I have a quick question about the details of running a model in JAGS and BUGS. Say I run a model with n.burnin=5000, n.iter=5000 and thin=2. Does this mean that the program will run 5,000 ...
0
votes
1answer
653 views

WinBUGS truncated normal distribution [duplicate]

I am estimating a stochastic frontier with a mixed model. So far the half normal distribution worked good but I need a truncated normal distribution. It does not work, and I receive the error ...
1
vote
0answers
800 views

How do I identify this trap error in WinBUGS?

I am currently working on my thesis and interested in estimating a Weibull conditional Hazard Frailty Model for recurrent event data using WinBUGS. I wrote the model but when am trying to use it with ...
3
votes
0answers
251 views

Calculate DIC from WinBUGS output

Has anyone ever tried to calculate DIC "manually" from WinBUGS output? I tried for the simplest linear model, and could get exact the same Dbar, but always a bit different Dhat. I think Dhat ...
4
votes
2answers
674 views

Variables that do not converge in winbugs

I conducted a Bayesian analysis in winbugs and then checked the convergence from the history plots. The regression coefficients look stabilized, but the variance parameters don't. I got the number of ...
1
vote
0answers
142 views

How do you code three summation functions in WinBUGS?

I recently learned how to use winbugs. I know nothing about coding and programming. My model is not included in the set of distributions that are already available in winbugs. This is a trivariate ...