[R-sig-ME] Extract individual likelihood from glmmTMB
Ben Bolker
bbo|ker @end|ng |rom gm@||@com
Thu May 9 19:33:17 CEST 2019
By "individual likelihood" do you mean you want the log-likelihood of
a particular model? logLik(fitted_model) should do it. (This is a
generic method that should work for most model types in R, although you
do have to be a bit careful in comparing log-likelihood values across
models from different packages as they may make different choices about
which constant terms to include/exclude.)
If you want something else, please clarify ...
cheers
Ben Bolker
On 2019-05-09 4:14 a.m., Suna Park wrote:
> Dear R-sig-mixed-models members,
>
> Hello, this is Suna.
>
> I'm studying the use of medical care for the elderly using zero-inflated
> Poisson and negative binomial mixed model with random intercept.
>
> And I'm considering how to extract the individual likelihood from the
> glmmTMB for Vuong(1989)'s test.
>
> Looking at the source code of the glmmTMB, it seems that the likelihood
> function is in the 'glmmTMB.dll'.
>
> But I can't open and see the function in the 'glmmTMB.dll'.
>
> Could you let me know about how to modify the likelihood function or
> extract individual likelihood?
>
> Thank you,
>
> Suna Park
>
>
> *Suna Park*
> College of Economics & Finance, Hanyang University,
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
More information about the R-sig-mixed-models
mailing list