Exponential Power Model
model
{
for( i in 1 : N )
{
x[i] ~ dexp.power(alpha, lambda)
}
# Prior distributions of the model parameters
alpha ~ dunif(0, 20.0)
lambda~ dunif(0, 1.0)
}
Simulated data set with alpha = 2.5 and lambda = 0.25
The MLEs are alpha = 2.5920853 and lambda = 0.2042697
Data
( click to open )
Inits for chain 1
Inits for chain 2
( click to open )
Results