Print

---

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