En R, dada una salida de optim con una matriz de arpillera, ¿cómo calcular los intervalos de confianza de los parámetros utilizando la matriz de arpillera?

Question

Dada una salida de optim con una matriz de arpillera, ¿cómo calcular los intervalos de confianza de los parámetros utilizando la matriz de arpillera?

 fit<-optim(..., hessian=T)

hessian<-fit$hessian

Lo que más me interesa es el contexto del análisis de máxima verosimilitud, pero tengo curiosidad por saber si el método se puede ampliar más allá.