Skip to content
Prev 385423 / 398506 Next

How to obtain individual log-likelihood value from glm?

John Smith 
Dear R-help,

The function logLik can be used to obtain the maximum log-likelihood value
from a glm object. This is an aggregated value, a summation of individual
log-likelihood values. How do I obtain individual values? In the following
example, I would expect 9 numbers since the response has length 9. I could
write a function to compute the values, but there are lots of
family members in glm, and I am trying not to reinvent wheels. Thanks!

counts <- c(18,17,15,20,10,20,25,13,12)
     outcome <- gl(3,1,9)
     treatment <- gl(3,3)
     data.frame(treatment, outcome, counts) # showing data
     glm.D93 <- glm(counts ~ outcome + treatment, family = poisson())
     (ll <- logLik(glm.D93))

Thread (8 messages)

John Smith How to obtain individual log-likelihood value from glm? Aug 25 Bert Gunter How to obtain individual log-likelihood value from glm? Aug 25 Peter Dalgaard How to obtain individual log-likelihood value from glm? Aug 25 John Smith How to obtain individual log-likelihood value from glm? Aug 27 John Smith How to obtain individual log-likelihood value from glm? Aug 28 Mark Leeds How to obtain individual log-likelihood value from glm? Aug 28 Peter Dalgaard How to obtain individual log-likelihood value from glm? Aug 29 Peter Dalgaard How to obtain individual log-likelihood value from glm? Aug 29