Como posso estimar intervalos de confiança de 95% usando a criação de perfil para parâmetros estimados maximizando uma função de probabilidade de log usando optim em R?
Eu sei que posso estimar assintoticamente a matriz de covariância invertendo o hessiano , mas estou preocupado que meus dados não atendam às suposições necessárias para que este método seja válido. Eu preferiria estimar intervalos de confiança usando outro método.
O método de probabilidade de perfil é apropriado, conforme discutido em Stryhn e Christensen , e no livro MASS de Venables e Ripley, §8.4, pp. 220-221?
Em caso afirmativo, existem pacotes que podem me ajudar a fazer isso no R? Caso contrário, como seria o pseudo-código para esse método?