search

contactmodel / COVID19-Japan-Reff

217 stars 35 forks source link

Definition of CI function in B-MLE #8

Closed taniokay closed 4 years ago

taniokay commented 4 years ago

In calculation of 95% CI, the function CI is like:

CI <- function(par_CI){return(2*(-llk(par_CI)+llk(opt_est$value)))}

But opt_est$value is the value of llk corresponding to the optimal Rt (opt_est$par), so I suppose the function should be like :

CI <- function(par_CI){return(2*(-llk(par_CI)+ opt_est$value))}

aakhmetz commented 4 years ago

Thank you for noticing - it was corrected

taniokay commented 4 years ago

thanks for your replying and correcting!