Gompertz model



   [gompertz_ex1]
   
   model
   {
      for( i in 1 : N )
      {
      x[i] ~ dgpz(alpha, theta)      
      }      
   # Prior distributions of the model parameters         
   # Gamma Prior      
         alpha ~ dunif(0.001, 10)
         theta~ dunif(0.00001, 2)               
   }
   
The data set is taken from Murthy et al. (2004, pp. 119 ).

Murthy, D. N. P., Xie, M., Jiang, R. (2004). Weibull Models, Wiley-Interscience.

MLE's are alpha = 1.968482 theta = 0.019104

Data
list(N=20, x=c(0.481, 1.196, 1.438, 1.797, 1.811, 1.831, 1.885, 2.104, 2.133, 2.144,
2.282, 2.322, 2.334, 2.341, 2.428, 2.447, 2.511, 2.593, 2.715, 3.218))
Inits for chain 1
list(alpha=1.0, theta=0.1)
   
Inits for chain 2
list(alpha=5.0, theta=1.2)



Results



[gompertz_ex2]

On this page