**Chapter 2 Exercises**

**Coin tossing**

**Solutions**

1. Write your own BUGS code for calculating the probability of obtaining at least 15 heads in 20 tosses of a fair coin.

model {

Y ~ dbin(0.5, 20)

P <- step(Y - 14.5)

}

` node mean sd MC error 2.5% median 97.5% start sample`

P 0.021 0.1434 0.00158 0.0 0.0 0.0 1 10000

Y 10.02 2.244 0.02181 6.0 10.0 14.0 1 10000

2. Adapt your code to calculate the probability that no more than 6 patients in a group of 30 will fail to respond to a drug with a success probability of 0.7.

Modifications to the code are shown in red below. We first change the success probability from 0.5 (fair coin) to 0.7. Requiring no more than 6 failures is the same as requiring at least 24 successes. Alternatively we could calculate the number of failures, numfail , as 30 – Y, and then compute the required probability directly, as shown below ( P.alt ).

model {

Y ~ dbin( 0.7 , 30 )

P <- step(Y - 23.5 )

numfail <- 30 - Y

P.alt <- step(6.5 - numfail)

}

` node mean sd MC error 2.5% median 97.5% start sample`

P 0.1616 0.3681 0.003231 0.0 0.0 1.0 1 10000

P.alt 0.1616 0.3681 0.003231 0.0 0.0 1.0 1 10000

Y 21.01 2.517 0.0246 16.0 21.0 26.0 1 10000

numfail 8.992 2.517 0.0246 4.0 9.0 14.0 1 10000