model
{
for (i in 1 : N) {
x[i] ~ dbern(q) # incidence of HSV
logit(p[i]) <- beta0C + beta * x[i] # logistic model
d[i] ~ dbern(p[i]) # incidence of cancer
x1[i] <- x[i] + 1
d1[i] <- d[i] + 1
w[i] ~ dbern(phi[x1[i], d1[i]]) # incidence of w
}
q ~ dunif(0.0, 1.0) # prior distributions
beta0C ~ dnorm(0.0, 0.00001);
beta ~ dnorm(0.0, 0.00001);
for(j in 1 : 2) {
for(k in 1 : 2){
phi[j, k] ~ dunif(0.0, 1.0)
}
}
# calculate gamma1 = P(x=1|d=0) and gamma2 = P(x=1|d=1)
gamma1 <- 1 / (1 + (1 + exp(beta0C + beta)) / (1 + exp(beta0C)) * (1 - q) / q)
gamma2 <- 1 / (1 + (1 + exp(-beta0C - beta)) / (1 + exp(-beta0C)) * (1 - q) / q)
}