Dear All,
I am running Bayesian meta-analysis using the code below. I will like to calculate a proability that drug A and drug B is equivalence (log odds ratios is between -1 and 1). I will be glad if some can show me how to use 'and' function in WinBUGS
pr1 <- step(d-1)
pr2 <- step(1-d)
model
{
for( i in 1 : Num ) {
rc[i] ~ dbin(pc[i], nc[i])
rt[i] ~ dbin(pt[i], nt[i])
logit(pc[i]) <- mu[i]
logit(pt[i]) <- mu[i] + delta[i]
mu[i] ~ dnorm(0.0,1.0E-5)
delta[i] ~ dnorm(d, tau)
}
d ~ dnorm(0.0,1.0E-6)
# Choice of priors for random effects variance
#tau ~ dgamma(0.001,0.001)
#sigma <- 1 / sqrt(tau)
tau<-1/(sigma*sigma)
sigma~dunif(0,10)
delta.new ~ dnorm(d, tau)
}
Regards
Frank Peter
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail [log in to unmask]
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
|