model

{

for (i in 1:30){

y

*<- (alpha+1)/(alpha*pow(x*

*,2))*beta*exp(-beta/x**)*(1-exp(-alpha*beta/x**))*

}

alphaI<-1/alpha

betaI<-1/beta

alpha~dgamma(1,0.5)

beta~dgamma(1,0.3)}

alphaI<-1/alpha

betaI<-1/beta

alpha~dgamma(1,0.5)

beta~dgamma(1,0.3)

*}*

list(x=c(12.545023 ,12.284369 ,3.987655, 10.038319 , 8.627547 ,7.946265 ,26.129204 , 7.166539 , 6.628683,5.151729 ,28.152538, 29.142671 ,20.168677 , 5.866670, 28.280793 , 2.937401 ,27.683897, 24.756131, 8.070980, 6.646225 ,21.090805, 16.476565 , 3.456984 , 5.266947, 12.236634, 28.287604, 25.203209, 23.420769, 10.336561, 29.234815))list(x=c(12.545023 ,12.284369 ,3.987655, 10.038319 , 8.627547 ,7.946265 ,26.129204 , 7.166539 , 6.628683,5.151729 ,28.152538, 29.142671 ,20.168677 , 5.866670, 28.280793 , 2.937401 ,27.683897, 24.756131, 8.070980, 6.646225 ,21.090805, 16.476565 , 3.456984 , 5.266947, 12.236634, 28.287604, 25.203209, 23.420769, 10.336561, 29.234815))

*list(alpha=1,beta=1)*