model{
for (i in 1:n)
{ y[i] ~ dpois(lambda[i])
log(lambda[i]) <-  b[1] + b[2]*log(i)
}
b[1] ~ dflat()
b[2] ~ dflat()
}

list(n=20, 
y=c(1,6,16,23,27,39,31,30,43,51,63,70,88,97,91,104,110,113,149,159) )

list(b=c(1,1))

model{
for (i in 1:n)
{ y[i] ~ dpois(lambda[i])
log(lambda[i]) <-  theta*log(i)
}
theta ~ dflat()
# theta ~ dnorm(0,0.001)
}

list(n=20, 
y=c(1,6,16,23,27,39,31,30,43,51,63,70,88,97,91,104,110,113,149,159) )
list(theta=0.5)