Olá Sérgio, Altere o inits para inits <- function(){list(p=rbeta(N,1.001,1.001),r=rgamma(N,0.01,0.01))} Abraço Márcio Nicolau Send from mobile Em 11/05/2012, às 11:55, Sérgio Henrique almeida da silva ju<sergio.edfisica em gmail.com> escreveu: > inits <- function(){list(p=rbeta(N,1.001,1.001),r=rgamma(1,0.01,0.01))}