Burr XII model


   [burrxii_ex1]
   
   model
   {
      for( i in 1 : N )
      {
            x[i] ~ dburrXII(alpha, beta)
   
      }
         
         alpha ~ dunif(0.01, 10.0)
         beta~ dunif(0.01, 10.0)   
         diffAlpha <- diffLC(alpha)
         diffBeta <- diffLC(beta)                         
   }

Simulated data set with alpha = 2.0 and beta = 5.0
The MLE’s are using 'maxLik' package in R alpha.mle = 2.11996 beta.mle = 4.03369

Data
list(N=30, x=c(1.3104183, 1.0720485, 0.3748349, 0.6313278, 1.3653709, 0.7758386, 0.5586983,1.0547609, 0.8484189, 1.1714541, 0.6496190, 0.4721445, 0.6450944, 1.1540506, 0.6946741, 0.7307015, 0.7611312, 0.6225298, 1.0603977, 0.2639594, 0.5878855, 1.2542625, 0.6167935, 0.6843502, 0.6144164, 0.9843655, 0.9506997, 1.0561631, 0.7506981, 0.7833595))
Inits for chain 1
list(alpha=1.0, beta=4.0)
   
Inits for chain 2
list(alpha=5.0, beta=8.0)



Results

[burrxii_ex2]

MAP estimates are

[burrxii_ex3]

On this page