MODEL{
a~dgamma(r,mu)
c~dgamma(d,e)
x~dpar(alpha,c)}
data list (
n=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39)
x=c(0.4998870.5004730.5040090.5213740.5264050.5346220.5553860.5596610.571820.6030740.6223970.6342640.6886210.7200430.7227180.7334670.7402360.7421210.7536540.7652220.7921360.7962290.8173770.8585270.8669270.874310.9064660.9165620.9416860.9596180.9638710.9705111.0095151.2167781.3134761.7793431.8661931.9257272.504857)