Uma vez que você tenha a verossimilhança pode usar também a função mle2() do pacote bbmle, que usa também a optim() mas arruma a saida e facilita para a construção de intervalos de confiança, via aproximação quadrática e também perfis de Verossimilhança.<br>
<br><br><br><br>-- <br>Wagner Hugo Bonat<br>LEG - Laboratório de Estatística e Geoinformação<br>UFPR - Universidade Federal do Paraná<br>