logLik.plm function

logLik.plm.R

logLik.plm <- function(object){
  out <- -plm::nobs(object) * log(2 * var(object$residuals) * pi)/2 - deviance(object)/(2 * var(object$residuals))
  
  attr(out,"df") <- nobs(object) - object$df.residual
  attr(out,"nobs") <- plm::nobs(object)
  return(out)
}

Shouldn't plm::nobs(summary(object)) read plm::nobs(object)?

Spot on, thanks @Helix123!
I'll change it! It should, of course, give the same result!

This is great! Thanks for coming up with this!

thank you!